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1 Introduction 

In aerospace applications, magnetic suspension systems may be required to operate over large variations 
in air-gap. Thus the nonlinearities inherent in most types of suspensions will have a significant effect. 
Specifically, large variations in operating point may make it difficult to design a linear controller which 
gives sat isfactory stability and performance over a large range of operating points. 

One way to address this problem is through the use of nonlinear compensation techniques such as 
feedback linearization. In the last decade, this area has seen a great deal of research effort among system 
theorists [6], [7], [8], and has been successfully used in flight control applications [9], but has yet to be 
widely used in the area of magnetic suspensions. 

Nonlinear compensators have received limited attention in the magnetic suspension literature. In 
[11] t he control system design for suspension of a 1-ton, 4-passenger vehicle is presented. A flux sensor 
is located on the pole-face of the suspension magnet. A minor feedback loop is closed on flux which 
linearizes the dependence of magnetic force on position. The force goes as the square of flux, so an 
additional square root linearization is required. The idea of using an analog multiplier to compute the 
ratio i/x is mentioned, which would allow elimination of the flux sensor However, this technique is 
reported as prone to drift and noise, and was abandonded. This flux- feed hack scheme is described in 
more detail in [12]. In [13] a nonlinear correction law is used to correct the inverse square law magi et 
behavior in a flywheel suspension. The nonlinear compensation is implemented with analog multiplier 
and square root circuits. Microprocessor-based linearizing transformations are reported in [14] in t he 
context of a demonstration system. 

In recent years, progress has been made in the theory of nonlinear control systems, and in the sub- 
area of feedback linearization. Here, [6] is of fundamental importance in that it presents the conditions 
under which a system may be linearized. In [7] these results are globally extended, and in [8] the theory 
is developed for multi-input /multi-output systems. In a subsequent section the results of [C] are applied 
to a third order suspension example. However, for a simple system, it is often possible to construct 
the linearizing transformations by inspection. We start then by demonstrating the idea of feedback 
linearization using a second order suspension system. In the context of the second order suspension, 
sampling rate issues in the implementation of feedback linearization are examined through simulation. 
The system which is studied is presented in the next section. 


2 Nonlinear Suspension Model 


In this section, the open-loop dynamics for a simple one-degree-of-freedom suspension ere presented, his 
system exhibits the essential issues faced in the design of tractive type suspensions, that is, suspem ms 
which operate as variable reluctance devices The example is drawn from [10] pgs. 22 23, 84-86, md 
193-200. The only change is that the system is inverted such that gravity acts to open the air-gap. 1 Ids 


system is shown in Figure 1. 

The details of the electromagnetics are worked out in [10]; for th 
details are the coil voltage 
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and the force on the plunger 
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whore the first term is the electromagnet force, the second term is the gravitational force on the plunger, 
and the third term is a disturbance force acting on the plunger in the direction the electromagnet force. 
The nonmagnetic sleeve is assumed to exert no frictional forces on the plunger. 

If we define 


u = v c 
C = wdtioN 2 

X[ = X 4 00, (*i > 9o) 

x\ - x 2 

then the state equations for the open-loop suspension are 
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2.1 Second order system 

If the coil current is assumed to he the control input, then the suspension state equations are reduced 
to second order. 


= x 2 



These equations will adequately model the system if the coil current is controlled by a bigh-bandwidth 
current loop with sufficiently high voltage-drive capabilities. In applications, it is most typical to drive the 
coil with such a current loop, as this essentially eliminates the dependence of position-loop performance 
upon the electromagnet coil resistance and inductance. 


3 Linearization of second-order suspension 

The basic idea of feedback linearization is to define transformations on the states and input(s) such 
that the nonlinear system appears linear and operating-point invariant in terms of I lie transformed 
representation. Then a controller can be designed for the transformed variables. 1 Ins allows the closed- 
loop system stability to be made independent of operating point. 

For the second-order equations (4), a transformation on the input is all that is required to linearize 
the system. This transformation may be derived by inspection without using any formal mathematical 
machinery; this is the approach taken in [11], [12], [13], and [14]. That is, if the coil current i is made to 
vary as 

i = x '\f^r (5) 

then the suspension is globally linearized in terms of the new input v. 'The notation for the auxiliary 
input, v has been chosen to match Ihe notation in [(>]. 

Specifically, substituting from (5) into (4), the system stale equations become 

X} - X'l 
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These equations are linear, with an input v } and disturbance terms g and fd • 

Here, v is a signal internal to the compensator which may be thought of as a setpoint for acceleration 
in the direction of increasing airgap. In operation, the signal v will be computed within the compensator, 
and constrained to be less than or equal to zero. Since the magnet can only supply accelerations in the 
direction of decreasing air gap it would not be physically meaningful to ask for acceleration in the 
direction of increasing air gap by setting v greater than zero. Thus the term — v in (5) will always be 
greater than or equal to zero, and the square root will yield a real number. 

The plant appears linear in terms of the new input v. This compensation of the nonlinear term does 
not however stabilize the plant. To stabilize the system, the nonlinear compensator is preceded by a 
linear compensator. The resulting closed-loop system is shown in Figure 2. The compensator may be 
thought of as having two parts, a nonlinear compensation section and a linear compensation section. It is 
the function of the nonlinear section to implement (5) in order to adjust t as a function of xi and v } such 
that the acceleration of the plunger is equal to v . It is the function of the linear section to specify the 
value of v as a function of the error between the position setpoint and the measured position such that 
the linearized plant is robustly stabilized and has good disturbance rejection and settling time properties. 
The signal v forms the connection between the linear and nonlinear sections of the compensator. 

This combination of linear and nonlinear compensation sections stabilizes the plant such that the loop 
dynamics are independent of operating point. Such operating point independence is the main advantage 
of using a nonlinear compensator. Note that as viewed from the input to the nonlinear section, the 
incremental relationship between v and x\ is equal to 1/s 2 , independent of operating point. I hus, the 
linear compensator can be designed to control a double integrator via standard linear techniques. If it is 
desirable to reject static disturbance forces with no position error, then the linear compensator can be 
designed to include an integral term. This integral term will adjust the value of v to balance gravity and 
any low-frequency components of the disturbance fd • 

In applications where large excursions or disturbance forces arc anticipated, the additional complexity 
of the nonlinear compensation approach is justified. The major caveat is that we are assuming that the 
suspension model is accurate. For the electromagnetics an accurate model can readily be developed, and 
thus nonlinear compensation techniques are applicable. The nonlinear compensation technique was used 
in the construction of a class demonstration system which is described below. 


4 Classroom demonstration implementing linearization 

In the Spring of 1988, the author constructed a single degree of freedom levitation system for use as 
a classroom demonstration which implemented the nonlinear compensation technique described in the 
previous section. As developed there, if the plant state equations are given by (4), then applying the 
nonlinear compensation law (5) results in a system which appears to be linear in terms of the intermediate 
signal v. The demonstration system uses a high-bandwidth current-drive to regulate the electromagnet 
current, and thus (4) is applicable. 

In the demonstration system, a one inch steel ball bearing is suspended below an electromagnet 
consisting of 3100 turns of #22 magnet wire wound on an 1 inch diameter by 4 inch length steel core. 
The coil current is controlled by a Bose-type switching regulator, with a half-scale current switching 
frequency of 10 kHz, and a full scale current of 2 Amperes. The operating point current is about 0.4 
Amperes at a typical operating point air gap of 1 cm. The system is digitally controlled by an 8088/8087- 
bascd single-board computer and data acquisition system at a 400 Hz sampling rate. The control law 
for the nonlinear compensation section uses (5) to linearize the magnetic force relationship. 4 his allows 
the stability of the closed-loop system to be essentially independent of the operating point. The control 
law for the linear compensation section is then developed via classical techniques applied in the discrete- 
time domain. The position of the ball is sensed optically, and nonlinearities in the sensor output versus 
position are compensated for in software. 

In order to apply the nonlinear compensation technique, an accurate model of the plant is required. 
For the classroom demonstration, this model is developed by measuring the force on the ball as a function 
of current and position. This measurement is accomplished by using a balance beam for measuring the 
magnetic force on the ball. A 1 inch ball bearing is glued into one end of an aluminum balance beam 
of rectangular tubular cross-section measuring 1 inch wide by 0.75 inches deep by 12 inches long. The 











beam pivots at the center on a thin wire which is held by fixed side supports. Balance pans are hung 
from the beam on both sides of the pivot midway between the pivot and the ends of the beam. T ese 
pans are used to add or subtract weight carried by the suspension. At the end of the beam opposite the 
ball, a micrometer is positioned to push against the beam and thus provide a position reference. 

The idea here is that the ball glued into the end of the beam can be placed into suspension and the 
beam thus provides a handle by which the force applied to the ball can be varied. This is accomplished 
by putting weights into the balance pans on either side of the pivot. The beam is made of aluminum, an 
thus does not interact with the electromagnet at low frequencies. At the ball end, the beam is made thin 
so as not to interfere with the optical measurement. The ball is attached to the bottom of the thinned 
end of the beam such that it interacts with the optical sensor in the same fashion as a freely suspended 

The force relationship (2) was well fit by the experimental data with the parameters C = 4.43 x 
lO -4 Nm 2 /A 2 and g 0 = 0.25crn. The mass of the ball is 67 grams. These parameters are used in the 
nonlinear compensation law (5). The only deviation from the relationship (2) is at high currents (> 1A), 
where the effects of magnetic saturation are apparent. 

The optical position sensor is constructed as follows. A 24 volt, 5 watt incandescent lamp is used as 
the source, and a piece of cadmium sulfide photo-cell is used as the sensor, in what is a standard position 
sensor for magnetic suspensions. Using the balance beam described above, the sensor output is measured 
for a number of ball positions. When the shadow-line cast on the sensor is in the central region of the 
sensor, the sensor output is essentially linear with ball position. However, as the shadow- line approaches 
the upper or lower edge of the photo-cell, the sensor sensitivity begins to decrease. This nonlinearity in 
the relation between ball position and sensor output is corrected in software in the section of code which 
inputs the sensor voltage. The corrected position measurement is then linear with actual ball position. 
It is this corrected position measurement which is passed to the rest of the control loop. 

The position sensor was found to have several defects which limit the system performance. Fits , 
the incandescent bulb output decreases significantly as a function of time. This is believed to be due 
to the evaporation of the filament. Material driven off of the filament is deposited on the inside of the 
glass envelope, thereby decreasing the bulb brightness. The second problem is that the cadmium sulfi Ic 
sensor is sensitive to any light falling on its surface, independent of the source. 1 hus ambient hg iting is 

indistinguishable from the light emitted by the bulb. 

Both of these effects cause problems in the nonlinear compensation law (5) and in the correction 
of the sensor nonlinearities. First, the decrease in bulb intensity and any changes in average ambient 
light act as offset terms which drive the system to incorrect points on the sensor correction curve an 
in the magnet nonlinearity correction law (5). This offset deteriorates the system stability. Secondly, 
the ambient ligty has a large component at twice the power line frequency, especially in rooms with 
fluorescent lighting. This signal at 120 Hz acts as a large noise source which causes error motions in the 
ball position. 

The above problems can be solved as follows. First, the light source needs to be made more constant 
with time. This can be achieved by using a more specialized incandescent bulb, or by switching to a 
semiconductor light source such as an infra-red light emitting diode. The ambient lighting offset and 
noise problems can be solved by either or both of two approaches which are classical. The first is to make 
the system narrow-band. Commonly available IR diodes emit a relatively narrow-band optical s.gnal; 
laser diodes are narrower. In this case, an optical hand-pass filter can be placed in front of the sensor, so 
that only the emitted frequencies are sensed, and the ambient lighting is greatly attenuated. The second 
approach is to switch the light-source on and off at a high frequency and use synchronous detection to 
reject signals which are not at the same frequency and phase as the source. The frequency of switching 
must be made much higher than the cross-over frequency of the position control loop, perhaps on the 
order of 10 kHz switching frequency. This rate is easily within the capabilities of available electronics. 

The results derived in the previous section for the nonlinear compensation laws assume that these are 
implemented in continuous time. For discrete-time implementation, the issue of sampling rate becomes 
important. This issue is investigated in the next section. 
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4.1 Sampling rate issues 

Duo to the complexity of the transformations it is most likely that a linearizing compensator will be 
implemented in discrete time. As an introduction to one issue involved in discrete-time implementation, 
the effect of sampling rate on the second-order suspension system (4) is investigated by simulation. For 
this example, the suspension parameters have been been given the values developed for the class demon- 
stration system described above. These values are M = 67 grams, and C = 4.43 x 10“ 4 Nm 2 /Amp 2 . 

The system is simulated assuming a nonlinear compensation law of the form (5). The four graphs 
shown in Figure 3 indicate the system behavior when a net 0.05g acceleration ( v = —0.05 in (5) ) is 
specified. The lines labelled ‘ideal’ show that if the nonlinear compensation was perfectly implemented, 
the force on the ball would be constant, and the graph of velocity vs. time would be a straight line. 
However, with any finite sampling rate this is not the case. The system is open-loop unstable, and 
uncontrolled between sampling instants. Thus it ‘runs away’ during the interval in which the control 
current is held constant. The graphs show the result of this process for sampling rates of 1 kHz and 200 
Hz. To get reasonable behavior, it can be seen that a sampling rate on the order of 1 kHz is required. 
In the class demo, due to computational speed limitations a 400 Hz sampling rate is used. This is found 
to be adequate as long as the ball is not allowed to approach too close to the pole face. 

Another way to look at the effect of sampling rate is to examine the system behavior under closed-loop 
position control. To this end, a linear proportional plus lead compensator is designed in discrete-time to 
stabilize the nominal plant which would result if the nonlinear compensation were perfect. That is, in the 
ideal case, the nonlinear compensated system appears as a double integrator independent of operating 
point. In the finite sampling time implementation, the quality of this approximation deteriorates as the 
air gap closes. This can be seen in Figure 4 which displays simulated step responses for the closed-loop 
system at four nominal operting points and for the two sampling rates. Note that the system with 200 Hz 
sampling goes unstable at the 0.5 cm and smaller air gaps, whereas the behavior of the 1000 Hz sampled 
system only begins to deteriorate when the air gap approachs 0.3 cm. The unstable response for 200 Hz 
sampling is not shown for the 0.3 cm air gap. 

This example shows that the practical implementation of linearizing transformations may require 
very high sampling rates. Also, what may be considered a satisfactory sampling rate depends on the 
range of operating points which are encountered in system operation. Certainly, the issue of discrete-time 
implementation merits further study. 

Experience with this simple nonlinear compensation system provided the impetus toward an under- 
standing of feedback linearization techniques in more generality. A description of the application of 
feedback linearization to the third order suspension system is given in the next section. 

4.2 Linearization of third-order suspension 

For more complex plants it may be difficult to develop linearizing transformations by inspection. The 
results of [6] provide a general approach to this problem. A good introduction to these ideas is presented 
in [4] and [5]. These references assume no more than an undergraduate background in control theory 
and are thus a good place to start for someone new to this area. 

Without reviewing the results from the above references, if the plant satisfies a controllability con- 
dition and a condition on the existence of solutions to a set of partial differential equations, then trans- 
formations z\ — Ti(x) t .,.,z n = T n (x),v — 7' n +\(x, u) can be constructed such that in the z-v space 
the system appears linear. Here, x is the state-vector of the nonlinear system, 2 is the state-vector of 
the linearized system, and n is the system order. Under these transformations, the nonlinear system is 
mapped to the controllability canonical form 



For the system (3), the required conditions are satisfied, and the results of [6] yield the linearizing 
transformations z\ — xj, z 2 — x 2) z 3 = — (C/M)(i/x[) 2 } and z 4 = — u). Thus the system 
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Figure 3: The open-loop system with nonlinear compensation showing its performance with sampling 
periods h of 1 and 5 milliseconds. Position x is in cm separation from the pole face, velocity t; is in 
cm/sec, force on the ball is in dynes, and current I is in amps. 
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Figure 4: The closed-loop system with nonlinear compensation showing its performance with sampling 
periods h of 1 and 5 milliseconds as the operating point position is moved towards the pole face. Position 
x is in cm separation from the pole face. The nominal operating points are 2, 1, 0.5, and 0.3 cm 
respectively. 



appears linear in terms of state varibles 2l , z 2 , and z 3 , and with a properly redefined input v The 
states z, and z 2 are simply the original position and velocity. State z 3 is the acceleration applied to 
the suspended member. Thus it makes physical sense that the suspension will appear linear in z 3 . 1 e 
suspension force happens to vary nonlinearly with the untransformed state and input, but Newton s law 
guarantees linearity in terms of a transformed state varible which is proportional to acceleration. In an 
implementation, the voltage drive u must be computed in terms of v: 

u = _Mx } v +iR (8) 

i 

Since v drives the derivative of z 3 , we can think of v as being a setpoint for the slope of the acceleration. 
Note that the coil resistance voltage drop iR is directly added to the input u. 

Thus we have found a set of linearizing transformations. However the transformations are not unique. 
Direct substitution will verify that the transformations Ti=x\,T 2 - 2x|X 2 , 


T 3 = 1x\ - 


2C’r 

Mx\ 


( 9 ) 


and 


6Cr 2 J 2 2t v 

Ta = — rr-r + - «) 


Mx i 


M' 


( 10 ) 


though more complex than the first set., do indeed globally linearize the system. Actually, there are 
an infinity of such transformations which linearize this system. This is a consequence of the natur • ol 
nonlinear systems. It is clear however that the first set has the greatest physical meaning, and thus 
would be chosen in any practical context. Note also that in the first set the transformed state z 3 need 
never be computed. This is so because the system is linear between the transformed input v and the 
original position state variable x,. Further, note that the input transformation (8) depends only upon 
position xi and current i. Both of these quantities may be readily measured. 


5 Conclusions 

As we have seen in the magnetic suspension examples, the technique of feedback linearization is of 
great utility in designing control loops for nonlinear systems such that the closed-loop systems are well- 
behaved despite large variations in operating point or disturbance forces. Sampling rates for discrete- ime 
implementations have been shown to be critical, especially at small air gaps. For practical applications, 
the most important area which we have overlooked is that of robustness with respect to plant mode ing 
errors. This is an area which has also been a topic of current research (5]. In [3] it is noted that the 
electromagnet nonlinearity results in nonlinear cross-coupling terms in the control of five degrees o 
freedom of a precision linear bearing suspension. Thus it will be advantageous to implement nonli lear 
compensation laws for this and other multivariable suspension systems. 
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Abstract 

The control of systems with unknown dynamics and unpredictable disturbances has raised some challenging 
problems. This is particularly important when high system peformance is to be guaranteed at all times. Recently, 
the Time Delay Control has been suggested as an alternative control scheme. The proposed control system does 
not require an explicit plant model nor does it depend on the estimation of specific plant parameters. Rather, it 
combines adaptation with past observations to directly estimate the effect of the plant dynamics. 

This paper formulates a control law for a class of dynamic systems and then presents a sufficient conditiorf for 
control system stability. The derivation is based on the bounded input-bounded output stability approach using 
Loo function norms. The control scheme is implemented on a five degrees * I- freedom high speed and high precision 
magnetic bearing The control performance is evaluated using step responds, frequency responses and disturbance 
rejection properties. Tiie experimental data show an excellent control performance despite the system complexity, 

1 Introduction 

Some classical control methods deal with well known linear time- invariant systems. In many applications, however, 
some relevant part of the system maybe unknown, time varying, or nonlinear. Controlled systems are thus often 
limited to operating in only a small portion of their available range. For example, servo motors must operate in 
the linear part of their range for accurate control. Restrictions such as these have led to the developrrn nt of control 
techniques that deal with such complexities. 

Several types of modern control strategies have been developed to deal with nonlinear, time- variant systems. One 
of the first methods to accommodate nonlinear systems was Model Reference Control. This technique employs a 
model of the system and uses the difference between the model response and the plant response as the input signal to 
the plant [18]. The model is either a physical model or a simulated system on a computer. Although it has no variable 
parameters, it is very useful for either specifying desired performance or for the observation of unarcessabie states. 
A drawback in tins technique is that it requires knowledge of the full dynamic model and system limitat ions. When 
perfect cancellation of the system nonlinearities is not achieved due to imperfect modeling or inaccurate parameter 
values, the dynamic performance of the plant may be degraded to the point of closed loop instability [22]. 

Another advanced technique is Adaptive Control. An adaptive system measures a certain index of performance 
which is a function of the inputs, states and/or outputs of the system. From the comparison of the measured index of 
performance with a set of given ones, the adaptation mechanism modifies the parameters of the controller or the set of 
given ones [4,14,16,21]. i here are several classes of adaptive control. A very common variation uses a desired reference 
model as a basis for comparison and is termed Model Reference Adaptive Control (MRAC). In the direct MRAC, no 
attempt is made to identify the plant parameters. Controller parameters are directly updated. In Self-Tuning control, 
plant model parameters are identified/modified and the controller action is automatically updated according to a 



|,* ( ,| regulator design Another approach generates <1., control a. lion in part by an adaptive ^forward «*!" II 
which “behaves” a* the “inverse” of ll.e plant [22]. All adaptive controllers share the distinguishing feat.' re of s>s jm 
identification followed by variation of parameters to maintain desired performance. A drawback of adaptation 
£ b generally slow and computationally intense. Often the environment changes faster than ^*§ ^2 41 
performance degradation or even instability. Other references on adaptive control include [8,9,1 1,12,19 20,24 . 

Other control methods, such as Variable Structure Controllers, take totally different strategies to achieve stai.lt y 
in nonlinear systems. This type of controller utilizes state feedback in a control law which switches the structure 
of the closed Lop system between trajectories which may themselves be unstable or marginally stable but when 
combined by the control law in a switching technique, result in a system which is stable. A method of sw. - hing 

elites high frequency izjng the f am iliar features and patterns of a situation and which use past 

in [13] and l 17 l proposed in references [30,31,32,33], depends neither on estimation 

tint II nil II I iiiM i l 1 IMM f 

This paper uses the concepts developed in references [30,31,32,33] to explore tne poie 

~ w-j- * 

TlJ^^rCuh^r Z'p to robot manipulators »<l 

under large system parameter variations and disturbances [30,31,32,33,35]. Mammy 

also performed for linear SISO systems [34]. . t with non |inear input action and then presents 

This paper formulates a control law for a eta ol the bJnded input-bounded output 
a sufficient condition for control sys em s a l y. , • ; mD i eme nted on a live degrees-of-freedom high 

stability approach using L,„ ^ p.^^'iro/perrrorm.rtce is evaluated using step responses, frequency 

Sp e o d ^,"s and di.uZceTe]!"o„ properL, The «p.„n«„tal data show an excellent control performance de.p.te 

the system complexity. 


2 Time Delay Control 

In this paper we are concerned with a class of systems described by the following differential equations, 


i(0 = F(x,<) + G(x.u.O + D(<) (1) 

where x(t) € 7v n and u (/) € are the system state vector and control input vector respectively. F(x,<), G(x.u,f) 
and D(<) are vector functions with appropriate dimensions and represent respectively known dynamics, unknown 
dynamics and disturbances. The variable t represents time. In order to transform the syjtem into a familiar form, 
Equ.(l) can be written as 

x(t) = F(x,0 + H(x,u,t) + Bu($) (2) 

where the new term H(x,u,*) is defined as 

H(x, u ,/) = G(x,u,0 -I- D(<) - Du(<) (3) 


and B is a matrix to be selected by the designer. A reference model that generates the desired trajectory is chosen 
as a linear time invariant system, 

x m (0 = A m x m (t) + B m r(f) (4) 

where x m (0 G Tv" is a reference model state vector, r (t) € K r is a reference input. A m and B rn are constant 
matrices with appropriate dimensions. 

The class of systems considered in this paper includes systems that satisfy a matching condition. It was shown 
in references [30,32,33] that systems in a special canonical form satisfy the matching condition. These systems can 
be partitioned as follows 


x = 



X r 


F(x, t) 


x, 

F r (x,<) 


0 * 


0 


; B = 


, H r (x,u,<) _ 


. Br 


where the partial states are x q G Tv" _r , x r G Tv r , x, = [x r + i , x r +2 • ■ • > x n ] i 
following dimensions F r (x,/)> H r (x,u ,/) € 7£ r ,B € Tv riXr and B r € T£ rx 
the reference model of Equ. (4) are also partitioned in the same manner, 


G 7v n-r . The vector functions have the 
is of rank r The matrices involved in 



' o | I, ' 


0 

A-m — 


, B m = 



A mr 


. Bmr 


where I, G )x ( fi “ r > , A fnr G ft rx " 

matrix K of the form, 


B mr G 7£ rXr , B m G 7v nxr and r (t) G Tv r . Furthermore assume a feedback 


K = 


0 

K r 


where K G H nxn and K r G TZ rx ". 
according to 


The objective is to generate a control action u that forces the error to vanish 


e = (A m + K)e = A f e 


(5) 


The control action that combines past observations with adaptation for systems describ> d by Equ (2) is given by 


u(t) = B+[-x(f - L) + F(x,< - L) - F(x, f ) 
+A m x(t) + B m r (<) + B.i(f - L) - Ko(<)] 


:i!)2 


( 6 ) 



where the parameter L represents the time delay [30,32,33], the error vector e ,s defined as the difference between 
the plant and the reference model state vectors, 

e = x m - x < ' > 

The term B + is the pseudoinverse matrix defined as B+ = (B T B)-'B T . For the spec.al canonical form considered, 
B + is given by 


" o 


B + = (B T B)- l B r = < 


B r 


0 

B r 


[o ; Bj) 


= (B^Br) -1 [° ■ B r] = B r ’[0 : Ir] = [° : B r ‘] 


( 8 ) 


The control action now reduces to 

u(t) = B7‘[-Xr(< - IA + Fr(x>< - L ) - F r( x -0 + A mr*(0 

+ B r u(i — L) + 13 mr r(!) — K r e(<)] 

Note that this control law is a special case of a general algorithm which uses convolut.ons for estimating unknown 

SyS ' Thelb“! If this research is to be able to control such systems and guarantee performance despite the presence 

° f MGSS ““ .W control matrix D„ « 

.. t -V(* t) -x(t - L) + F(x t- L) + Bu(t - I) attempts to cancel the undesired known nonlinear dynamics 

Ffx o'lhe unknown nonlinear dynamics and the unexpected disturbances H(x,t). (3) the term A m x + m r inserts 
b ( x,t ), tne unknown nu feedback term -Ke adjusts the error dynamics. Thus 

useVin the time delay’control is reminiscent of numerical methods used to solve differen t equaUons. 


3 Stability Analysis 

3.1 Error dynamics 

-SSSSSSSSsr 


x(0 


F r (x ,0 + Hr(x,u,<) + [-Xr(t- i)+F r (x,/ L) T r(x <) 

+B r ii(< - L) + A mr x(<) + B mr r(<) - K r e(<)] 


H r (x,u,<) - Hr(x,u,( - L) + A mr x(f) + B mr r(t) K r e(t) J 

The previously defined error e of Equ.(7) is now governed by 

o(t) = (A m + K)e(t) + H(x,u,< — L) - H(x.u,i) (9) 

w |,ere the second M d third terms are foreleg functions due to the unknown „..em dynamics and „n,. red, r. able 
disturbances. 


Rewriting Equ. (9) as 


where 


e(<) = (A m + K)e(i) + p(0 


0 


P(0 = 


Pr(0 


( 10 ) 


and 


p r (0 = H r (x, u,t - L) - H r (x,u,0 

One may ask the question: what conditions does the vector p(/) have to satisfy for I lie system to be stable? A 
sufficient condition for stability will be derived in the next section. 

3.2 Sufficient conditions for stability 

This section presents a general solution to this multi-input multi-outp it control problem We will use the bounded 
input-bounded output approach based on [ JOO norms in order to derive sufficient conditions for stability. We now 
consider the governing differential of the error as given by Equ. (10) an I its correponding time response, 

e(0 = e (A "* + K)< e(0) + f c ( A m +K), t-T» p ^ r )dr 

Jo 

We will use ||(.) r || to indicate the norm of the time trancated function (.) and ||(.)l!» f° r induced i iatrix norm. 
Taking the norm of the error [28], 

||er|| < || (e (A ™ + K)*) r ||, ||e(0)||+ S up <€lo , T1 f' ||- (A - K,t '- r) IU|p(r)||dr 

< || (e (A " +K) ‘) T ll< l|e(0)|| + ||PT||sup (€[0 . r] .f o ‘||,‘ A - +K - ( ‘- T )||,dr 

The desired error dynamics given by (A m 4 - K) are always chosen to be asymptotically stable. This implies that 
there exist finite positive constants m, A such that 

|| c (A m +K)(<-r)||. < mc -A(t-r) Vr, t > r 
which implies 

su P<€[o,T] fo ||e (A -+ K)( ‘-e || f dr < su Pte(0 T ] f (1 - e" A ‘) = % 


|| ( e(A-+K). )T ||. 


s »P(€[0,T] ll(« (A “ + K) ')ll* 


< sup l€ j 0 7 ] me xt = m 


Therefore, the norm of the error is bounded 


||e T || < « + ^||pt|| t 1 *) 

where 

a = m|| u{0)\\, 0= y 

In order to be more specific on these stability conditions, we need to expand the forcing term p(0 We can rewrite 
IV (f ) as 



Pr(<) = 


H r (x(< - L),u(t — L),t — L) - H r (x(<),u(< - L),t) 

+H r (x(<),u(t - L),t) - H r (x(t),u(t)-0 

||[p(0]r||= ||[Pr(«<l|[H r (x(t-L),u(i-I),<-I) (12) 

-Hr(x(0,»(< - O.OItII 

+ ||[H r (x(0,u(<- L),t) — H r (x((),u(<),<) T \\ 


Assuming that the function H r (x,u,<) is continuous and difTerentiable, the Mean Value Theorem yields [6] 

||[H r (x(i - L),u(t — L),t - L) — H r (x(t),«(< - L).0]r|| 

< II [(tt* + WO.'* - l)X) ] t I|l 


(13) 


||[H r (x(0, «(< - L),l) - H r (x(t), «(0. 0 ]t|| 

< II [^ t W<). < 4( u ( ( )- u l<- L ))] r ll 


where E ( t — L : t) and 6 E (u(J L),u(t)). 

Equations (13) and (14) involve terms in x and u. 
used, 


To express them in terms of e, the following expressions are 


Q H r dH r . 

— + lhT x 


OHr 

Ot 




(15) 


u (t) - u(t - L) = B 7 ‘(-x r (< - L) + F r (x,t- O - F r (x,t) 


+A mP x(<) + B mr r(0 - K r e(<)] 

= B; x [-x m r(t - L) +e r (t - L) + F r (x,< - L) - F r (x,0 


(1C) 


+X m r(0 — (A mr T K r )®(01 

since e r = x mr - x r and x mr (0 = A mr x m (() + B mr r(i). Substitution of (13), (14), (15) and (16) in ( 
algebraic manipulation yields, 


2) and some 


llprll < ||2g* + tfeimWl' + ll^ll^ ll^ll 

+||^i t B r - l (Amr + K r )||,||e T || 

+ll^ t '®r Ml' l|(*"*r (0 — *mr(l — 0)11 

+ ll^ tB r _1 ||. IMI 

+ ||^B r - 1 ||. l|[Fr(x,t - L) - F r (x, i )] T II 


|| [F r (x, t — L) — F r (x, <))]r|| 



(18) 


where £ e (t - L,t). Substituting Equ. (18) in Equ. (17) yields 
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IIptII < c i + c 2 ll e T|| + C 3 ll®rll 


(19) 


where the constants Ci,C2 and C3 can be identified as 


c\ - 


1 *T B ' !li 


Ap Ap 

L + -^Xm|| + ll*mr (0 - X mr (t - I ) | 


, ,,dH r dH r . .. 

+ '- || -sr + '5r ,t ” 11 

||2tD; , (A„, + K r )|| i 


„ dH r . , ,, 5H r „_! „ 

C3 = 11 i^r 11 ' 1 + 11 '^r Br 1,1 + 


l ll^ll / 


We will assume later that the terms in the right hand side of the three equations listed above are bounded and hence 
the parameters , c 2 and c 3 will be bounded. The norm of the error in Equ.(ll) can now be evaluated and is found 
to be, 

||e T || < a + 0 HptII < « + 0c\ + 0 c i\\ e T II + /?ca||eT|| 


ll°r|| < 


or + 0c | + 0C 3 ||<!t|| 


if 0C ■< < 1 


(20) 


( l -0c,) 

The above equation relates the norms of the error and the error derivative. To obtain absolute bounds on the error 
another equation of this form is needed. The norm of the derivative of the error e, can be found from Equs. (10) 
and ( 19), 


Ml < ||A m + K||j Ml + IIptII 

< [ ||A m + K||, + c 2 ] ||er|| + e t + c 3 ||er|| 


This condition can be stated as, 


II ■ || ^ . [l|A m + K||i + e 2 ] it .I :<• ~ I 

IN II < jrzz ) + (T^T) l|eHI ,f C3 < 1 


( 21 ) 


( 22 ) 


and substituting Equ.(20) one obtains 


IMI < 


[ ||A m + K||, + c 2 ] (a + 0c i + /?c 3 || er|| ) 

+ ( 1 - 0c >)e l 


Ml < 


(l-c 3 )(l-0c 2 ) 

II A m + K||i(a + 0ci) + Cl + for* 
[1 — 0c 2 — c 3 — 0c 3 ||A m + K||i] 


if r ., 4. 0( r . 2 4 f3 || A m + K||< ) < 1 is satisfied. Using the expressions for the constants C! ,c 2 and r 3 , we have 


:t9<i 



( 23 ) 


(ll# : ll^ + ll^ D r 1 |l< + II^Br-'IMI^II.^) 

[1+ a ||A m + K||« ] 

+ ||^n L D" 1 (A mr + K r )||i x < 1 

Based on the foregoing analysis and relations between the vector functions H r , G r , and B r given by 

<9H r — dG r 0H r _ dG r ^ ^ 
dx ~ dx an du “ <9u 


the following sufficient condition follows, 


Theorem 1 , 

If the functions F(x,t),G(x y u ^),D(0 are continuous and differentiable, and 

If the Jacobian matrices and the vector functions SQz., 

If the eigenvalues of the matrix (A m + K) are in the left-half plane, and 

If the following stabthly condition is satisfied 

[1+ f IIAm + KIUtll^lliL 

+ ||^B r - 1 -I||i(l + II^IM)] 

+ IK^B- 1 - I)(A mr + K r )||i f < 1 




and 


(24) 


Then the time delay controller is stable, and 

the resultant bounds on the norms of the error, e, and its derivative e are 

a + fdci - 

" C " “ [1 - pC'i — c 3 - 0c 3 || A m + K||, ] 

„ ||Am + K||< (o + ) + ci+r 2 a 

l|G|1 " [l-/?C2-C3-/?c 3 ||A m +K|| i ] 


(25) 


The constants C\ , c? , c 3 a and 0 are, 


r i 


(£*' - 1)11, |t|2t+ *»ll +||X„„(0 - x,„,(I - 0)111 


c 2 = 


^3 


r t| aG r t dG r . 

+ t “-5r + ir x " 11 

Hl^D ,- 1 - I)(A„, + K r )||, 




du 


a = m||e(0)|| 


/? = T 


rn 

J 
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Figure 1: Region of Stability in the ||er|| — ||er|| space 


Condition (ii) implies that the rate of change of the functions F r , G r and D r with respect t, x and u are bounded. 
Condition (iii) implies that the desired error dynamics are chosen to be stable. Condition (iv) relates the time delay 
L, the rates of change of F r and G r with respect to x and u, and the desired error dynamics specified by (A m + K). 

The bounds on the norms of the error, e and its derivative e are shown graphically in Figure 1. Equ. (20) yields 
a straight line as the boundary between stable and unstable regions. Similarly, Equ. (22) yields another straight 
line. The common region bounded by the two straight lines is the region of stability as shown. The absolute bounds 
||e||mor and ||e|| max shown in Figure 1 correspond to those given by Equ. 125). 


3.3 Special cases and discussions 

Condition (iv) stated in the proposition 1 can be rewritten in a more convenient and usable form. This simplified 
and convenient condition, however, is more conservative. Starting with condition (24) and after some algebraic 
manipulations we obtain, 


'^ b ,- 1 - i ||. < 


{ 1 - U + T ||A m + K[|j] | 

/ it ! 11 dF r II r \ ^ 


( 26 ) 


(! + Ilex'll* L 

+ x [ ||A mr + K r ||i 

+||A m + K ||<( 1 + ll^-lli L )] 

The condition (26) implies a bound on the variation of ^ relative to the controller gain matrix B r The size 
of this bound is dependent on the delay time L, the norms of ||^ t ||» and and the desired characteristics 

of error dynamics given by (A m + K). The smaller the delay time is and the smaller the bounds on H^Hi and 

||^Lc-|| t are, then the larger the allowable size of bound is on the range of relative to B r . 

^n the case of first order SISO systems the vector functions reduce to scalar functions G r = g, B r = 6 and F r = /. 

Some interesting results are stated below for this class of systems. For first order SISO systems, we have 

|U(A™ + K)(r-r)||. _ | e « m +*)(t-r) | _ 




in — 1 

A = + k) 



||A mr +K r ||, = ||A m + K||, = \a m +k\ = — (a m + k) 
the stability condition (26) becomes 




dg_ 
i Oil 


< 




(27) 


The above condition implies 


l - 


- 2 »» , (I) 


3 + 


2 "s" i ) 


< i + 


('- 2 »I" l ) 

( 3+2 »i l|t ) 


(28) 


with the lower bound on 
becomes, 


dg_ 

du 


being always positive. 


For the case where f = 6 the stability condition (28) 
ou 


(29) 


ll|» L < 5 

This implies that as |||||| becomes large, the .IretcTdt 

intuitive sense. From (28), it can be observed that for sufficiently small l yu h 

the following limiting case, 

dg_ 

iif- - mi. 

dg 

2 


1 

< 3 


(30) 


3 6 


4 

< 3 


A 

The result indicates that stability is maintained for a variation of 66% of £ with respect to b. 

When the control distribution matrix # is a constant and known, the controller gain matrix B, may bech- 
1 ,h 2<L. - B This enables exact cancellation of the known dynamics F r and approximate cance 
unknown demies ind disturbances [33], The stability condition (26) then reduces to abound , on , de ay im 
terms of 2& and (Am + K). The known dynamics F r , the control distribution matrix and g 

matrix b7<Jo not enter the stability condition because of exact cancellation. 

4 Application: Control of a high speed and high precision magnetic 
bearing system 

The magnetic unde, 

environments such as integrated circui man g . * s un an induction motor. In order to 

The pump action is produced once the rotor, wi a es magnetically in the X Y and Z directions shown 

minimize impurities and particle generation, t e r0t ° r * SU *f summarized in Table 1 This system has five degrees 
in the figure. Some information relevant to this design are summarized in lable 

which may be described by a differential equation of the form, 

T r n 1 

(31) 


d 

' x v 


X r 

+ 

0 

+ 

0 

dt 

•• 1 

x r J 


_ F r (x, t) _ 

G r (x, u, t) 


D r (f) . 






Figure 2: A five-axis magnetic bearing 


Rotor Mass = 2.2 Kg 

Air gap for thrust, bearing = *100 ;rm 

Air gap for radial bearing = 250 /im 

Maximum current to bearings = 10 Amp. 


Maximum rotational speed = -15,000 RPM 

Table 1: Relevant system parameters 


where x<j E and x r E 7£ 5 represent displacements and velocities of the rotor with respect to the bearing stators 
respectively. The current inputs to the electromagnets are represented by the vector u € ft 5 . The control objective 
is to levitate the rotor and maintain stability. Also, the control system must reject disturbances under spinning and 
nonspinning conditions of the rotor. 

This plant is multi-input, multi-output with all five degrees of freedom unstable open loop. Disturbances and 
coupling include forces due to gravity, magnetic actions, unbalance and gyroscopic efFects. All of these effects will 
show up in the vector function G r . Note that since the magnetic force is proportional to the current squared and 
inverseley proportional to the gap distance squared the function G r (x,u) is a nonlinear function dependent on the 
state x and the control action u. 

The variation of the component of G P in the Z direction in terms of the gap x and control current u are shown 
in Figures 3, 4, and 5. Figure 3 shows that for a gap of 0.15 mm, the rotor acceleration corresponds to 30 m/sec 
and 1 10 m/sec 2 for control currents of 1 and 2 amps respectively. Also, for the same gap opening and with currents 
levels of 1 and 2 amps, §£ changes from about 0.22 x 10 6 ^ to about 0.88 x lO 6 ^, and changes from about 32 
m/amp-sec 2 to about 112 m/amp-sec 2 . It is clear that this particular device experiences drastic dynamic changes. 
Therefore, such dynamic information would be necessary if a conventional control system is used; otherwise the 
system performance may be acceptable only for some specific operating conditions. 

The control algorithm was implemented on a DSP chip as shown in Figure 6. In this experimental setup, we have 
the option of controlling the system using either a linear analog controller which resides in the compensation block, 


10(1 





or a time delay controller implemented digitally in the DSP board. The position signal for the Time Delay controller 
is obtained through the test points TP2 and/or TP3 and this signal is then sent into an analog to digital converter 
(A/D converter) which is linked to the DSP board. The A/D board has an adjustable built in low pass filter where 
the position signal can be filtered. The control voltage signal is sent out through the D/A converter which has a low 
pass filter with adjustable cutoff frequency. 

The program was written in C and assembly language. The sampling frequency used was around 5 KHz. The 
computation time for the control algorithm was about 70/isec. The cutoff frequency of the filter for the position 
signal was kept at 7.2 KHz since the signal from the sensor is already bandlimited to 1 KHz. The ^pos.tmn jgnaU 
are then obtained by the DSP board and the control actions representing currents are sent out through the D/A 
converter The D/A converter has a low pass filter with a cutoff frequency set to 1.5 KHz ; Although a samp mg rate 
of 5 KHz was adequate for this system, the sampling frequency could be increased up to lo KHz. This sampling ra e 
could further be increased by optimizing the program code and hence reducing the computation time for the Time 

DCl ?n section, we will use the control procedure described earlier to maintain a desire. I performance. ^ e r nodel 
reference for the thrust bearing position was chosen as a second order system w.th a natur d frequency of 200 rad/sec 
and a damping ratio of .707. The experimental data shown in Figure 7 indicate that tl e actual position response 
tracks the reference model response very closely. In this case, the position of the rotor moves from 200/mi to 0/im 
which corresponds to the suspended configuration. The error between the desired and ad ual position trajectories is 
shown in the same figure and has about less than 10% maximum error. The control current necessary to produce 
this response is also shown in that same figure with a maximum current of about 1.75 amps. This is an , . xce 
performance considering that the controller has no detailed information about the system F.gur ' . 

closed-loop frequency response between the reference position and actual position of the thrust bearing In 
it7s clerthat the magnitude and phase characteristics are close to those of the reference model -Acted . In orde 
to check the disturbance rejection properties of the control system, an additional current ,s injected through the 
drive amplifiers (Auxiliary input 2 in Figure 6) to create an intentional disturbance force. The frequency of this 
input is then varied from 0.1 Hz to 10 KHz (sine sweep). In order to check the disturbance rejection proper ie&, we 
measured the frequency response from the additional current to the position of the rotor. The disturbance rejec. ion 
properties of the thrust bearing are shown in Figure 9. This curve repres, -its , a compliance curve J^e controHer 

rejects disturbances up to the bandwidth which is again around 200 rad/se. I he static stiffness is a > • / 

and the minimum stiffness is about 300 KN/m at the frequency of 200 rad/sec. 

Figure 10 shows the closed loop frequency response for a radial bearing. Again, this response ,s between he 
reference position and actual rotor position. This is very similar to that of the reference modeh Figure 11 show 
disturbance rejection of the radial bearing when the rotor is at rest and while it n 

RPM 30400 RPM and 34800 RPM. When the rotor is not spinning, the static stiffness is about 200 MN/m and 
fhe m, n, ml stiffness is about 500 KN/m. It is clear that the disturbance rejection d/se T 

for these different operating conditions. Figure 12 shows the effect of using a lower a 
this case, the disturbance rejection curve moves up indicating a lower stiffness. These data demonstrate 

control scheme possesses excellent robustness properties. 


5 Conclusion 

an ] UnP t rtion The result of stability analysis performed based on the bounded input-bounded output 

The control performance, evaluated using ,„p response, an, I rhsturbanc, reject, on properues, 
shown to be excellent despite the complex nonlineanties in the system. 
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ABSTRACT 


Certain experiments contemplated for space plat forms must be isolated from 
the accelerations of the platform. In this paper an optimal active control is 
developed for microgravity vibration isolation, using constant state feedback gains 
(identical to those obtained from the Linear Quadratic Regulator [LQR] appn ach! 
along with constant feedforward (preview) gains. 

The quadratic cost function for this control algorithm effectively weights 
external accelerations of the platform disturbances by a factor proportional to 
(I/lC)*. Low frequency accelerations (loss than 50 Hz) are attenuated by go atei 
than two orders of magnitude. The control relies on the absolute position and 
velocity feedback of the ex oriment and the absolute position and velocity 
feed— forward of the platform, and generally derives the stability robtis tiess 
characteristics guaranteed by the LQR approach to optimality. 

The method as derived is extendable to the case in which only the relative 
positions and velocities and the absolute accelerations of the experiment and space 
platform are available. 



1. INTRODUCTION 

A space platform experiences local, low frequency accelerations (0.01 — JO Hzj 
due to equipment motions and vibrations, and to crew activity [lj. Certain 
experiments, such as the growth of isotropic crystals, require an environment in 
which the accelerations amount to only a few micro— g's [2]. Such an environment is 
not presently available on manned space platforms. 

Since the experiment, and space platform centers of gravity do not coincide, a 
means is needed to prevent the experiment from drifting into its own orbital motion 
and into t he space platform wall. Additionally, some experiments require umbilicals 
to provide power, experiment control, coolant flow, communications linkage, oi 
other services. Unfortunately, such measures also mean that, unwanted plat lotto 
accelerations will lie transmitted to the experiments. This necessitates experiment 
isolation. Passive isolators, however, cannot compensate for umbilical stillness, nor 
can they achieve' low enough corner frequencies oven il umbilical.'' are absent. 
Active isolation is therefore essential. 

The problem, then, is to design an active isolation system to minimize these 
undesired accelerat ion transmissions, while achieving adequate stability matgins and 
system robustness. Spatial and control energy limitations must also b< 
accommodated. 
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2. mathematical model 

The general problem has three translational and three lotational demur s of 
freedom. For simplicity, however, this analysis will consider only the 
one-dimensional problem. The general problem could be treated in an analogous 
manner. Let the experiment be modeled as a mass m. with position x(t). Assume 
that the space station has position d(t), and that umbilicals with stiffness k and 
damping c connect the experiment and space station. Suppose further that a 
magnetic actuator applies a control force proportional to the applied current i(t). 
with proportionality constant a. Such a model is shown in 1* iguie 1. 

The system equation of motion is 

mx r c( x — <1 ) + k( x — (i) 4- oi — 0 1 J 1 


Division by m and rearrangement yields 




m 


In state space notation this becomes 


x = A x + b u + f 


where 




■ “ >1 

, i) — * 

0 

> 

k e_ 


a 

m in 


m 


lift 



u 


The objective i 




f = 


in m 


s to minimize the acceleration x(t). 


3. OPTIMAL CONTROL PROBLEM 

The optimal control problem is that of determining the ontrol current 

u(t) — j which minimizes a suitable performance index 


for the system dcscriljed by Eqn. (3) subject t.0 the state variable conditions 


x(0) x () 

I i in x(l) - 0 
t— »oc 


( ')<i f 
\ oii i 


Another reasonable assumption is that f(t) is Ixmnded. and it will be h-ui.d 
mathematically advantageous (and only minimally restrictive) to assume that i( : ) t- 
also a dwindling I unci ion: 


1 i m 1(1 ) = 0 
1 — X 1 


A quadratic performance index 


.1 =i/ X [x 1 VV X + W..U-) dt 
“() 

has I non chosen, as one that lends itself well to the variational approach to optimal 
controls, since an analytical solution is desired. Tin* upper limit of the definite 
integral has been selected so as to yield a time-invariant controller. Here Wj is a 
square 2x2 constant weighting matrix while w. { is a weighting constant 



Although, Wj 


could lie a full 2x2 matrix, for this problem a diagonal lorm 


has been employed tor the sake of simplicity. 


W l = 


' la 

l) 


w 


lb 


The performance index consequently reduces to 


so that each state is weighted independently. 

If sinusoidal motion of i he experiment is considered, so that 


( 7 ) 


N I 


x(t) = H sin uX 


uid x(t) 


"x(t ), t he 


and cont rol as 


.1 




U/ 


) 

T \v. ? ir] dl 


It is apparent, that this performance index conveniently weight 
frequencies much more tin, at higher frcquennes. 


accelerat ion 


at a.. 
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1. SOLUTION 

binding the optimal control to minimize Eqn. (4) is a variational problem of 
Lagrange, for which the initial steps of the solution are well-known fe.g.. Elbert 
[4]). The variational approach is outlined Itelow, following which the complications 
added by the non homogeneous term f(t) will I)e addressed. Current optimal 
controls texts either assume that f(t) = 0 (e.g., [4], p. 2(12) or recpiire that, it have a 
restricted range space (e.g., [6], p. 238). The solution that follows provides an 
analytical optimal without imposing such restrictions. 

The argument of the cost function .) from Kepi. (4) is augmented !>v the 
Lagrange' multiplier A limes t-lu* system equation of motion Eqn. (3) where 


A = 


M 


( I IJ J 


The result .J can be expressed as 
■ x> 


II dt 


0 


where the Hamiltonian II is 


H = ^ (x ^ \V. x + w.jtt - ) 4- A* (x — Ax — bn — f) 


It is desired to obtain an optimal solution it = u which minimizes .J. 
I’he first variation of J(x, u. x) is 


M = 


■oo . 


[f Ax + «1 *| <1,. 

(I W Ok 


420 



wliicli is sot, ('(|ual to zero to minimize .J. 


However, intonating by parts, 


rOO 
• 0 


{^ 1 4 ) <it, = 

t>X 



ftX (ll 


0 


so that the ah >ve expression for W becomes 


« = 


•00 
• 0 



im 


Both Ax and hi are arbitrary variations, so M - 0 only il 


on _ ; t 

Ux~ ~ 


on _ 

W\- 


lllo 


The conditions given by Eqn. (5) still apply. 
Solving Etpis. (Ha) and (Mb) yields 

A = VV t x - A A 


(15a) 


(151) 
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* 

Temporarily eliminating n produces the result 



( 16 ) 


where 


A = 


A 


W, 


i_ 

w. 


I) I * 1 


3 

-A 


If K<|». (1(1) is now solved for A in terms of x and f. 1- tpi. ( lob) will then furnish an 

f 

expression for the optimal control n . 

As noted before. optimal control texts generally treat the homogeneous 
problem (where f(t) = 0), but they do not provide an analytical solution to the 
nonhomogeneous system described by (a) and (16). Salukvadze has treat ’d the 
non homogenr ‘Ot is problem [-1,5], but his dilficult treatment seems largely e> hate 
remained either uneomprohended or under— appreciated. 1 his me hod is especial:;. 
\v<‘1 1 — suited to low— fre<|ueucy disturbance rejection, and has been applied below ’o 
the present problem. 

Tlie homogeneous solution to Eqn. (15), where i = U, 



ATI 
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The four eigenvalues of A may be found to be, in ascending order of real parts. 


2 1/9 - l/“ 

-/?! + (^i - 4/y 1 '- 


(ISa) 


-i\ - - 4/y J 


U8b) 


fl i ~ fl l 


(18c) 


/'■.( = ~t l 9 


(18d) 


vliere and J., are defined as follows: 


A = 2k _ £ “ _ ' V w lb 
1 m m- m w 3 


(19a) 


M> = 'V 


a w. , i 

la, k 

9 + •> 

m“vv.j m J 


The eigenvectors of A corresponding to the respective eigenvalues may be chosen 


to be 


* - i + Jil + ■' l7 - > + ^ 

\ Vk 


(20a) 
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where 7^, 7.), 7^, and 7^ are defined below: 
k 


(20b) 


n 

7 3 ~ 2 

m w. 


( 20 ( 1 ) 


7 t = w la 


l 'sing E(|iis. (18) through (20) with (17) the solution to the homogeneous sysi.Mii is 


x c.e 



1>1 + (: 2° ** H*> +c f ” Hi 


‘l - “i " 1 * 1 1 

/y iLyi -/y -lb 1 

Cj e j>^ + e 2 e “ ^ + c.,e li L> 


| Uk l 1 , k 

1 K J 


1 and where e I are arbitrary constants. 


Application of the variation of parameters method with terminal conditions 
(Eqns. 5b. c) leads to the general solution of the noil-homogeneous system, with two 
constants of integration yet undetermined. 

If the two constants of integration are eliminated by solving for A in terms of 
x and f, the general solutions for A. and A., become: 


\5 l T '6 2 T m 


-lb' 


+ V “ 

(22a) 



+ ^8 e " 

(22b) 
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in which the £ } 's are functions of the eigenvalues and eigenvectors of A. and of the 
disturbance f(t). 

The Solution Form 

Using the fact that 


* 1 T 

U (t) = i r A l ll 


(cf. Eqn. (15b)] 


3 


(23) 


the optimal control is found to be 


— // . t, /i 1 t ~lb l . r , , 

U (t) = X | + 7/. 2 x. 2 + Ilf I <> U' ) (lt + 'If “ / 0 “ f 2 U)dt 




(24a) 

where 

'll = “7i ( m -, V 1 2 ) 


(24b) 

h = A? < iTT + "l + "2> 


(24c) 



(2 Id) 



(24c) 

(It should be noted that the feedback gains n { and ;/ 2 are those which would 
result from applying standard LQR theory to the homogeneous system equation 
x = Ax + bu). In Eqns. (24) p v /t, are the eigenvalues of A with negative real 


parts, [see Eqns. (18a,b)] and 
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y t > = £ a + ra a 


(24f) 


By repeated application of the method of integration by parts, the control may be 
re— expressed in terms of an infinite sum: 


» (t) = '/[ x i + , /-) x 2 + 


oo 

V 


r = 0 




+ \ 


oc 

V 

r = 0 


-D r 




/'•> 


r+1 


(25) 


Rewriting f., in terms of d and d, the cont rol function becomes 


u m = // , x ( t. ) + //._,x(t) + [jYf + ~i] *•(»■) 


r— I 


+ 


i = i 


M| M*) 

/.-I c ( 3 + 3 

/‘l 1 /'•> 


m ^ T + ' + ^ ^ m ^ i + l + i + P 

Ml M\> 


d (i) (l 


, [( - ' il ( — fi— _j 1'_. j| d 1 11 1 1) + liif’hcr order I e i ml- 


/* i /*2 


(2(i) 


This may be written in a more appealing lorm as 


„•(..) = C M x(t ) + C v x(t) + c (10 d(t) + c ()1 (1(t) + higher order term 


dl 


127) 


in which the constant coefficients c , c v , c (J() , and c (Jl may lie defined from lapis. 
(24) and (26). Clearly, if the infinite sums converge rapidly enough, the optimal 
control can be approximated by 


426 



For very low frequency disturbances the higher order terms in Kqn. (26) are 
negligibly small, and the control (Eqn. (2B)) closely approximates the optimal. II. 
in fact, the second- and higher-order derivatives of d(t) are identically zero, the 
approximation is exact. It can be shown that for the critically damped closed loop 


system the eigenvalues are real and equal, and that the convergence is 
than for the overdamped system. Further, as the closed— loop system 
become more negative the convergence speed goes up as well. 


more rapid 
eigenvalues 



5. CONTROL EVALUATION 
Physical Realizability of the Control 


The control, Eqn. (25), is physically realizable, if the states and sufficient 
derivatives of d(t) are accessible (or estimable by an observer), and if the higher 
order terms are negligible. It is not necessary that the eigenvalues be real, although 
the proof of this requires a more general linear— algebra or state— transition— mattix 
approach. 

If values are assigned to the system parameters, associated controller gains 
can be evaluated. Suppose that m = 100 Ibm, k = 0.3 Ibf/ft, c = 0 Ibf— sec/1 1. and 
o = 10 Ibf/ Amp. With w. ? arbitrarily set at 1 and w 1|} vari >d. associated integer 

values of w, can be found Ixilow which the eigenvalues /t. .md //., will always be 

I 

real. Such values are tabulated in Table I. Stated otherwise, the tabulated values 
of the weights W| and \v^ are those integer values (lor the oike ol simpli(ity) lot 
which the closed loop system is closest to being critically damped without being 
under damped. Corresponding controller feedback and leed— lorward gains (lor the 
first five derivatives) are also included. 

The states x(t) and x(t.) and the derivatives d l,l V). d f 1 *11 ' and d ^ ~ 1 1 1 ) are 
clearly available for an earth-based system. However, in space, tin* only absolute 

measurements which can be directly available are x(t) and d(t). itom width x(t). 
t’l(t) and x(t), d(t) are obtainable only by successive integration(s). Rearrangement 
of (28) into 


«*(.) = (< p + < <! 0 )x(t) + (c v + i',||)i(0 -<-,|„|xC' 

(2ft) 

or 

u*(t) = (c + c (I0 )d(t) + (c v + c <u )(1(t) + c p [x(t) — d( t )] + c v [x(t) -a(.)| 

( 30 ) 
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obviates the need for one accelerometer, but one accelerometer plus two integrations 
remain necessary for either the platform or the experiment. Since IxitHlftl] (or 
one of its integrals) has not been weighted in the performance index .1. experiment 
drift will be a problem that must be corrected either by another control loop or by a 
change of system states. The latter could be accomplished by incorporating an 
accelerometer attached to the experiment into the stale equation. Alternatively. 

one could append a tegtator to the plant, include the cum as a third state. 

and optimize the control di/dt. Hut for the sake of simplicity (i ".. fewer slates) the 
foilin'!’ lias been assumed (without development.) in this papei . 

The higher order terms of the control [Eqns. (’in) and -f>)] can be neglected, 
for low frequencies, if the eigenvalues and //^ are of sufficient. inodnln>. Hum 
eigenvalues, in turn, are under the control of the de: igner. determined by bis Hioice 
of weights w, a , w )b , and w.,. It is apparent from «qn. (25) that u (t ) «*se..tially 
reduces to two alternating power series. For a sinusoidal disturbance of Irequencv - 
the series form of the control converges for |u,'/yc| < 1 (1 = !--)• 11 (iU1 l,e >:ioul1 

that each alternating power series convi'rges like ^ ( -1 ) - j) • ^ lt ' 1 ' <>u 


frequency disturbances (i.e„ small relative to system closed loop eigenvalues) a 
control formed by series truncation very closely approximates the optimal. 

For example, suppose that the normalized frequencies | -///-, ! fur •> Hnusoid,.! 
disturbance are less than 1/5, and that only the feedforward control terms c i() du ) 
and e dl cl(t) are included with the feedback terms. Even so. the feedforward portion 
of the truncated control, at any time t, will be a current that is still wit bin •! . [m\. 
(1/5)“ of the feedforward |H)rt.ion ol the actual optimal. 11 tin noimaliztd 
frequencies are below 1/10. this approximation error will be less than 1%. 1 able 1 

shows that the gains c (]i of higher order derivatives d (,) (t) [see Lqn. (-’<>) h>i 
algebraic representations] are, in fact, quite small. 


O’ H ^ P ^ ^ pc ^ 
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In some circumstances there may be design constraints which prevent the 
designer from selecting weights that will lead to sufficiently rapid convergence. 
However, since convergence occurs rapidlv even for eigenvalues of relatively small 
modulus ( | ufa. | < 1/3), in a great many cases the designer will have much 
latitude in his choice of weights. For "low" frequency disturbances, in t hese cases, a 
control which includes only one or two feedforward terms will be "close" to the 
optimal. These frequencies will be well— attenuated. 

Higher frequency disturbances will also be well-attenuated, provided the 
input— to— output transfer function(s) are at least strictly proper in the I. a place 
Transform variable s. This will not be the case for the present problem if more than 
three feedforward gains (c^, c ( j|, C ( |.>) are included in the control. Practically, this 
means that only proportional and first— derivative leedtorward [I.\qn. (2-)j with r = 
0,1 or Hqn. (26) with n = 2] should be added to the feedback control terms. As will 
be seen shortly, however, adding even the proportional leediorward tertii(s) can 
dramatically improve the disturbance rejection over that afforded by LQR feedback 
alone. 


'Transfer Function and Block Diagram 

Neglecting the higher order terms, the transfer (unction Ix'tween input and 
our put accelerations or displacements is 


s-X(s) _ X(s) 

“ "W 


( — - c | , ) s + ( — - r |ii 
' o dl o dt) 


(— )s“ + (- 
'o' o 


V s 


)S + ( 


o 


and a block diagram of the controlled system can be drawn as in Figure 2. 



Control Stahilit.v. Stability Robustn ess, and General Robustness 

Since the control feedback gains are the same as those obtained by solution of 
the standard Linear Quadratic Regulator (LQR) problem, the closed loop system is 
stable and enjoys the stability robustness characteristics guaranteed by the (LQR) 
approach to optimality, viz., a minimum of 60° phase margin, infinite positive gain 
margin, and 6 dB negative gain margin [6]. Additionally, numerical checks indicate 
that it enjoys substantial insensitivity, or general robustness to uncertainties in k, < . 
and m, as indicated by Table 2 and Figures 3 through 10. By comparing the Bode 
plots of Figures 3, 5, 7, and 9 (corresponding to controls using both LQR F/B and 
proportional F/F) with those of Figures 4, 6, 8, and 10, respectively (corresponding 
to controls using LQR F/B only), one can see that adding feed-forward 
substantially improves disturbance rejection at low frequencies. For example a 
comparison of Figures 3 with Figure 1 indicates that the optimal control method 
described above can lead to acceleration reductions ol greater than tour orders ol 
magnitude for all frequencies. This reduction is more than two orders of magnitude 
below that afforded by LQR feedback alone at the lower frequencies, i.e.. those most 

heavily weighted in the performance index. 

The order of the reduction is eventually limited by control cost, ol course, 
probably in terms either of actuator-related limitations (such as heat-removal oi 
f orce _g encra tion requirements) or of power limitations (especially in a space-station 
environment). The control also leads to displacement reductions of the same 
magnitude, limited in this case by actuator-stroke or spatial limitations. Providing 
a unit transinissibility for very low frequencies and weighting (x-d) and/or / (x— cl) 
in the performance index J would be steps toward addressing these latter 

limitations. 


Computational Aspects 

A significant amount of algebra was required to solve the two-state problem 
of this paper, and the labor involved increases dramatically with each additional 
state. However, such symbolic manipulators as MACSYMA may be used to ease 
the workload if a symbolic solution is desired. Further, well-known numerical 
methods exist (i.e., Potter's method [7] or Laub's method [8]) for solving the 
solution to the homogeneous system. These can readily provide the feedback gains 
in numerical form, even for problems with many states. It, might be anticipated, 
then, that a numerical method also exists for finding the desire 1 feed— forward gains. 
Such is the case, as will l>e shown in a later paper. 



6. CONCLUSIONS 

This paper has applied an existing method for obtaining an optimal control 
to the microgravity platform isolation problem, for which the disturbances to be 
rejected are low-frequency accelerations. The system was assumed to be 

representable in the form x = Ax + bu + £, with quadratic cost function 

j = * °°(x T W.x + w..u 2 )dt and diagonal weighting matrix W The. resultant 

J 0 1 

control law was found to be simple, stable, robust, and physically realizable. 
Further it was shown to have excellent acceleration— and displacement— attenuation 
characteristics, and to Ik' frequency -weigh ted toward the low end of the acceleration 
spectrum. 

The method is extendable to the case for which only relative positions and 
velocities, and absolute accelerations, are available; and can be applied >o .<s to 
weight relative displacements in the performance index. 

The approach as presented is algebraically intensive, but symbolic 
manipulators can be used to ease the algebraic >abors. further, siiue the method 
produces feedback gains identical to those obtained by the LQR approach to 
optimality, numerical computation of those gains is easily accomplished, even lot 
large systems. The feed-forward gains can be found numerically with comparable 

ease. 
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Tab lo 1. Optimal F/F .and F/B Gains for So looted State 
Variable and Control Weightings. 
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Table 2. Closed loop transfer functions for system with design 

parameter values of k = 0.3, c = 0.000622, and m = 100; but 
with actual parameter values as shown. Gl, G3, 65. and Gi 
include both LQR F/B and proportional F/F; G2, G4. G6. and 68 
include LQR F/B alone. Weighting parameters used were 
W j = 258, w lb = 10, w 3 = 1 (see Table l). 


System Parameters 


,. / lbf \ 
k { TT> 

, lbf-sec^ 
c <-ft— 1 

m (Ibm) 

0.3 

0.000622 
(C=0. 17.) 

100 

0.3 

0.000622 

100 

0.45 

0.000622 

100 

0.45 

0 . 000622 

100 

0.3 

0.00622 

100 

0.3 

0.00622 

100 

0. 15 

0.00622 

90 

0.45 

0.00622 

90 


Closed Loop Transfer Function 



^X(s) 


s 2 D(s) \ 

i 

— — — i 

Gl(s) 

0.0000622s + 0.0001 S 

0.31056s 2 +4 ,4675 s + 16 .062 1 

C2(s) 

0.0000622s + 0.0300 
0.3 1056s' 2 +4 . 4675 s+ 16 . 062 i 

G3(s) 

0.0000622s +0.0151 
0.31056s^+4.4675s+16.07 < 1 

G4(s) 

0.0000622s + 0.0450 
0 . 31056s 2 +4 . 4675s+ 16 . 077 1 

G5(s) 

0.000622s + 0.0001 
0 . 3 1056 s 2 +4 . 4680s +1 6 . 0624 

G6(s) 

0.000622s + 0.0300 
0 . 3 1056s 2 +4 . 4680s+16 . 062 1 ' 

G7(s) 

0.000622s + 0.0151 
0 . 27950s”+4 . 4680s + 1 6 . 0 < > 1 

(i 8 ( s ) 

0.000622s + 0.0150 
0 . 27950s"*- 1 . 4080 s *»- 1 0 . 0 m 1 
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Figure 4 
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Periodic Accelerations 

Rel 
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Crew Motion 


2E-4 

17 
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Gen eral rrooiem : 

In the upcoming space 
station, planned by NASA 

tor completion in the I^^Os, 
minimize the low frequency 


accelerations traps mi tt cd 


from the space station to an 
experimental platform contained 
on (inside) the space station. 

“minimize": reduce to ~ 10 

• L I 

it possible 

“low frequency": o.oot to 20 Hz 
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O bj e c Hr i ve : 

Find the << bcst**i(+) l 

to minimize x(t). 
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Sta te Eq nations : 



o I 

-k -c 
m m 



O R 


x = Ax + b u 

where x = 


Js-d 

m 


x(-t) 

x(t) 


A*f-I -il 

# m m J 


‘•(5I-U4 
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Problem Statement: 


ROMAC 


Determine the control u(t) 

which minimizes the performance 
indey ^ 

1 J * *2 j (x T W,S +U T Wu)olt 


for 


SU 


the system 

X * A X + 3 u 


hject to the conditions 
\ X (o) * X o 
I I i m X (t) — O 


I i m f (t) = O 
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ROMAC 


S elution AA e thod (Differentia I 

Equations Approach) : 

1. Augment "the performance index S 
with the state equations using 
Lagrange multipliers. 

2. Take the I variation 8 J of the 
augmented performance index T and 
Set i *ir ec^ual to zero: 


s J = / «- f£- Sx]d+-o 

where H = (x T W, x + u T W 3 y) + A T (* - A X- B u - f j 

3. Integrate the third term of the integrand 
by partSj combine ttrms J and set coefficients 
of the arbitrary variations £* and Su 
ec^uat to zero. 

Result; ^ = W, X - A A 

a - W,-' B T * 

4. Substitute for u in the state equations, 
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r 


Th« solution (!■*; the optimal control U = “*) 
is now u* = W,"'B T X whore X is ■found by 
Solving the i yiteim 

(i\ u bw-'b t i r- j rfj 

W“h -a t \U] N 

subject to 

X (©) » X* 

| i m X (t) = o 

i i m fW" 2 


* 

c Find the solution of the homogeneous system 

■ 'Ym 


Of 5 


W, -A T 


i. Use -the variation of parameters method 
■to find the general solution of the 
nonhomogeneous system, i.e., of 

m_[ A Bw 3 -B T im + ; 

Urk -A T JIM l«I 

7 Apply the terminal conditions on 

conclude that n of the In arbitrary 
constants in the general solution are 
(•dual to zero. 
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9 . Sol yt -for X M in ’term * ©*f X (t) in a manner 

Such that the remaining n arbitrary 
constants are eliminated. 

Use the equation 

U** w 3 -'b t * 

to -final u* in terms of* X. 

Result: 

u*(+) = (W 3 '' 8 T X lt x;')x 

+ (w;' b t C) e ' At x^'fw d+ 


where f ^ ^ 1 is the Jordan 


L°-aJ 

Canonical Form of* tine Hamiltottwe 

matrix [A BW;'B t 1 

J ) 

w. -a t 

• 1 4 

t 

where A contains only the negative 
eigenvalues of the Hamiltonian 
matrix^ corresponding to the 
eigenvalues of the closed-loop 
System (assuming | control! able 
where X- X « is the eigenvector 

matr tx which lceicl* To The 

ctbov & J* C* F» « 

' / 
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where. X “ 


c-o 

x‘*° 1 

II 

11 

<-•) 

V <*° 

11 

An J 


an< J where integration «"f "the indefinite 

integ ra 1 requires const ants 
of integrate on that a re all 
i denti cally zero, 

lO. Integrating repeatedly by parts and using 
the fact* that 

ca t c t - *- a-bw/'b t p 

where Pa-X^X ^ 1 [P is the solution to the 

well-known Algebraic 
Riccati Equation] j 

develop equivalent "Forms "for **M * 

U*W- -W 3 -'B t Px + w," B t 


.-W,-B T P 8 - W J -B T X^- , Z o (-A-)' 

<r) 


= - w:‘ B r P X - Wj 1 B T z(-A- T ) Pf 

3 rso 

= - w," b t p * - w;‘ b t x,; t (rl 

Note*. A state transition matrix approach yields 

U*ft) » -W 3 " B t P x -W 3 " , 8 T X‘;°"e' A+ / « AT X“f«dr 
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Solution : 

u*tt) = -W;'B T P x + Wj" B T X‘' I ""e- At I e A+ X';°f <•*) dt 

= - w," b t p x - w;‘ b t xJJV**/ e AT x{;‘V w dr 
= -w:'b t p x - wr b t C i (-A-rX;Y rt 

= -w;'b t p x - w : 1 b t i (-A- T ) r+, p? (r) 

J J r»o 

» - w ; 1 b t p x - w:' B T x‘ T £ (-a-)"' xj f ' fr) 

3 — j 11 r»o 41 


= - w:' b t 


x; f <r) 


Dropping higher order -terms (for r > o) : 


u*(t) = -w:'b t p x + Wj 1 8 T Xjj ” 1 A” i 


= - w;'b t p x + w; , B T xi;°- , A- x~’Pf 

= - w; 1 8 T p x w;' b t x,; t a _t xj p f 

= -w:'b t Px * w;'B T x- T A' T x 1 T ,f 


3*(t) > - W," 1 B T P X + W 3 ''8 T A- T Pf 

where A' T = - P ( PA + W,r^ (a- B W^ p)' 

n(x u Ax;;)- T = x:;""a-'x 1 < ;" 


ese are 


several forms for the control law. 
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where A*A-BW."B T P 
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roblem : 



i(t) = C p x(t) +• C v x (t) 

+ c d j(t) + c d a(t) 

o u i 

where Cp , C V| Cj Q f and C d| 
are constant gains. 
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Table I Optimal F/F and F/B G, 
Variable and Control 


System Parameters: 

m = 100 lbm k = 0.3 lbf/ft 

c = 0.000622 lbf-sec/ft ( ( = 0.1%) 
a = 10 lbf/amp 
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ins for Selected State 
Weightings. 


K/K Cains 
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Table 2. Closed loop transfer functions for system with design 

parameter values of k = 0.3, c = 0.000622, and m = 100; but 
with actual parameter values as shown . G 1 » G3 , 05 , and l»i 
include both LQR F/B and proportional F/F; G2, G4, 06, and G8 
include LQR F/B alone. Weignting parameters used were 
10, Wo = 1 (see Table 1). 


la 


= 258, 


'lb 


Closed Loop Transfer Function 


System Parameters 


k ( lbf 1 

k ( Tr } 

, lbf-sec-, 

c <-ir— > 

m (lbm) 

0.3 

0.000622 

(C=0-17.) 

100 

0.3 

0.000622 

100 

0.45 

0.000622 

100 

0.45 

0 . 000622 

100 

0 . 3 

0.00622 

100 

0.3 

0.00622 

100 

0.45 

). 00622 

90 

0.45 

0 . 00622 

90 


s 2 X(s) 

s 2 D(s) 


G1(S) 

0.0000622s - 0.0001 

0.31056s 2 +4.4675s+l6.0624 

G2(s) 

0.0000622s + 0.0300 

” 0 . 3l056s 2 +4 . 4675s+ 16 . 0624 j 

G3(s) 

0.0000622s + 0.0151 j 

0.31056s 2 +4.4675s+l6.0774 1 

G4(s) 

0.0000622s + 0.0450 

0.31056s 2 +4 .4675s+l6 .0774 

G5(s) 

0.000622s * 0.0001 
0 . 31056s 2 +4 . 4680s+ 16 . 0624 

G6(s) 

0.000622s + 0.0300 
~ 0.31056s 2 +4.4680s+l6.062 1 

i 


0.000622s * 0.0151 
0 . 27950s" ->-4 . 4680S+16. 077 1 


0.000622s + 11.0450 

0.27950s" -4. 4680S+ 16. 077 1 
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C onclusions i 

1. An optimal con+rol ha* b««n dttcrmmeiJ 
for file nonkowojentou s LQ R problem. 

2. An approximat ion to this optimal control 
Via* been -found which uses constant 
feedback and feedforward gains. 

3. The optimal control has the -following 
a d v a nt ages: 

a. The gains can be easily determined. 

b. The control is very robust (CO phase 

margm J infinite positive gain margin ) 

* * « \ 

£ dB negative gam margin;. 

The control is applicable to a wide 
range of problems. 

d. The control offers substantial 

improvement in disturbance rejection 
over -that afforded by LQ R feedback 

at l© n e . 

e The control can be easily implemented. 
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Abstract 

Numerical calculations and experimental measurements of forces from an actuator 
of the type used in active magnetic journal bearings are presented. The calculations are 
based on solution of the scalar magnetic potential field in and near the gap regions. The 
predicted forces from a single magnet with steady current are compared with experimental 
measurements in the same geometry. The measured forces are smaller than calculated ones 
in the principal direction but are larger than calculated in the normal direction. This 
combination of results indicates that material and spatial effects other than saturation play 
roles in determining the force available from an actuator. 
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1. Introduction 

Recently, there has been much interest in the use of active magnetic bearings to replace 
or augment traditional bearings in turbomachinery. Magnetic bearings offer a number of 
potential advantages, including low power loss, suitability for harsh environments, and the 
ability to change the bearing characteristics to minimize vibration. 

This paper presents data on calculation and measurement of forces from a magnetic 
actuator similar to those that are used in magnetic journal bearings. The data are for open- 
loop, steady conditions only and are part of the development of models for active bearings. 

Much research has recently been devoted to magnetic bearings. No extensive literature 
survey is attempted in this paper, but from a sampling of the published papers on the topic 
[1-9], and the proceedings of the two international symposia on magnetic bearings [ 10, 1 1 ], 
it is clear that the concentration is on control aspects. Although the development of robust 
control strategies is important in optimizing magnetic bearing characteristics, particularly 
stability characteristics, it is also necessary to understand the forces arising from actuators 
in order to optimize fully a magnetic bearing system. This is expected to become 
increasingly important when the requirements for peak force and for force to weight ratio 
are stringent, which may be expected in aerospace applications. Better understanding of 
the characteristics of the journal bearing actuator is necessary in order to take advantage of 
all available parameters, including gap size, shape and material selection in a system 
optimization. This paper presents preliminary results of work toward this goal. 

2. Theory 

A computer program has been written that calculates the force exerted on the journal 
by a magnet having a steady current in its coils. The force is found by calculating the 
energy stored in the air gaps between the magnet and the journal, then performing a 
numerical perturbation to obtain a central difference of the energy change per unit position 
change. This gives the force in the direction of the perturbation. The force for a magnet at 
an arbitrary location can be calculated. 

The algorithm includes the following assumptions: 

i. The permeability of the metal is infinite compared to that of the gaps, which is 
assumed equal to that of free space. This implies that all the energy is stored in 
the gaps. 

ii. All magnetic flux closes the path through both metal parts. Expansion or 
fringing of the field near the gaps is allowed, however. 

iii. The coil current, therefore the MMF, is constant over a perturbation. 


In an isotropic domain not containing currents, where time variations are only of 
low frequency, the magnetic field can be represented as the gradient of a scalar field <Kx,y). 
The energy contained in the domain is 


O = B*BdV 

Jv 

where the flux density B is given by 

B = - V(J) 

and the potential (|> satisfies the governing equation 


( 1 ) 


( 2 ) 


V<}> = 0 


with the boundary conditions 


at> 


= 0 on free boundaries 


(3) 


an (4) 

and, because of assumption (i) above 

4> = 4>1 on pole face 1 

= 4>2 on pole face 2 (5) 

4> = 0 on journal surface, 

as shown in Figure 1 . 

Initially the boundary values of <J> on the pole faces, <t» j and d>2 , are not known, 
but must be determined in relation to the datum of<J>=0 on the journal surface. The problem 
is made tractable by the fact that the governing equation is a linear one, so that the values of 
<J> internally are determined within a multiplicative constant even for an arbitrary choice of 
boundary condition values. The fact that the flux must be the same through the two gaps 
allows the ratjo 

$2 


K = 


$1 


( 6 ) 


to be determined. Then the fact that the difference between (he two potentials is the 
magnetomotive force, 

<J>l-<t>2 = C (7) 

allows determination of the actual surface potentials. 


480 



2.1. Computer Program 

The algorithm above is embodied in a FORTRAN computer program, GAPFOR1 , 
which uses the finite element method for calculating the magnetic potential in two 
dimensions. For a given journal position the program calculates the gap height as a 
function of angular location and generates a finite element mesh for each gap. Flux 
fringing is allowed by extending the finite element domain beyond the edges of each pole 
face. Then the journal position is perturbed four times, first with positive dx and negative 
dx, then with positive dy and negative dy. At each step the mesh is regenerated and the 
energies are recalculated. The central difference analog to the derivative of energy is then 
computed, which is equal to the force in the perturbation direction. 

To achieve rapid computational speed and efficiency, a dedicated finite element 
program was written for this application. It includes a grid generation routine as well as a 
banded gauss elimination solver for the assembled equations. 

Forces from one magnet of a bearing 

The computer program has been used to predict the forces from one magnet acting on 
the journal at various positions of the journal within the clearance space. The geometry 
corresponds to that of the experimental apparatus described in Section 3. Half of the entire 
clearance space is mapped, since all positions of the journal with respect to a single magnet 
can be represented in terms of positions in this half plane. Figures 2 and 3 show a map of 
force versus x,y position. The magnet is the upper vertical magnet, and a steady current of 
1 ampere through the coils is used. The dimensions and other parameters are the same as 
those of the experimental apparatus described below. The figure indicates that the force in 
the y-direction (the principal force) varies between 6 and 180 N as the journal is moved 
along the y axis between -0.7 < y/c < 0.7. When the journal is also given an x-direction 
eccentricity, the y-force increases somewhat. 

Except at x=0, there is also an x, or normal, component to the force, as shown in 
Figure 3. This normal force increases rapidy as x is made larger. At x'c = 0.7 it is 12% of 
the principal force. 

3. Experiment 

3. 1 Design of Apparatus 

The apparatus for force measurement is shown in an exploded view in Figure 4. Each 
magnet is independent and is wound with 400 turns capable of carrying current of 2.0 A in 
the steady state. In the steady force measurement mode, the rotor is held stationary by 
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Upper vertical magnet 



Figure 2. Attractive force from one magnet in principal direction. 
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Upper vertical magnet 



Figure 3. Attractive force from one magnet in normal direction. 


Exploded View of Experimental Apparatus for the 
Research: Vibration Control in Turbomachinery Using 
Active Magnetic Journal Bearings 
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Figure 4. Assembly drawing of force measurement apparatus. 
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pressure from six micrometer heads (three on each end) that are in turn held by cantilever 
arms instrumented with strain guage bridges. Thus all mechanical force on the rotor passes 
through the micrometer pushers and the strain guage arm transducers To minimize any 
tangential force on the pusher ends Teflon ball sockets with steel spheres were used 
between the pusher micrometers and the support disks. A full assessment of the 
uncertainty due to friction has not been completed. 

3.2 Results of Measurements 

Forces were measured at several locations and for several values of steady current. 

The figures referred to below display dimensional data as measured, with forces in 
Newtons plotted against y/c, the eccentricity ratio in the vertical direction. All of the forces 
presented are from the lower vertical magnet, so the vertical forces are in the negative y- 
direction. The eccentricites in the x-direction are all positive. Three traverses of the y- 
direction were made, at x/c positions of approximately 0.0, 0.24, and 0.45. Assessments 
of the errors in measurement are not complete; however, it is expected that the error in 
position measurement is no greater than plus or minus 0.05 in y/c and x/c, and that the 
error in force measurement is no greater that plus or minus 5 N. Errors in current level 
control are within 0. 1 A. A larger series of measurements that were made before the 
addition of the balFsocket contacts was eventually discarded because the measurement error 
due to friction appeared to be significant. 

The data support some of the anticipated relationships among the position, current and 
force variables but appear to disagree with other aspects of the present theory. Figure 5 
shows the vertical force as a function of y/c for several values of current. The force tends 
to increase roughly as the inverse square of the gap. The magnitudes of the forces, 
however, are considerably lower than those predicted either by the linear finite element 
theory or by the traditional theory based on assumption of uniform gaps, and the ratio 
between measured and predicted forces is not constant. Figure 6 is a comparison of the 
measured forces with those predicted by the finite element calculation. The results indicate 
that at large gap and/or small current the ratio between the measured an ! predicted forces is 
about 1 .5, but at smaller gaps and/or higher currents this ratio increases, eventually 
exceeding 2.0 for all the three values of current that are plotted. 

Several mechanisms may be operating to cause these discrepancies, including flux 
leakage, non-uniform permeability of the materials and magnetic saturation. Some part of 
the disagreement is likely the result of measurement errors, but the differences appear to be 
significant even after allowing for reasonable experimental error. These disagreements 
reinforce the need for additional work, already planned, on force calculation. 
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Fy(N) 



y/c 


Figure 6. Measured forces and forces predicted by linear FEM calculation. 



The linear finite element theory predicts the existence of forces from a magnet that are 
normal to its axis of symmetry when the rotor is displaced from this symmetry axis, but the 
forces that are measured are considerably stronger than those predicted by calculation. 
Figure 7 shows the x component of force when the rotor is placed as closely as possible on 
the y-axis. The normal force appears to be somewhat stronger at higher current levels but 
all these forces are small, on the order of 5 % or less of the principal force, so it is difficult 
to attribute much significance to this ratio in view of the experimental uncertainty. At 
higher values of x/c, however, the normal force becomes much more significant. Figures 8 
through 1 1 show the vertical and horizontal components of force when the x/c value is 0.24 
or 0.45, and Figure 12 shows the value of the x force as a function of position for several 
values of x/c while the current is held constant at 1 .0 A. In general it appears that the 
normal force increases significantly with increasing x/c, and at x/c = 0.24 and 0.45 the 
horizontal force is about 10 % of the principal force. Theory predicts a ratio of about 2 % 

to 4 %. 

Figure 13 indicates that within measurement uncertainty there are not significant 
differences in the y-components of force at the three different values of x/c. This is in 
general agreement with the theoretical predictions, which show some increase in the y- 
component of force as the x-eccentricity is increased. 



Fx (N) 



y/c 


Figure 7. Horizontal force from lower vertical magnet at different values of current. 






Fy (N) 



Figure 10. Vertical force from lower vertical magnet when x/c - 0.45. 
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Figure 1 1. Horizontal force from lower vertical magnet when x/c = 0.45. 
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x/c~0.24 


x/c*0.45 


Lower vertical magnet 
i - 1.0 A 



Figure 12. Horizontal force at different x positions with 1.0 A current 





Figure 1 3. Vertical force at different x positions with 1.0 A current. 
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4. Conclusion 


Numerical calculations and experimental measurements of forces in a magnetic journal 
bearing actuator are presented. 

In summary, the general trends of the measured principal forces agree with the 
predictions of the theory while the magnitudes of forces are somewhat smaller than those 
predicted. The measured forces in the normal direction appear to be significantly larger 
than those predicted by theory when the rotor has an x eccentricity. The accuracy of the 
measured results has not been firmly established, and these conclusions should be 
reexamined as later measurements are made. It appears, however, that these effects will be 
significant even after considering experimental uncertainty, and both of these phenomena 

warrant further study. , 

Additional work is planned that will include force calculations considering the finite 
permeability of metals, and further force measurements using other materials. A new 
apparatus is being constructed that will allow more precise position control and force 
measurement. 
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Abstract 


A 12-state lumped-element model is presented for a flexible rotor supported by two 
attractive force electromagnetic journal bearings. The rotor is modelled as a rigid disk with 
radial mass unbalance mounted on a flexible, massless shaft with internal damping (Jeffcott 
rotor). The disk is offset axially from the midspan of the shaft. Bearing dynamics in each 
radial direction are modelled as a parallel combination of a negative (unstable) spring and a 
linear current-to-force actuator. The model includes translation and rotation of the rigid 
mass and the first and second bending models of the flexible shaft, and is unique in that it 
simultaneously includes internal shaft damping, gyroscopic effects, and the unstable nature 
of the attractive force magnetic bearings. 

The model is used to analyze the dependence of the system transmission zeros and 
open-loop poles on system parameters. The dominant open-loop poles occur in 
stable/unstable pairs with bandwidth dependent on the ratios of bearing (unstable) 
stiffnesses to rotor mass and damping dependent on the shaft spin rate. The zeros occur m 
complex conjugate pairs with bandwidth dependent on the ratios of shaft stiffnesses to 
rotor mass and damping dependent on the shaft spin rate. Some of the transmission zeros 
are non-minimum phase when the spin rate exceeds the shaft critical speed. 

The transmission zeros and open-loop poles impact the design of magnetic bearing 
control systems. The minimum loop cross-over frequency of the closed-loop system is the 
speed of the unstable open-loop poles. And for super critical shaft spin rates, the presence 
of non-minimum phase zeros limits the disturbance rejection achievable at frequencies near 
or above the shaft critical speed. Since non-minimum phase transmission zeros can only be 
changed by changing the system inputs and/or outputs, closed-loop performance will be 
limited for super-critical spin rates unless additional force or torque actuators are added. 

This paper reports work performed as part of the authors masters thesis 
[McCallum 1988], which was completed in the winter of '87 and spring of '88. This work 
was funded under Draper Lab I.R.&D. 
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Motivation 

The use of modem control system design methodologies for 
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variation of rotor dynamics with spin rate 
unstable nature of attractive-force magnetic bearings 


The application of modem, multi-input multi-output (MIMO) control system design 
methodologies to the control of flexible rotors is the focus of current research. These 
methodologies include: linear quadratic gaussian (LQG) control; LQG/LTR (Linear 
Quadratic Gaussian/Loop Transfer Recovery); and H-infinity optimal control. 

In many magnetic bearing applications, such as aircraft engines, the suspended 
rotor spins at speeds exceeding the first and second critical speeds. The application of 
modem design techniques to the control of these high speed rotors requires a linear system 
model that is simple yet includes all dynamic effects that arc important at high speeds. 

When rotors spin above their critical speed(s), rotor flexible body modes can be 
excited. When the rotor flexes, internal damping can serve as a mechanism of instability 
[Crandall 1980, Johnson 1986]. Thus, the model used for control design should include 
rotor flexibility (first and second bending modes) and the effects of internal shaft damping. 
Angular dynamics are also important at high spin rates since gyroscopic coupling between 
input axes can be large for high speed rotors. Gyroscopic effects are particularly important 
for jet engines, where the rotor-bearing system can be subject to relatively large rotation 
rates. Coupling between angular dynamics and translational dynamics can be significant, 
and should be included in the model. Finally, the unstable nature of attractive-force 
magnetic bearings should also be included. 

This paper presents a model that simultaneously includes first and second mode 
bending, gyroscopic effects, the effects of internal shaft damping, and the unstable nature 
of attractive-force magnetic bearings. The model is used to determine system open-loop 
poles and transmission zeros and their dependence on spin rate. The paper concludes with 
a brief discussion of the implications for control system design. 

The model presented here is intended to bridge the gap between models available in 
the literature when the author's thesis research was performed and models used for detailed 
rotor dynamic analysis. In practice, the model presented here should be augmented with 
finite element or experimental analyses in a magnetic bearing control system development 

effort. 
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The remainder of this paper is composed of three sections and two appendices. A 
description of the lumped-element model is presented first. A discussion of the open-loop 
eigenvalues and transmission zeros, and their variation with shaft spin rate, is then 
presented. The paper concludes with a brief discussion of the implications of the system 

dynamics for control. 

Linearized equations of motion for the model are presented in the appendix. A 
complete derivation of these equations, and a more thorough analysis of the dynamics, can 
be found in [McCallum 1988]. A bibliography also appears in an appendix. 



Model Description 
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Parameters used in the modelling effort arc from a testbed designed at Draper Lab in 
1987-1988. A testbed similar to the one shown above is now under construction. 


The testbed consists of a flexible rotor supported at its ends by two attractive-force, 
permanent magnet biased, electro-magnetic journal bearings. Inputs to the testbed arc the 
bearing currents in each axis. Measurements of the shaft-end positions are used for 
control. 


Although model parameters are from the Draper testbed, the model and results 
presented here can be generalized to other actively-controlled rotors. 

The model presented here considers only radial dynamics (radial translation and 
rotation about shaft radii). Axial dynamics decouple from the radial dynamics to first 


order. 



Model Description (cont.) 
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The rotor is modelled as a rigid flywheel with radial mass unbalance mounted on a 
flexible, massless shaft with internal damping. In the model, the flywheel is allowed to 
translate in radial directions and rotate about its radii. Each end of the shaft is modelled by 
a rigid rod of appropriate length, with a parallel combination of a spring and a damper 
acting between the rod end and the shaft wall in each of two perpendicular radial directions. 
The shaft ends are allowed to move (independent of flywheel and of each other) in shaft 
radial directions. 


As desired, this rotor model includes gyroscopic effects, first and second mode 
bending, the effects of static mass unbalance, and internal shaft damping. Note that, since 
the rigid rods are not necessarily of equal length, the effects of coupling between angular 
and translational dynamics are also included. 

Assumptions used in modelling the rotor are - 

• the rigid flywheel is dynamically balanced 

• the shaft ends are axisymmetric, so that stiffness and damping are the same in all 
radial directions 

• shaft stress is linearly dependent on strain-rate (linear internal shaft damping) 

• axial dynamics decouple from the radial dynamics 


Dynamics of each bearing in each radial direction are modelled by a parallel 
combination of an unstable spring and a linear current to force actuator. Losses and high 
frequency roll-off are ignored in the bearing model. Past magnetic bearing designs have 
demonstrated input current to output force frequency responses that are flat to high 
frequencies [Traxler 1984, Ulbrich 1984, Maslen, 1988]. This bearing model has 
appeared often in the literature [ex. Downer 1986]. 

System inputs are the currents in coils for each axis of the two bearings. It is 
assumed that (only) shaft end positions at the bearing are available. The measurements and 
bearing forces are assumed to be in the same radial plane of the shaft - the effects of 
noncollocation are not considered here. 
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Model Description (cont.) 



The figure above shows an end view of the model. Note that the shaft equivalent 
springs and dampers rotate with the shaft. This model is sometimes called the Jeffcott rotor 
in the literature [Johnson 1986]. Rotation of the dampers with the shaft is important for 
modelling shaft whirl modes, since the equivalent damping of forward and backward whirl 
modes is different for non-zero spin rates. In fact, the combination of shaft spin and 
internal damping can add energy to the system for supercritical spin rates. 



Linearized equations of motion are derived in detail in IMcCallum 1988], and are 
presented in the appendix of this paper. The model that results for this system has 12 
states, and predicts dynamics that are a function of shaft spin rate. 


The linearized equations of motion were used to analyze the dependence of system 
eigenvalues (poles) and transmission zeros on system parameters (stiffness, damping, spin 
rate). Only the dependence on spin rate is discussed in this paper. 

The model predicts 12 eigenvalues. 4 pair of eigenvalues correspond to rigid body 
motion of the suspended rotor. 

The first 4 rigid body poles correspond to angular motion of the flywheel. These 
poles occur in 2 stable/unstable pairs. Since the shaft is much stiffer then the bearings, the 
speed of these poles is dependent on the ratio of the equivalent rotational stiffness of the 
bearings to the radial moment of inertia of the flywheel. For zero spin rate, these poles he 
on the real axis. As spin rate increases, the imaginary component of these poles increases 

due to gyroscopic coupling. 

4 poles correspond 10 translational motion of the flywheel. These poles occur in 2 
stable/unstable pairs, each pair corresponding to displacement of the flywheel in the 
horizontal or vertical direction. The speed of these poles is dependent on the ratio of the 
bearing (unstable) stiffness to the rotor mass. For zero spin rate, these poles lie on the real 
axis. As the spin rate increases, the translational poles have a small imaginary component 
due to coupling with the angular rotational mode. If the flywheel were in the shaft center, 
the translational motion poles would be independent of spin rate. 

As the flywheel moves away from the center, coupling between the angular and 
translational open-loop modes increases. In these cases, all rigid body modes will be a 
combination of angular and translational motion. 



Analysis - open-loop eigenvalues (cont.) 

2 eigenvalue pairs correspond to shaft-end bending. 



REAL 


The 4 remaining poles correspond to shaft-end bending. These poles occur in 
pairs, each pair corresponding to bending of one shaft end. The speed of these poles is the 
ratio of the equivalent shaft end stiffness to the shaft end damping. For zero spin rate, 
these poles lie on the real axis. As spin rate increases, the shaft-end bending poles have an 
imaginary component equal to the shaft spin rate. This is because the projection of the 
shaft-end deflection in the lab frame changes as the shaft rotates. 
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Transmission zeros occur when, for some initial condition, a natural motion of the 
system exists for which the output vector is zero for a non-zero input vector. 

For our model, the outputs arc the shaft-end positions and the inputs arc the bearing 
currents. For a transmission zero to exist, the shaft-ends must be stationary. And if the 
shaft ends are stationary, bearing current is directly proportional to bearing force. Thus, 
for this model a transmission zero occurs if and only if there is a natural motion of the 
system where the shaft-ends remain fixed for non-zero bearing forces. 

These natural motions correspond to motions the rotor would have if mounted in 
perfectly rigid bearings. Thus, the frequencies and (part of) the directions of transmission 
zeros are the same as the eigenvalues and eigenvectors of the rotor with both ends fixed. 

The extensive literature on the dynamics of rotors in stiff bearings can be used to 
predict transmission zeros. Any rotor that is unstable in stiff bearings will have non- 
minimum phase transmission zeros. 



Analysis - transmission zeros (cont.) 
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The model predicts 4 pair of transmission zeros. The first 2 pair correspond to 
translational whirl (first-mode bending). 1 pair corresponds to forward translational whirl. 
1 pair corresponds to backward translational whirl. Forward whirl is defined as a circular 
rotor motion in the same direction as the rotor spin. Backward whirl is defined as circular 
rotor motion in the direction opposite the rotor spin. Forward and backward translational 

whirl are illustrated in above. 

The second 2 pair of transmission zeros correspond to rotational whirl (second- 
mode bending). Again, 1 pair corresponds to forward whirl, the other pair correspond to 
backward whirl. For the rotational whirl case, forward whirl corresponds to motions 
where the ends of the equivalent rigid rods move in circles in the same direction as the shaft 
spin. Backward whirl is the opposite motion. 



Analysis - transmission zeros (cont.) 
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The above figure shows the transmission zeros as a function of spin rate. 

The translational motion transmission zeros have a frequency equal to the first 
critical speed. The frequency of these zeros is only a weak function of rotational speed. 
However, the damping of these zeros is a strong function of spin rate. Damping of the 
backward whirl zeros increases with spin rate, while the damping of the forward whirl 
zeros decreases with spin rate. The forward translational whirl zeros become non- 
minimum phase for all spin rates exceeding the first critical speed (the translational 
resonance frequency). 

At zero speed, the rotational motion transmission zeros has a frequency equal to the 
second critical speed. The frequency of these zeros is a strong function of spin rate, with 
forward whirl zeros asymptotically approaching the gyroscopic whirl frequency. The 
damping of the forward whirl zeros is a decreasing function of rotational speed - for spin 
rates exceeding the second critical speed, the forward whirl zeros are non-minimum phase. 
The damping of the backward rotational whirl zeros increases with spin rate, with the 
frequency of these zeros asymptotically approaching zero. 


Conclusions - implications for control 
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that cannot be overcome through feedback. 


The open-loop system dynamics impact the design and achievable performance of 
magnetic bearing control systems. 


The minimum loop cross-over frequency of the closed-loop system is the speed of 
the unstable open-loop poles associated with rigid body motion. This result has appeared 
in the literature [Groom 1979, Downer 1986]. 

For spin rates at or above the rotors critical speed, the transmission zeros become 
non-minimum phase. As a result, the achievable disturbance rejection of the closed-loop 
system is limited at frequencies near or above the shaft critical speeds (the frequencies of 
the zeros). In addition, when the transmission zeros become non-minimum phase, each 
complex conjugate pair contributes 180 deg. of phase lag rather than phase lead. This 
transition explains stability problems encountered when spinning a magnetically suspended 
rotor through its critical speeds (as for systems with notch filters). 

Non-minimum phase zeros cannot be changed through active control. Canceling 
right-half plane zeros with compensator poles results in a system that it is internally 
unstable. Zeros can only be changed by changing the system inputs or ouputs, i.e. by 
adding additional actuators or sensors. The limitations imposed by non-minimum phase 
zeros are physical limitations that cannot be overcome through feedback. 



Appendix A - Equations of Motion 


This appendix presents linearized equations of motion for a magnetic bearing - 
suspended rotor system. Included are a summary of the assumptions used in deriving the 
equations, descriptions of reference frames, and a list of symbols. A detailed derivation of 
the equations of motion can be found in [McCallum 1988]. 

A.l Summary of assumptions 

The magnetically-suspended rotor is modelled as a rigid flywheel with radial mass 
unbalance mounted on a flexible, massless shaft with internal dampi ig. In the model, the 
flywheel is allowed to translate in radial directions and rotate about its radii. The shaft ends 
are allowed to move (independent of the flywheel and of each other) in radial directions. 
The mass is not assumed to be at the midspan of the shaft, so that the forces at the shaft 
ends are not equal. The rotor model includes: gyroscopic effects; the first and second 
bending modes; the effects of internal shaft damping; and differences in bearing loads. The 
rotor model excludes the effects of axial torque and axial translation. Further, the shaft is 
assumed to be axially symmetric so that the shaft stiffness and damping are the same in all 
radial directions. Also, the rotor is assumed to be configured such that the spin axis is 
aligned with the axial principal axis of the rotor (i.e., the rotor is dynamically balanced) and 
the radial moment of inertia is the same in all directions. 

The magnetic bearings are modelled as a parallel combination of an unstable spring 
and a linear current- force actuator in each direction 1 , i.e., 

Fyi = KyiYi + Ki Si Y i (A. 1.1) 

Bearing losses and high frequency dynamics are ignored. 

It is assumed that (only) measurements of the shaft end positions are available. 
Sensor dynamics are excluded. In addition, it is assumed that the bearing forces and the 
sensor measurements are at the same points on the shaft - non-colocation effects are 
ignored. Justifications for the assumptions listed here can be found in [McCallum 1988]. 


'This model appears often in the literature. 
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Reference frames 


A. 2 


The system state is described by the position and velocity of the flywheel's 
geometric center, the positions of the shaft ends, and rotations of the flywheel about its 
radii. Positions and velocities are described in the cartesian reference frame XYZ shown in 
Figure A.l. This frame is fixed to "ground" with its origin at the force-free equilibrium 
position of the flywheel's geometric center. The X-axis is horizontal and coincident with 
the centerline of the shaft when the system is in force-free equilibrium. The Y-axis of the 
frame is vertical. The Z-axis is horizontal and coincident with a radius of the flywheel at 
equilbrium. Figure A.l also shows the frames xyz and 123, which are used in [McCallum 
1988] to derive the equations of motion. The rotational motions of the flywheel are 
described in terms of the gyroscopic coordinates illustrated in Figures A.2 and A. 3. 



Figure A. 1 - Reference frames 
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Figure A. 3 - Gyroscopic coordinates 


A. 3 Equations of Motion 

The system states (x) are chosen as the flywheel translational and rotational 
positions and velocities and the positions of the shaft ends 


x = Y r Y r Z, Z, 5, ^ Y, Z, Y 2 Z 2 


(A. 3.1) 


r»2f> 


(A.3.2) 


System inputs are the changes in bearing coil currents from their equilibrium values 
u = [ iyi izi *Y2 1Z2 ] T 

System outputs are the positions of the shaft ends 

y = [Y t Z! y 2 z 2 ]T 

Linearized equations of motion for the rotor-bearing system can be written in sta 
form as: 

x = A((o) x + B u + d 
y=Cx+Du+n 


where 
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[e-space 


(A. 3.4) 
(A.3.5) 


(A. 3.6) 
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(A. 3.9) 


D = [0(4x4)] 


(A.3.10) 


Note that, in this representation, radial mass unbalance is modelled as an input force 
disturbance. The effects of dynamic balance could be included similarly. The effects of 
measurement center offset could be added to this model as additive measurement noise. 


bi 

d2 

d$, 

d^ 2 

dYj 

dYr 

dZj 

dZr 

la 

Ir 

>Yi 

*Zi 

k i 

Kyi 

Kzi 

li 

M 

n 

u 

u 

to 


List of Symbols 

= equivalent damping of shaft end i 
= vector of input disturbances 

= torque disturbance acting on the rotor in the -direction 
= torque disturbance acting on the rotor in the ^-direction 
= force disturbance acting on shaft end i in the Y-direction 
= force disturbance acting on the rotor in the Y-direction 
= force disturbance acting on shaft end i in the Z-direction 
= force disturbance acting on the rotor in the Z-direction 
= axial moment of inertia 
= radial moment of inertia 
= current in the Y-coil(s) at shaft end i 
= current in the Z-coil(s) at shaft end i 
= equivalent stiffness of shaft end i 

= equivalent (unstable) stiffness of bearing at shaft end i in the Y-direction 
s equivalent (unstable) stiffness of bearing at shaft end i in the Z-direction 
= equivalent length of shaft end i 
= rotor mass 
= sensor noise 
= control inputs 

= radial mass unbalance distance 
= shaft spin rate 


2 Bold denotes a vector. 
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x = state vector 

= gyroscopic coordinate, approximately equal to rotation about the +Y-axis 
%2 = gyroscopic coordinate, approximately equal to rotation about the +Z-axis 

y = system outputs 

Yj = displacement of shaft end i in the Y-direction 

Y r = displacement of the rigid rotor in the Y-direction 
Zj = displacement of shaft end i in the Z-direction 

Zj = rotor displacement in the Z-direction 

t = time 

0q = orientation of the radial unbalance vector with respect to the +1 axis 
0, <)), \|/ = Euler angles describing rotor orientation 
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ABSTRACT 

The prediction of critical speeds and forced response of 
active magnetic bearing turbomachinery is of great interest due 
to the increased use of this new and promising technology. 
Calculating the system undamped critical speeds and forced 
response is important to all those who are involved in the design 
of the active magnetic bearing system. This paper is the first 
part of a two part paper which presents the theory and results of 
an investigation into the influence of sensor location on the 
undamped critical speeds and forced response of the rotor and 
bearing system. 

Part I of this paper concentrates on an extended Jeffcott 
model which was used as an approximate solution to a more 
accurate transfer matrix procedure. Theory behind a 
two-degree-of-f reedom extended Jeffcott model will be presented. 
Results of the natural frequency calculation will be shown 
followed by the results of the forced response calculation. The 
system response was predicted for two types of forcing. A 
constant magnitude excitation with a wide frequency variation was 
applied at the bearings as one forcing function. The normal 
unbalance force at midspan was the second source of excitation. 
The results of this extended Jeffcott solution gives useful 
design guidance for the influence of the first and third modes of 
a symmetric rotor system. 

NOMENCLATURE 

A shaft relative motion max amplitude for 1st mode (cm) 
a mass eccentricity of imbalance (cm) 

B shaft relative motion max amplitude for 3rd mode (cm) 

C ratio of bearing damping to shaft damping (dim) 

damping of AMB (N-s/cm) 
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C 2 damping of flexible shaft (N-s/cm) 

F <p constant magnitude force applied to journal mass (N) 

1 square root of -1, complex variable (dim) 

K ratio of bearing stiffness to shaft stiffness (dim) 

Kj stiffness of AMB (N/cm) 

K 2 stiffness of flexible shaft (N/cm) 

Kfi row 1, column 1 of stiffness matrix in Jeffcott solution 

K 12 row x > column 2 of stiffness matrix in Jeffcott solution 

K 21 row 2, column 1 of stiffness matrix in Jeffcott solution 

K 22 row 2, column 2 of stiffness matrix in Jeffcott solution 

L bearing span (cm) 

M ratio of bearing journal mass to rotor midspan mass (dim) 
M-^ equivalent bearing journal mass (kg) 

M 2 equivalent rotor midspan mass (kg) 

R shaft absolute displacement (cm) 

X shaft maximum displacement in the X-direction (cm) 
r x shaft deflection at bearing journal location (cm) 
r 2 shaft deflection at midspan mass (cm) 
r s shaft deflection at AMB sensor location (cm) 

2 axial distance along rotor (cm) 
z s axial distance to AMB sensor (cm) 

a sensor relative position to midspan (dim) 

0 normalized shaft motion at sensor location (dim) 

u> angular velocity of shaft (rad/s) 

natural frequency normalized to rigid bearing 
critical speed (dim) 
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INTRODUCTION 


The evaluation of critical speeds and forced response for 
turbomachinery with fluid-film and antifriction bearings is now 
standard design practice for many manufacturers. The standard 
transfer matrix solution technique (Myklestad, 1944; Prohl, 1945) 
is the current industry standard for evaluation of rotor response 
and undamped critical speeds. More recently interest in improved 
forced response and stability of high pressure compressors and 
pumps have forced designers to consider active magnetic bearings 
(AMBs) for either retrofit or new machinery application. The 
initial application of magnetic bearings to centrifugal 
compressor was evaluated using standard critical speed codes 
without consideration for sensor location (Hustak et al . , 1987; 
Schoeneck and Hustak, 1987) . The comparisons of predicted 
response and critical speed placement to actual test and field 
results (Hustak et al., 1987; Schoeneck and Hustak, 1987; Kirk 
et al., 1988) have drawn attention to possible improvements in 
the analytical representation of the magnetic bearings. 

This paper is the first of a two part paper which presents 
an evaluation of the effect of sensor location on the predicted 
undamped critical speeds and forced response of turbomachinery. 
This paper concentrates on the solution of a 

two-degree-of-f reedom model developed by extending the original 
Jeffcott model to include bearing stiffness and damping, journal 
mass and accounting for non-colocation of bearing and sensor. 

The second paper will discuss the evaluation of a modified 
transfer matrix solution and will present results of a typical 
rotor bearing system analysis. 

The extended Jeffcott model will be considered to have 
sensors either inboard or outboard of the bearing centerline. 

The system response is calculated for two different types of 
forcing functions. The first excitation force is the usual 
unbalance located at the midspan mass. The second is a constant 
magnitude excitation applied at the journal mass while the 
excitation frequency is varied. The second type of rotor 
excitation is available in an actual active magnetic bearing and 
rotor system. 


PRINCIPLE OF ACTIVE MAGNETIC BEARING OPERATION 

The AMB is composed of two major mechanical parts, the rotor 
and the stator. Both are made of ferromagnetic laminations. The 
rotor laminations are placed on the machine shaft at the selected 
journal location. The stator laminations are slotted and include 
windings to provide the magnetic levitation and position control. 
For each degree of freedom, two electromagnets are required since 
they operate by attraction only. Figure 1 shows the stator 



lamination construction of a radial bearing and sensor with the 
rotor lamination sleeve in the background. 

Rotor position is monitored by sensors and this signal is 
compared to a nominal reference signal with a closed loop 
controller which supplies a command signal to the power 
amplifiers. These amplifiers provide power to the electromagnets 
to resist rotor movement from the nominal position. The design 
of the control loop gives the option to select the effective 
bearing stiffness and damping. The details of this design 
procedure are not the subject of this paper but the values of 
stiffness and damping must be carefully selected to give the 
rotor system the desired optimum dynamic response and stability. 

Before power is applied to the bearings, the rotor is 
supported on two auxiliary ball bearings located in close 
proximity to the AMB. The clearance between the rotor and the 
inner race of the ball bearing is selected to prevent rotor 
contact with the AMB pole pieces or the internal seals of the 
compressor while the rotor is at rest or during an emergency 
shutdown. When power is applied to the electronic controls, the 
electromagnets levitate the rotor in the magnetic field and 
rotation by the driving source, such as a motor or turbine, can 
be started. The sensors and control system regulate the strength 
and direction of the magnetic fields to maintain exact rotor 
position by continually adjusting to the changing forces on the 
rotor. Should both the main and redundant features of the AMB 
fail simultaneously, the auxiliary bearing and rotor system are 
designed to permit safe deceleration. 

When the turbomachine is running the rotor shaft may take a 
dynamic mode shape such that the displacement at the sensor 
location may not be the same as the magnetic bearing centerline 
displacement. The command signal is taken from the sensor 
location but the actuator applies the force through the coil such 
that the average force acts at the bearing centerline. This 
variation in command signal and actuator location is unique to 
the active magnetic bearings and may be used to the advantage of 
the designer to help place critical speeds. The performance of 
the AMB supported machinery may be more accurately predicted if 
proper account is taken for sensor location. This requires a 
modified, iterative solution strategy for current standard 
state-of-the-art computer codes for critical speeds, forced 
response, and stability. To initially evaluate the influence of 
the sensor placement, a modified extended Jeffcott rotor model 
will be developed with an assumed deformation to study the 
sensitivity of rotor bending modes and response to sensor 
location . 
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AMB EXTENDED JEFFCOTT MODEL 


The original rotor model developed by H. H. Jeffcott consists 
of a single mass on a flexible shaft supported by rigid bearings. 
Kirk and Gunter (1972) modified this model to study the effect of 
support flexibility and damping on the synchronous response of 
the single mass flexible rotor. This paper extends the original 
Jeffcott model by assuming the existence of AMB supports. The 
extended model adds journal mass, bearing stiffness and damping 
at bearing locations, and assumes rigid bearing pedestals. The 
AMB extended Jeffcott model is shown in Figure 2. 

To develop the extended Jeffcott model the disk mass plus the 
two center quarters of the shaft mass are lumped at midspan, M 2 . 
The journal and shaft end quarter masses are lumped at bearing 
locations and modeled as M^/2. The model is assumed to be 
symmetric; therefore, it can be simplified to a 
two-degree-of-f reedom (2D0F) system as shown in Figure 3. An 
unbalance force is shown at M^, and a constant magnitude 
excitation force is shown acting on M^ . 

The equations of motion (EOM) for the 2D0F system are written 
as follows; 


M 2 r 2 = M 2 w 2 ae iwt - C 2 (r 2 -r 1 ) - K 2 (r 2 ~r 1 ) [1] 

M iri = F*e iwt + C 2 (r 2 - ri ) + K 2 (r 2 - ri ) - C x r s - K x r s [2] 

In equations [1] and [2] the deflections at M^ and M 2 are defined 
as r^ and r 2 respectively. The deflection at the AMB sensor 
location is defined as r s . It is indicated by the EOM that the 
bearing forces are proportional to the sensor location 
deflection - not the bearing location deflection, as would be the 
case with conventional fluid-film or antifriction bearings. 

The sensor location deflection is calculated after assuming 
mode shapes of a half-period of a sine wave. These mode shapes, 
modeling the first and third modes, are shown in Figure 4. Using 
Figure 4, the equation for the sensor location deflection is 
written as follows: 

r s = r i + (r 2 - r^) sin(7ta/2) [3] 

where , 

a = z ? /(L/2) . [4] 

Equation [4] defines the value a as the ratio between the sensor 
offset and the shaft half-span. 
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After substitution of equation [3], equations [1] and [2] 
can be written in matrix form as follows: 


'm-l o' 

I~ r l 


[Cl 01 

" r l 


fKi 01 

~ r l 


~F<J>e iwt 

M 2. 

L r 2. 

+ 

1 

CM 

u 

0 

1 

_ r 2. 

+ 

CN 

0 

I 

_ r 2- 


M 2 u 2 ae iwt 


By assuming a solution of r = Re iwt , the matrix equation [5] 
can be written in the following form: 


K 11 

Ki 2 ' 

R 1 

= 

F <p 

[ 6 ] 

_ K 21 

K 2 2 . 

R 2 . 


M 2 u 2 . 



where , 

Kn = ( K-l( 1-0) + K 2 - Mjw 2 ) + iw( C 2 + C-^l-/?)) 
k 12 = ( K l£ " k 2 ) + iw ( c l& “ C 2 ) 

^21 ~ ™ 

k 22 = ( k 2 _ m 2 w 2 ) + iuC 2 
(3 = sin na/2. 


INFLUENCE OF SENSOR LOCATION ON UNDAMPED NATURAL FREQUENCY 

The influence of sensor location on the critical speeds of 
the AMB rotor system is initially investigated by calculating the 
natural frequencies of the extended Jeffcott Model. . The sensor 
location is varied inboard and outboard of the bearing centerline 
by as much as 20% of half-span. 

Eliminating damping terms, and C 2 , from equation [6] , and 
solving for the determinant of the resulting stiffness matrix, 
results in the following equation for the natural frequencies: 


Mi M 2 


in the following non-dimens onal 

form: 


w 


4 - (1 - 0 ) — + 


K x (M!+M 2 ) 


K-)W Z + 


M x M-l M 2 


Equation [7] can be written 


«2 - 


[ 8 ] 


K K 

(1 (1 + K ( 1 - 0 ) )Q + - = 0 

M M 


where , 


M = M 1 /M 2 


K = K 1 /K 2 
0 = w 2 /(K 2 /M 2 ) 


Results from Natural Frequency Analysis of AMB Extended Jeffcott 
Model 


To show how sensor location influences the first and third 
natural frequencies of various geometries of the extended 
Jeffcott model, the solution of equation [8] was graphed for 
alpha values ranging from -2.0 to 2.0. This exemplifies sensor 
separations of 20% of half-span both inboard and outboard of the 
bearing centerline. The results are shown for mass ratios, M, of 
1.0 (Figure 5), and 0.25 (Figure 6). The stiffness ratio, K, 
varies from 0.1 up to 10 in each analysis. 

The results are similar for both mass ratios. The 
sensitivity to non-colocation of bearing and sensor is increased 
in two different situations. An increase in sensitivity occurs 
as stiffness ratios increase. This is attributed to the fact 
that when the bearing stiffness increases relative to the shaft 
stiffness there is more bending energy in the rotor. This causes 
a greater difference between bearing and sensor deflection at 
higher stiffness ratios. The sensitivity also increases as 
sensor-bearing separation increases. This also results in 
greater differences between bearing and sensor deflections. 

The direction in which the criticals move depends on whether 
the sensor, at an inboard or outboard location, gives more or 
less response than the normal bearing centerline. For the first 
mode, the inboard sensors have a greater deflection than the 
bearings; therefore, the critical increases due to higher bearing 
forces. Outboard sensors detect less deflection than at the 
bearing, thus decreasing bearing forces and lowering the critical 
frequency. The opposite occurs at the third mode. The inboard 
sensors detect less deflection, while the outboard sensors detect 
more deflection than at the bearing centerline. This lowers the 
third critical frequency for inboard sensors and raises it for 
outboard sensors. 
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INFLUENCE OF SENSOR LOCATION ON FORCED RESPONSE SOLUTION 

The response is calculated for two forms of excitation 
applied independently to the AMB extended Jeffcott model 
equations. In order to calculate the response, the matrix 
equation [6] is solved for R x and R ? using Cramer's Rule (Anton, 
1984) . The solution has the following form: 


F<p K 2 2 “ k 12( M 2 aij2 ) 


K 11 K 22 ” K 12 K 21 
K i;l ( M 2 aw 2 ) - F <p K 21 

r 2 

K 11 K 22 “ K 12 K 21 


The resulting response is complex, thus there exists an amplitude 
and a phase angle associated with both R^ and R^ • The phase 
angles can be obtained from the following equations: 


e l 


cos 


-1 


real (R^) 

sqrt [ { real (R^) ) 2 + { imag(R 2 ) ) 2 ]_ 


[ 11 ] 


® 2 


cos 


-1 


real (R 2 ) 

sqrt [ { real (R 2 ) ) 2 + (imag(R 2 )} 2 ] 


[ 12 ] 


It must be noted at this point that the unbalance force and 
the constant excitation force are never applied at the same hi™ 6 
as this introduces additional complexities not accounted for by 
the Jeffcott model. Physically, the constant magnitude excitation 
force is applied in one plane only, therefore the motion of the 
masses is in one plane only. However, for computational and 
analytical simplicity, it will be assumed that the constant 
magnitude excitation force acts in two mutually perpendicular 
planes in a manner similar to the unbalance force. The correct 
solution for the constant magnitude excitation force plane, is 
then, simply the real part of the deflections R x and R 2 shown 
above. However, since the motion is assumed to be circular, the 
maximum amplitudes in the X-direction, X 1 and X 2 , will be the 
same as R-i and R 2 , for masses M x and M 2 . Similarly, the phase 
angles calculated from Rj^ and R 2 are also valid for motion in one 
plane. Thus, regardless of whether an unbalance force or a 
constant magnitude excitation force is applied, the solution 
technique and the solution itself will remain the same. 
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Results from Forced Response Analysis of AMB Extended Jeffcott 
Model 


The results that are to be shown are of the bearing journal 
response. Similar results occur for the mid-span mass response; 
therefore, they are omitted. 

The response amplitude at the bearing location is plotted 
versus shaft frequency in Figure 7. The excitation causing this 
response is due to an imbalance at mid-span resulting from an 
eccentricity of 0.076 mm at M 2 . In this case M = 1.0, K = 2.0, 
and the bearing damping is set to 0.283 N-s/mm. The shaft 
damping is assumed negligible. The value of a, being varied as 
in the critical speed solution, ranges from -0.2 to 0.2. 

In Figure 7 it can be seen that the first mode peak resonance 
frequencies increase from the colocation case, a = 0, when 
inboard sensors, a > 0, are used. The peak frequencies are shown 
to decrease with outboard sensor, a < 0, use. For the third 
mode, the peak frequencies are lower with inboard sensor use, and 
raised for outboard sensor use. These results are very 
consistent with the results shown from the influence of sensor 
location on natural frequencies. The same reasoning can be used 
to explain both sets of results. 

The results from the case using a constant magnitude 
excitation at the bearing location shows the same tendencies as 
the unbalance case. Shown in Figure 8 are the results of the 
bearing location excitation case. 


CONCLUSIONS 

This preliminary investigation into the effect of sensor 
location on the rotor dynamic performance of AMB turbomachinery 
gives very useful results. The natural frequency and forced 
response results from the AMB extended Jeffcott model could give 
the rotor-bearing system designer greater confidence in the 
proper selection of sensor location. 

From the test run of the AMB extended Jeffcott model the 
following specific conclusions were made: 

1. For inboard sensors, as the sensor is moved away from 
the bearing, the first mode critical frequency goes higher. For 
outboard sensors, as the sensor is moved away from the bearing, 
the first mode critical frequency goes lower. 

2. For inboard sensors, as the sensor is moved away from 
the bearing the third mode critical frequency decreases. However 
for outboard sensors, as; the sensor is moved away from the 
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bearing the third mode critical frequency increases. 

3 The vibrational characteristics of the two mass rotor 
system did not vary when a constant excitation force was used 
instead of a rotating unbalance. 

4. The effect of the sensor position on the critical 
frequencies was considerable, when the stiffness ratio was high 


5. Higher mass ratios led to increased sensitivity of third 
mode critical frequency to changing sensor positions. Vice 
versa, low mass rations led to increased damping effects in the 
third mode, making the rotor relatively insensitive to changing 
sensor locations. 


RECOMMENDATIONS 

The AMB extended Jeffcott model is a simple but 
very useful approximation of a much more complex rotor-bearing 
system. The following recommendations are made for 
extending the current approximate analysis: 

1. The addition of pedestal stiffness and damping should be 
included in this model. 

2. A method for accounting for sensor location should be 
included in the transfer matrix codes that are used to calculate 
forced response. 

3 This research should be extended to a stability analysis 
in order to improve prediction capability for AMB turbomachinery . 

4. Experimental test results must be generated for 
comparison and verification of the analyses developed. 
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Figure 2. The modified JefTcott model 
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Figure 4. Mode shapes for the modified Jeffcott model 
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LUENCE OF AMB SENSOR LOCATION ON UNBALANCE RESPONSE PREDICTION 
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ABSTRACT 

The need for better performance of turbomachinery with active magnetic bearings has necessitated 
a study of such systems for accurate prediction of their vibrational characteristics. This is the second 
part of a two part paper on the effect of sensor location on the forced response characteristics of 
AMB turbomachinery. This paper presents a modification of existing transfer matrix methods for 
rotor analysis, to predict the response of rotor systems with active magnetic bearings. The position 
of the magnetic bearing sensors is taken into account and the effect of changing sensor position on 
the vibrational characteristics of rotor systems is studied. The modified algorithm is validated using 
a simpler Jeffcott model described in part I of this paper. The effect of changing from a rotating 
unbalance excitation to a constant excitation in a single plane is also studied. A typical eight stage 
centrifugal compressor rotor is analysed using the modified transfer matrix code. The results for a 
two mass Jeffcott model are presented in part 1 of this paper. The results obtained by running this 
model with the transfer matrix method have been compared with the results of the Jeffcott analysis 
for purposes of verification. Also included, are plots of amplitude vs frequency fo the eight stage 
centrifugal compressor rotor. These plots will demonstrate the significant influence that sensor lo- 
cation has on the amplitude and critical frequencies of the rotor system. 

NOMENCLATURE 


B 

I 

t 

K 


A. 

A, 

M 

.!/„ 

y. t 

V„ 

& 

e 

i 

1 

m i 

xc 

xs 

o, c 

F4> 

K 

a 


Modulus of elasticity (N /mm 2 ) 

Moment of inertia of massless shaft ( mm 4 ) 

Polar moment of inertia (m 1 ) 

Transverse moment of inertia (mm 4 ) 

Stiffness ratio, k x lk 2 (dim) 

Stiffness in transfer matrix equations (N/mm) 

AMB stiffness (N/mm) 

Shaft stiffness (N/mm) 

Mass ratio, (dim) 

Moment component in transfer matrix equation (N-/?im 2 ) 

Moment component in transfer matrix equation (N -mm 2 ) 

Shear component in transfer matrix equation (N) 

Shear component in transfer matrix equation (N) 

Eccentricity component in transfer matrix equation (mm) 

Eccentricity component in transfer matrix equation (mm) 

Damping (N-scc/mm) 

Exponential constant = 2.7 1 828 18 (dim) 

Complex constant (dim) 

Length of massless shafts in transfer matrix equation (mm) 

Point masses in transfer matrix equations (Kg) 

Deflection component in transfer matrix equations (mm) 

Deflection component in transfer matrix equations (mm) 

Slope component in the transfer matrix equation (rad) 

Slope component in the transfer matrix equation (rad) 

Constant excitation force in Jeffcott model (N) 

Constant excitation force component in the transfer matrix equations (N) 
Constant excitation force component in the transfer matrix equations (N) 
Ratio of the distance between bearing centerline and sensor to half-span 
of the rotor (dim) 

Frequency of shaft excitation (rad/scc) 


560 



s 


Frequency of shaft rotation (rad/sec) 


Subscripts 

i Rotor station number in transfer matrix equations 

INTRODUCTION 

The use of active magnetic bearings in turbomachinery is a comparatively new development but 
one which has shown great promise for better control of rotating equipment. I he idea behind these 
bearings is not new however. The use of magnetic attraction to levitate the rotor shalt free of the 
bearing had been tried before, but because the system is inherently unstable unless a real-time 
control system is used, the use was not successful. The first actively controlled beanng was devel- 
oped in the 1950's. Since then the use of active magnetic bearings has gained widespread acceptance 
particularly in North America and Canada. VVeise [2] has given some examples of the varied uses 
to which active magnetic bearings have been applied. Kirk [5] lists a number of turbomachmcry 
installations where active magnetic bearings have been used. Magnetic bearings possess a number 
of advantages compared to conventional bearings. I hey give an almost unlimited control over roto 
vibrational characteristics due to adjustable stiffness and damping. Automatic balancing is possible 
by allowing the rotor to spin on its inertial axis. This leads to decreased vibrations and noise. Active 
magnetic bearings do not require lubrication, and since they are non-contact bearings, they elimi- 
nate the possibility of wear and tear of the stator and rotor surfaces. We.se [2] demons rates he 
tolerance of magnetic bearings to a wide range of temperatures and also their insensitivity to hostile 
environments. Zlotykamicn [ l] gives a good description of the various advantages of active mag- 
netic bearings. 

Most of the research in active magnetic bearings has been in the control systems used. Schweitzer 
[7] shows a method for controlling an elastic rotor so that it can be represented by a low order 
model amenable to control techniques. Williams, Keith and Allaire [6] have developed theoretical 
relationships to relate the characteristics of a controller transfer function to the stiffness and damp- 
ing properties of an active magnetic bearing. Burrows and Sahinkaya L_8] have evaluated various 
strategies for applying a magnetic bearing to control the synchronous vibration of a flexible rotor. 
Kirk et al [5] have presented results of shop tests on a high speed eight stage centnfugal compressor 
supported by active magnetic bearings along with some design recommendations. Keesee UJ has 
examined the effects of sensor position on the cnhcal frequencies of rotors with active magnetic 
bearing. This work is an extension of Kcesec's work to include sensor position effect on forced re- 
sponse vibration amplitudes, using the modified transfer matrix method. 

RESEARCH OB JECTIVE 

As stated before in the introduction, the sensors are not located at the place where the attraction 
forces are applied on the rotor shaft, but at some distance away along the axis of the shaft. Due to 
this “non-colocation” of the sensors from the bearing position, the deflection sensed by the sensors 
is not the same as the actual deflection at the bearing but differs from it by some magnitude dic- 
tated bv the mode shape of the rotor shaft. Because of this, the stiffness and damping forces of the 
active magnetic bearing depend not on the deflection at the beanng location, but on the deflection 
at the sensor location, l or such cases, the vibrational charactenstics of the rotor system is different 
from that obtained using conventional analysis programs. I he objective of this research is to take 
into account, the effect of sensor non-colocation on the vibrational charactenstics of rotors with 
active magnetic bearings. 

This research is an extension of the work done by Keesee [3] and involves the modification of an 
existing transfer matrix code to account for sensor non-colocation. But whereas Keesee s research 
was limited to studying the effect of sensor non-colocation on critical frequencies, this work also 
considers sensor non-colocation effects on forced response amplitudes. I he other objective of this 
research was to compare the vibrational characteristics of rotor systems, when they arc subject to 
unbalace excitation with circular synclironous shaft rotation, and constant excitation in one plane 
with no shaft rotation or whirling. The effect of changing mass ratios and stiffness ratios was also 

studied. 

The modification of the transfer matrix program was validated by comparing its results for the two 
mass model with the results obtained from a simple program written to specifically analyse the two 



mass model described in part I of this paper. A typical eight stage centrifugal compressor rotor 
model was also analysed using the modified transfer matrix method and its icsults were compared 
with the results of the Jellcott model to verify the trend of behaviour of ihe rotor system with 
varying sensor locations. 

THE MODIFIED TRANSFER MATRIX METHOD 


The first analytical study of flexible rotors using the transfer matrix met! od was presented by 
Myklestad and Prohl. 'I he rotor is divided into several discrete masses cal ed stations and these 
masses are joined by massless flexible shafts. The response of the system is computed by using in- 
fluence coefficients, and formulating a set of equations. The equations are solved and a final sweep 
is made to obtain the solution. J. VV. Lund analysed the equations involved in the transfer matrix 
method for the case of elliptic non-synchronous response of the rotor system and wrote a program 
using these equations to study the vibrational characteristics of rotor systems. The program was 
subsequently simplified to analyse circular synchronous response of rotor systems. This paper de- 
scribes the modification of precisely this simplified transfer matrix program written by E. J. Gunter 
Jr. and R. G. Kirk at the University of Virginia, Charlottesville, Virginia. The modification was 
done to take into account the non-colocation of sensors in active magnetic bearings. 

To understand the classical transfer matrix method, consider a typical rotor section element as- 
sumed to be composed of a point mass and a massless elastic shaft to its right. 


Consider the forces acting on the mass to formulate the equations required for equilibrium. Refer- 
ring to figure 1 
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T he solution in the Y direction can be obtained from the solution in the X direction since it is as- 
sumed that the motion of the shaft is circular. Hence equations in the Y direction are not required. 
Now consider the equations for the massless elastic shaft of station i. f rom figure 2 
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These equations can be presented in matrix form and the malnces arc called transfer matrices As 
can be seen from this matrix, a correction to account lor the shear delormation effect has also been 
included. The terms of this correction factor are explained below 

™ 1 ‘ E 17] 

a t = area of section i 

G { = shear modulus of section i 

[(7. + .6jQ(l. + dS)] ± (20- + 12-mWa- 2 ] [is] 

[6.(1. + /t)(ll. + dr 2 ) 2 ] 
fi is poissons ratio = EJ2G i — 1 

<ir ( is the diameter ratio = inner diameter /outer diameter 

MODIFICATION FOR SENSOR NON-COLOCATION 

Due to sensor non-colocation, at the station representing the bearing location, equations 
[3.1] and [3.2] are modified as follows 

= y'xlc + m i^ x ic - Ki* ic„ n - C i<° x ls,„ + FXC l + \ m i w2 

v xis = v xis + c f oX k sm + m f 0 2 x b - Kpc iSun + FXSj - ay^w 1 


[l-o 

[ 20 ] 


The bearing stiffness is multiplied by the deflection sensed at the sensor location instead of the ac- 
tual deflection at the bearing. 

ALGORITHM FOR MODIFICATION DUE TO SENSOR NON-COLOCATION 

The modification in the point matrix for the bearing station, due to the sens »r non-colocation has 
already been discussed. However a straight forward sweep of the rotor is p issible only in certain 
cases of sensor location. Upon examination, three cases ol sensor location dative to the bearing 
location can be listed. 

1. One sensor before the bearing 

2. One sensor after the bearing 

3. Two sensors, one each on either side of the bearing. 

For case 1, the sensor deflections arc saved in the sweeping process and then used at the bearing 
station. The sweeping process is straightforward. Refer figures 3. 

In this case, since the sensor comes after the bearing, the sensor deflection, are not known when 
the sweeping process reaches the bearing. Thus the sensor deflections arc assumed to be some ar- 
bitrary value. Generally, the deflections at the station belore the bearing are used as these arbitrary 
values. The sweeping process is then continued until the sensor location is reached. Here a 
comparison is made between the assumed sensor deflection and the sensor deflection calculated by 
the sweeping process. If the two quantities agree to within a certain margin of error, the sweeping 
process is continued from the sensor station onwards. If the two quantities do not lie within the 
error margin, the program iterates back to the bearing location and uses the sensor deflections cal- 
culated by the current sweeping process. These sensor deflections are used, as explained before, in 
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the bearing station point matrix calculations and the sweep process is continued. This leads to a 
series of iterations between the bearing station and the sensor station and these iterations are con- 
tinued until the sensor deflections used at the bearing station agree with the ensor deflections cal- 
culated at the sensor station by the sweeping process, i.e. convergence is obt; med. refer figures 3. 

Case 3 

With two sensors, one each before and after the bearing, the case can be split ip into two cases one 
resembling case one and the other resembling case two. Refer figures 3. W ten the sensor before 
the bearing is reached, the sensor deflections arc saved. These arc used, the irst time the bearing 
location is reached. The sweep process is then continued and the sensor aft r the bearing is dealt 
with in a manner similar to case two. 

DISCUSSION OF THE CONVERGENCE PROCESS 

To aid the process of convergence to the correct values of sensor deflcctio is, the Taylors series 
convergence technique in two variables was used. This method was the m >st suitable one since 
there is cross-coupling between the stiffness and damping terms. However, due to the very low 
magnitude derivatives involved, the convergence process fails and leads tr divergence from the 
correct solution. 

When the cross-coupling of the stiffness and damping terms was ignored, and the Secant method 
of convergence was used to converge on the sensor deflections along the two axes independently, 
the algorithm converged with diminishing oscillations. However the number of iterations required 
were more than those required, when no convergence algorithm was used. 

T hus simply using the sensor deflections obtained from the sweeping proce">, back at the bearing 
location, gave the fastest convergence. Refer table 1. 

MODIFICATION TO SEPARATE GYROSCOPIC STIFFNESS FROM TRANSVERSE 
STIFFNESS 

When the rotor is subjected to an external vibrational force assuming no ur balance to be present, 
the gyroscopic stiffness will depend only on the rotor spinning speed and n( t on the frequency of 
excitation. Considering the rotor spinning frequency to be “s”, equation > [3.3] and [3.4] are 
modified as follows 

Myic = Myic + - <»h )0 X ,c [ 21 ] 

M yis = M yis + oj{\I p - wl T )0 xis [ 22 ] 

Here w is the frequency of excitation. 

COMPARISON OF IIIE RESULTS OF THE 2 MASS ROTOR SYSTEM, FOR THE 
JFFFCOTT MODEL AND THE TRANSFER MATRIX METHOD 

The results of the 2 mass rotor system as obtained by the JefTcott Model program have already been 
shown and discussed in part I of this paper. The same rotor system data was used with the modified 
transfer matrix method program, so as to compare the results with the JefTcott model and thus 
validate the correctness of the modifications. T he JefTcott model is important, but because of its 
simplicity, its results are of limited use. Also it docs not model a complex rotor system composed 
of many disk masses and possibly different shaft cross-sections along the rotor length. Hence, it is 
the transfer matrix method that is more useful for application purposes, and the Jcffcott model will 
serve for the purpose of comparison only. 

Tables 2 and 3 give the comparison between the two programs. As can b • seen from the tables, 
there is a fairly close agreement between the results obtained from the two programs. The agreement 
in the critical frequency values is much better than that between amplitude values and again, am- 
plitude values agree better than phase angle values. This is because, the critic, 1 frequencies of a rotor 
system depends mostly on its mass and stiffness properties, both of which are accounted for in a 
similar manner in the two programs. The amplitude and phase angle values show greater disagree- 
ment due to the fact that the JefTcott code assumes a sine-wave shape for th mode shapes and this 
assumption is only an approximation of the actual mode shape. It can be si n from the tables, that 
the JeiTcott code underestimates amplitude values in most cases. Also, amplitude values show a 



higher disagreement when a constant force excitation is applied instead of an unbalance force 
excitation. 

THE EIGHT STAGE CENTRIFUGAL COMPRESSOR ROTOR SYSTEM MODEL 

To obtain a more realistic idea of the influence of sensor position, an eight stage centrifugal 
compressor rotor system was used with the transfer matrix code. I his rotor system is illustrated in 
figure 4 and its design parameters are given in table 4. 

The program was run with all the four cases of sensor positions, namely sensor colocation, inboard 
sensors, outboard sensors and two sensors on each side of the bearing, a values of -0.18, 0 and 0.18 
were used, as is indicated by the sensor locations. 

RESULTS OBTAINED FROM TIIE EIGHT STAGE COMPRESSOR ROTOR SYSTEM 

The results obtained from the eight stage compressor rotor system arc summarized below and il- 
lustrated in figures 5 to 11. 

1. The first mode critical frequency increases as the sensor is moved from the direction of the 
outboard location in the direction of the inboard location. 1'his is in agreement with the results 
obtained from the two mass rotor system. 

2. A significant difference between these plots and those of the two mass model, is the higher 
amplitudes exhibited by the compressor rotor system when the sensors arc moved inboard. 
To verify this deviation, an approximate two mass model of the compressor rotor system was 
run with the Jeffcott program. The results obtained with this approximate model are shown 
in figures 7 and 8, and show agreement with the results obtained by the transfer matrix program 
with the compressor rotor system data as input, i he reason for this behaviour will become 
clear when the mode shape shown in figures 9 is examined. 

The dcilection at the sensor location is less compared to the deflection at the bearing location. 
Due to this, a lesser stiffness and damping force is applied at the bearing location and this leads 
to higher amplitudes of the rotor system. 1'or a certain inboard sensor location, the deflection 
at the sensor location is reduced to zero and this condition will produce the largest amplitudes 
in the rotor system. The peculiar first mode shape that produces this phenomenon is similar 
to that observed in the third mode, and seems to be the result of the high bearing damping 
values. It has been observed that as sensors are moved inboard, the first critical frequency in- 
creases and the third critical frequency decreases. This is shown in figure 10. High bearing 
damping may bring the first and third critical frequencies together in such a case and thus 
produce a first mode shape similar to the third mode shape. 

3. Additionally, it was observed that when the mass ratio was increased, the maximum amplitudes 
occured with inboard sensor locations nearer to the bearing location. Refer figures 1 1. This 
can also be explained from the plots of the mode shapes. The maximum amplitude is observed 
when the deflection at the sensor location is zero. In such a case, the stiffness and damping 
forces at the bearing location are reduced to zero and the rotor system essentially exhibits 
free-free vibration. Due to this, the dcilection along the rotor longitudinal axis will depend only 
on the mass distribution of the rotor system and not on the bearing stillness and damping. A 
higher mass ratio means greater mass at the bearing location and therefore, lesser deflection in 
free-free vibration than that at midspan. In such a case, the point of zero deflection occurs 
nearer to the bearing location and a sensor placed at this point will produce the maximum 
amplitude of vibration of the rotor system. 

CONCLUSIONS 

The modification of the rotor dynamics codes to account for sensor non-colocation show a definite 
change in the vibrational characteristics when ihe sensors are moved away from the bearing lo- 
cation. The following conclusions can be drawn from this research: 

1. The first mode critical frequency increases as the sensor is moved from the outboard to the 
inboard direction. This is due to the fact that for the first mode, the sensors sense a greater 
deflection as they move inboard and away from the bearing location. This increases the effec- 
tive stiffness of the active magnetic bearing and results in higher critical frequencies. Because 
of this effect, it is possible to bypass the first critical by using the inboard sensors while starting 
the rotor and when the rotational frequency nears the first critical, switching to the outboard 


sensors. This has been suggested by Keesec [3]. The higher critical frequencies can be handled 
similarly. 

2. The amplitudes at the first critical will be higher with outboard sensors and decrease as the 
sensors are moved inboard and away from the bearing. However, this is not true in certain 
cases where the first and third mode coincide. In such a case, the amplitudes at the first critical 
increase as the sensors arc moved inboard upto a certain point and then decrease again. The 
reason for such behaviour can be traced to the presense of high damping values along with the 
condition of inboard sensors. The behaviour of the critical frequencies and amplitudes, with 
regard to changes in the sensor position, can be predicted by examining the mode shape of the 
rotor shaft at or near the critical frequencies. 

3. The results indicate a fairly close agreement between the Jeffcott model and the transfer matrix 
model. The comparison indicates greater differences in amplitude values compared to critical 
frequency values and greater deviation in phase angle values compared to amplitude values. 

4. The behaviour of the rotor system, with respect to changes in sensor location, docs not indicate 
any significant deviation when a constant force excitation is used instead of an unbalance force 
excitation. 

5. The effect of sensor location on first mode critical frequency increases with higher stiffness ra- 
tios. This has been explained in the second chapter. Higher mass ratios lead to increased 
damping effects in the third mode and hence lower the amplitudes considcrablv. 

The following recommendations can be made for future work in this area: 

1. The transfer matrix model modified for this work, does not consider the effect of pedestal 
stiffness and damping, ie it assumes a rigid foundation. The program can easily be modified to 
take this factor into account. 

2. The existing transfer matrix code can only handle circular synchronous rotation of the rotor 
system. It can be extended to analyse non-circular and non-synchronous motion of the shaft. 

3. The constant force excitation in the transfer matrix program is applied as a force function, di- 
rectly on the journal mass. The behaviour of the rotor system, when the constant force 
excitation is applied as a displacement function, and on the bearing or pedestal mass, needs to 
be investigated. 

4. fhe modified transfer matrix code assumes that no couplings are present in the rotor longi- 
tudinal cross- section. The code can be modified for the presence of a coupling, which then, 
would only transfer displacements and shears across the connection, but would not transfer the 
moments. 
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Figure I. Forces acting on the point mass 
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FREQUENCY, W (HZ) (*10 1 ) 


TYPICAL EIGHT STAGE CENTRIFUGAL AMB COMPRESSOR ROTOR, MODE SHAPES. AMPUTl 
K-22766.4 N/mm 0-91.77 Nsie/mm UNBAL-70.62 Nmm 1ST CRT AT MIDSPAN 



ROTOR AXIS ( MM ) (*10 2 ) 


INFLUENCE OF RUB PROBE LOCATION ON CRITICAL SPEED PREDICTION 
EIGHT STAGE CENTRIFUGAL COMPRESSOR. ALPHA - O.O. 0.18 
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2 MASS MODEL, UNBALANCE EXCITATION, M= 1, K = 2, ALPHA = 0.3 
FREQ = 130 HZ, CONVERGENCE TO XC, TRUE SOLUTION = 7.649442 


No. of 

T aylor's series 

% Diff. from 

Secant method 

% Diff. from 

Simple 

% Diff. from 

Iter. 

Convergence 

true Soln. 

Convergence 

true Soln. 

Iteration 

true Soln. 

I 

3.500004 

54.24 

3.500004 

54.24 

3.500004 

54.24 

2 

8.051951 

-5.26 

8.051951 

-5.26 

8.051951 

-5.26 

3 

7.737319 

-1.15 

7.737319 

-1.15 

7.737319 

-1.15 

4 

7.581234 

0.89 

7.93081 1 

-3.68 

7.628863 

0.27 

5 

7.849375 

-2.61 

7.605373 

0.576 

7.648726 

0.0094 

6 

3.625942 

52.60 

7.642738 

0.088 

7.650135 

-0.0091 

7 

69.54724 

-809.2 

7.655252 

-0.076 

7.649398 

0.00058 

8 

-959.629 

12645.1 

7.650030 

-0.0077 

7.649425 

0,00022 

9 

15109.55 

-197425 

7.649057 

0.0050 

7.649445 

-0.00004 

10 

-235777.3 

3082381 

7.649363 

0.001 

7.649442 

0.0 

11 

3681287.0 

-48124809 

7.649458 

-0.00021 

7.649442 

0.0 

12 


___ 

7.649447 

-0.000065 

7.649442 

0.0 

13 




7.649434 

0.0001 

7.649442 

0.0 

14 


_ __ 

7.649441 

0.000013 

7.649442 

0.0 

15 




7.649443 

-0.000013 

7.649442 

0.0 

16 


... 

7.649442 

0.0 

7.649442 

0.0 

17 


___ 

7.649441 

0.000013 

7.649442 

0.0 

18 


... 

7.649441 

0.000013 

7.649442 

0.0 

19 

— 

— 

7.649442 

0.0 

7.649442 

0.0 


Table 1 . Comparison of different convergence schemes 
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FIRST MODE, UNBALANCE EXCITATION, M = 1, K = 2, BEARIN(, /(MIDSPAN) 


JefTcott Code 

Transfer matrix 
Code 

% Difference 

« = -0.2 

Critical Frequency 
Amplitude 
Phase Angle 

32.167 /(32.167) 
2.285 / (4.425) 
103.48 / (102.57) 

32.167/(32.167) 
2.302 / (4.482) 
86/(85.1) 

0/(0) 

-0.74/ (-1.29) 
16.89 /(17.03) 

a = 0.0 

Critical Frequency 
Amplitude 
Phase Angle 

35.333 / (35.333) 
1.408 / (3.401) 
90.214 / (88.398) 

35.333 / (35.333) 
1.408 / (3.402) 
91.6 / (89.8) 

0 / (0) 

0 / (-0.03) 
-1.54 / (-1.59) 

a - 0.2 

Critical Frequency 
Amplitude 
Phase Angle 

40.333 /(40.333) 
0.539 / (2.281) 
92.571 / (86.391) 

40.167 / (40.000) 
0.571 / (2.310) 
104.1 / (81.8) 

0.41 / (0.83) 
-5.94 / (-1.27) 
-12.45/(5.31) 


Table 2. Comparison of the results of the 2 mass rotor system, for the JefTcott model and the transfer 
matrix method 




FIRST MODE, CONSTANT FORCE EXCITATION, M= 1, K = 2, BEARING/(MIDSP AN) 


Jcffcott Code 

Transfer matrix 
Code 

% Difference 

a - -0.2 

Critical Frequency 
Amplitude 
Phase Angle 

32.167 /(32.167) 
4.577 / (8.896) 
104.24 / (104.24) 

32.167 /(32.167) 
5.410 / (10.522) 
86.8 / (86.8) 

0 / (0) 

-18.20 / (-18.28) 
16.73 / (16.73) 

a - 0.0 

Critical Frequency 
Amplitude 
Phase Angle 

35.333 / (35.333) 
3.045 / (7.352) 
90.214 / (90.214) 

35.333 / (35.333) 
3.262 / (7.882) 
91.6/(91.6) 

0 / (0) 

-7.13/ (-7.21) 
-1.54/ (-1.54) 

<x • 0.2 

Critical Frequency 
Amplitude 
Phase Angle 

40.333 / (40.333) 
1.335 /(5.643) 
88.527 / (88.527) 

40.000 / (40.000) 
1.420 / (5.709) 
83.9 / (83.9) 

0.83 / (0.83) 
-6.37 / (-1.17) 
5.23 / (5.23) 


Table 3. Comparison of the results of the 2 mass rotor system, for the JefFcott model and the transfer 
matrix method 




ROTOR SYSTEM PROPERTY SI UNITS 


ENGLISH UNITS 


Total rotor length 1879.6 mm 74.0 in 


Distance to bearing I centerline 

304.8 mm 

12.0 in 

Distance to sensor at bearing I: 



Outboard sensor 

190.5 mm 

7.5 in 

Inboard sensor 

419.1 mm 

16.5 in 

Distance to bearing 2 centerline 

1574.8 mm 

62.0 in 

Distance to sensor at bearing 2: 



Inboard sensor 

1460.5 mm 

57.5 in 

Outboard sensor 

1 689. 1 mm 

66.5 in 

Mid- span diameter 

177.8 mm 

7.0 in 

Journal diameter 

177.8 mm 

7.0 in 

Journal length 

254.0 mm 

It 0 in 

Total rotor weight 

2.95 KN 

663.2 lb. 

Reaction at bearing I 

1.41 KN 

317.0 lb, 

Reaction at bearing 2 

1.54 KN 

346.2 lb. 


Table 4. Data for eight stage centrifugal rotor system model 
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REDUCTION in BEARING SIZE due to SUPERCONDUCTORS in MAGNETIC BEARINGS 
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Objective: Assess a Concept to Reduce the Size of Magnetic Bearings 
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Superconducting Bearings Are Smaller and, Hence, 
Suitable for Cryogenic Aerospace Applications 




o 


</> 

O) 

c 

IM 

c5 

o 

m 


o 

• mi 


<D 

O) 

O 

>. 


o 




< 


tr 

o 

o a 

C "D 3 
0) CO to 
0 ) 0 * 
0—0 

o 

O - 

-C o 

_ o '§ 

(U "2 O) 


CD 

o 


c 

o 


o 

>•* 


53 ® 

CL 5 

c O 

is 

C TJ 


cd ° 
c 


c 

CD 

E 

c 

o 

"D •=: 


2 CD 

O f 

O <= 

L. 

tC .CO 

Q_ 

CL 

13 .E 
CO ^ 

O <D 
~ xi 
o CD 

S 3 

■I ° 

_Q £3 

<D 


03 


<D 


s,*- 

1 ixi 

as s 

|8 

d Q) 

E ~ 

«1 

O) o 

— CO 

2 w 

CL 

^ E 

1 1 

5 -e 

O 3 


> 

C 

<D 

\ 

03 

CD 

i 

03 

c 

o 

o 

(D 


CD 

> 

O 

E 

CD 


T3 

CD 

O 

C 

CC 

CO 

a> 

CD 


CD 

O 

c 

CD 

X 


CO 

■# i 

xj 3 

c CO 

o >> 
a co 


<D 

03 

| $ 
as .© 

Ko 

c 3 

o o 

CO > 

« I 

w '« 

o nj 

O T3 
C O 

§>E 


_ E 

c jc 


“O 


si 

U- 

cd 

=5 CL 
^ 03 
CD ^ 

© o 

CO O 

jh 

E .b= 
c o 
o cr 

.i= CD 

> 

C o 

° £ 

03 CD 
CD > 

5 3 


CD 


O 

cd 


<d o 

> 

O CO ri 

E « o 

£ © Q. 

2 o| 

CO o CO 
d c _ 

o _c o 

-500 

CL CD _ 

O 

CD £ 

o 


CD 

$ 

O 


CO 

> — 

E ® 

03 O 
«M* C 
0 Cd 

0 E 
-c E a 

~ o c 

> t o 

3 s. a 
§ £ ** 
E 

CD 

5 3 

0-3 

CD 


O 
JO 

5 
o 

- 2> 

T3 X) 

rrf ^ 


cd 

o 

0 

CD 

E 

c 

cd 

o 


c 

5 


o 
03 

C 
O 

o c 

CO -* 

X co 
CD ~ 
CD -C 
c .E 

CD -t-j “E 
^ c 5 
rt © 


O ^ 
CD 
JD 

C 

o 

CL 
O 

O) C 
E o 

cd d) 

CD c 


CO 
"O 

<5 *0 


o . 

CD 


g 

c 

CD 
03 
O 

>N 

L. 

o 
o 

-a ^ 

.y cd — 

CL 


!; 

E 

(D 


*g 

’^3 

CD 


| H 

CO 

2 Cd 


cd 

o 

c 

g 

03 

g 

L. 

-Q 

JD 

CD 


X 

CD 


E -^ 0^3 


o 

CD 


c 

o 

o 


<0 

CD 


1 — X) CL 


cd 

cd g 

<0 CO CO 
CD >* C 

> r (1) 

> Q_ CO 


^ 03 


CD 

Q_ 

E 

CD 


58 <) 


Mechanical Technology Incorporated 


APPROACHES TO LONG-LIFE CRYOGENIC BEARINGS 
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Magnetic bearings are fundamentally limited by the magnetic, electric, and mechanical strength of 
materials used. The mechanism by which magnetic strength (flux density) limits the load capacity 
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Operating Principle of SUPERC Design Conce 
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ACTIVE SUPERCONDUCTING BEARING CONCEPT 

Superconductor 
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Effect of Current Density on Size Reduction 
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I. INTRODUCTION. 

This research is part of a NASA program to investigate the feasibility of applying 
high-T c superconducting ceramics to the development of a high-field (5-6 Tesla) magnet, 
with an ultimate objective of utilizing these materials in a liquid nitrogen cooled wind tunnel 
magnetic suspension balance. The current state of the high-T c superconducting ceramics 
inhibits their immediate application to high-field magnets for a number of reasons. Their 
superconducting properties are strongly anisotropic, with preferential critical performance 
perpendicular to the perovskite c-axis, i.e., in the a-b plane of the crystal. As a consequence, 
the bulk ceramics, which have largely random crystal orientations are not suitable for 
optimal engineering. Epitaxial films can be prepared in which the c-axis is parallel to the 
plane of the substrate, i.e., perpendicular to the substrate surface. However, since the 
design and demonstration of a high-T c superconducting high-field magnet will require 
preferential c-axis orientation over a large topological surface, many of the vacuum epitaxial 
deposition techniques are not suited for this application because of the limited surface areas 
which can be coated. 

These issues an be circumvented by advancing chemical spray pyrolysis as a means 
to deposit large area films on to a tape wire which can then be configured into the geometry 
necessary to drive the magnetic field. Chemical spray pyrolysis has a number of advantages. 
It is inexpensive. It can produce good quality films over very large areas. It is an open 
system deposition technique so continuous strips of film can be prepared, and it has already 
demonstrated feasibility in a number of electronics applications, particularly in the field of 
solar cell fabrication. 

This paper reports on the deposition of Tl-Ba-Ca-Cu-O thick films by spray 
pyrolyzing a Ba-Ca-Cu-O precursor film and diffusing thallium into the film to form the 
superconducting phase. This approach was taken to reduce exposure to thallium and its 
health and safety hazards. The Tl-Ba-Ca-Cu-O system was selected because it has very 
attractive features which make it appealing to device and manufacturing engineering. 
Tl-Ba-Ca-Cu-O will accommodate a number of superconducting phases. This attribute 
makes it very forgiving to stoichiometric fluctuations in the bulk and film. It has excellent 
thermal and chemical stability, and appears to be relatively insensitive to chemical 
impurities. Oxygen is tightly bound into the system, consequently there is no orthorhombic 
(conductor) to tetragonal (insulator) transition which would affect a component’s lifetime. 
More significantly, the thallium based superconductors appear to have harder magnetic 
properties than the other high-T c oxide ceramics. Estimates using magnetoresistance 
measurements, [1], indicate that at 77 K TEBajCaCuiOio will have an upper critical field, 
Hc 2, of 26 Tesla for applied fields parallel to the c-axis and ~ 1000 Tesla for fields oriented in 
the a-b plane. [2] 

Results to date have shown that superconducting films can be reproducibly deposited 
on 100 oriented MgO substrates. One film had a zero resistance temperature of 1 1 1.5 K. 
Furthermore, x-ray diffraction analysis of the films showed preferential c-axis orientation 
parallel to the plane of the substrate. These results have now made it possible to consider 
the manufacture of a superconducting tape wire which can be configured into a topology 
useful for high-field magnet designs. This paper will review the research which lead to the 
preparation of these films and plans for further development. 

II. EXPERIMENTAL RESULTS. 

In chemical spray pyrolysis soluble salts of the metallic components of the 
superconducting ceramic are dissolved in solution in stoichiometric amounts. In this work 
the nitrate precursors of barium, calcium, and copper, (Ba(NC>3)2, Ca(NC>3)2, C'u(N03)2). 
were used. Once thoroughly mixed this solution can be sprayed onto a surface which is hot 
enough to volatilize the transport fluid and then allow the salts to decompose, oxidize and 



form a polycrystalline ceramic. Ultrasonic techniques were used to nebulize the solution 
into droplets. The sprayed surface is thermally treated in the presence ^ *5“ 

suitable, to decompose the nitrates into their respective oxides Once a Ba-Ca-Cu-O oxide 
precursor film had been formed thallium was diffused into the film which was then thermally 
processed to produce the superconductive Tl-Ba-Ca-Cu-O phase. 

Stoichiometric control is crucial to the electromagnetic performance of high-T c 
superconducting ceramics. Thus a primary objective of the Phase I effort is to tune the 
spray process to produce Ba-Ca-Cu-O precursor thin/thick films with stolcl ? 1 °'j n ^J^ .n" : 
and uniformity. A second objective of the program is to fuse all the individual metallic ions 
into a uniform polycrystalline ceramic prior to thallium diffusion processing. Oxides of Ba 
and Ca fused into crystalline copper oxide structures are less likely to outgas during 
subsequent thermal annealing and thallium diffusion processing than would amorphous 
barium or calcium oxides loosely distributed in the film. Furthermore, a fused crystal 
structure will have better order and chemical uniformity. 

It is therefore a fundamental objective of this program to eliminate single species 
clustering or agglomerations in the films, enhance grain alignment, and to promote superior 
microstructure° This can be achieved by using a suitable chemical binder in the solution 
which favors thorough dissolution, intermixing and complexation of the different reactants 
prior to deposition. Ideally, this solvent should enhance mutual reactivity of the salts during 
the crystal fusion process. 


A. Precursor Development and Stoichiometric Control. . 

To develop stoichiometric control over the process and to avoid single component 
clustering in the films the effect of solution solvent and substrate temperature on metal ion 
incorporation, surface morphology and chemical uniformity were evaluated. The thermal 
behavior of each of the individual nitrate compounds was measured. 1 hese measurements 
were then used to evaluate the effect of a particular solvent by comparison to thermal 
analysis on a dried residue of that solution in which the nitrates were mixed in stoichiometric 

amounts. 


Since nitrate compounds are generally hygroscopic and their molecular weights can 
be altered by ambient humidity conditions, stock solutions of each nitrate were prepared. 
Stock solution concentration was measured using Atomic Absorption spectroscopy, (AAJ. 
by Northern Analytical Laboratory, Amherst, NH. 


1. Metallic Nitrate Precursor Denitration/Decomposition 

The denitration of the individual cation precursors was studied by drying Ba(N(J3)2, 
CafNChb and CufNChb water solutions on alumina substrates. The thin film material 
was scratched off the surface and analyzed using differential thermogravimetric analysis, 
(DTGA). These measurements were performed by Robert J. Andrews with the Solid State 
Sciences Directorate, of the Rome Air Development Center, Hanscom AFB, MA using a 
Perkin-Elmer Delta Series TGA7 Thermal Analyzer. 


a. Thermal DenUration/Decomposition of Cu(NOj) 2 in Water Solution. 

X-ray diffraction spectra on sprayed films shows that Cii(N03)2 in an at lV e ° l *. s 
environment breaks down at elevated temperatures to form Cai2(OI1)3' no 3- DKjA 

analysis of the dried Cu(N0 3 ) 2 -water complex, sh " w " r' All Thermal 

temperature for maximal denitration/dehydration of the Cu(NC>3)2 is 265 C. All thermal 
analysis measurements were done in oxygen atmospheres. It appears that between 10 and 
210° C residual water is driven from the film residue. A very sharp transition is observed 
between 210° C and 310° C which is attributed to the denitration/dehydrm 
Cu 2 ( 0 H) 3 'N 03 - The 33.7% percentage weight change between 210 C and 10 1 U c agrees 



favorably with the formation of CuO from Cu2(0H)3'N03, on the basis of the following 
fractional weight percents: 

Cu 2(OH)3'NC>3 > CuO 4- oxides of nitrogen. 

fractional 

wt% 100% 66.2% 33.8% 

The chemical stability of CuO in oxygen at higher temperatures is evident. 

b. Thermal Denitration/ Decomposition of Ca(NOj)2 in Water Solution. 

Analysis on Ca(NC>3)2, (see Figure 2), shows dehydration below 200° C, and a 
percentage weight change of 66.1% between 400° C and 750° C. Unlike the cupprous 
reaction, a series of weaker outgassing peaks are observed in the first derivative spectra 
prior to a very pronounced rate of weight loss at 733° C. This data agrees quite favorably 
with the formation CaO through the following path: 

Ca(N03)2 > CaO + oxides of nitrogen 

fractional 

wt% 100% 34.16% 63.84%. 

CaO appears to be stable in oxygen above 740° C, after denitration. 

c. Thermal Denitration/Decomposition of Ba(NOj)2 in Water Solution. 

The barium nitrates, Ba(N03)2, are the most stable of the starting compounds. (See 
Figure 3). First derivative spectra indicates pyrolytic phenomena starts above 600° C. A 
stable regime is noted between 800° and 900° C. This stability breaks down above 900° C, 
but first derivative spectra shows the film tends towards stability again above 1010° C. The 
percentage weight change between room temperature and 1010° is 40.4%, which is close to 
the fractional weight which would be lost in the formation of BaO: 

Ba(NC>3)2 > BaO + oxides of nitrogen 

fractional 

wt% 100% 58.65% 41.35% 

The stable oxynitrate phase observed in the Ba(N03)2 derivative between 800° and 
900° C can not be identified at this time. The Ba(N03)2 DTGA data suggests that unless 
the presence of the other nitrates, (Cu(NC>3)2 and 'its derivatives, or Ca(N03)2), can 
accelerate conversion of Ba(NC>3)2 to BaO at lower temperatures, annealing above 1000 C 
will be required to completely denitrate the BaCaCuO precursor films. 

2. Effects of Solvents and Chemical Binders on Nitrate Denitration/Decomposition in Mixed 
Precursor Solutions. 

To study identify optimal solution chemistry and chemical binder activity the nitrate 
precursors were mixed in stoichiometric ratios of 212 Ba:Ca:Cu and dissolved in: 

A) A water-only solution 

B) An ethanol-(20% vol)/water solution 

C) A glycerol-(20% vo!)/water solution 

Each solution was dried on a hot (120° to 230° C depending on solvent volatility) alumina 
substrate. DTGA scans were run to observe effects of the solvents and the other nitrates on 
chemical reactivity during denitration. Ethanol solutions were considered to reduce 


fi to 



Ba(NCH)2 precipitation from solution. Glycerol solutions were considered to reduce single 
nitrate species agglomerations in solution and to enhance the chemical reactivity between 
the different nitrate compounds. The glycerol molecule provides a longer chain with 
hydrogen bonding sites which could conceivably complex different nitrate species in close 
physical proximity to one another. 

The fractional chemical weight distributions for solutions tested were determined 
using a 212 ratio: 

Ba(NC>3)2 Ca(NC>3)2 Cu(NC>3)2 

Ratio 2 1 2 

wt % 49.12% 15.49% 35.39%. 


a. The Effect of a Water-only Solvent. 

Figure 4 shows DTGA spectra for the nitrates when dissolved in water only. 
Thermogravimetric changes below 200° C are being attributed to water absorbed in the 
nitrate film by exposure to air. Copper denitration proceeds as before, unaffected by the 
presence of the other nitrate complexes. Its characteristic spectra in the first derivative is 
once again observed at 265° C. The percentage weight change observed between _28 C 
and 369° C was roughly 11%, and is attributed to copper denitration. This is in general 
agreement with an 11.9% weight differential that would be expected using the fractional 
weight percentage change observed during Cu2(OH)3'N03 denitration/dehydration 
between 220° C and 300° C. 

In mixed-water solutions Ca(NC>3)2 and Ba(NC>3)2 denitration appears to be shifted 
to lower temperatures. Principal lirst derivative spectra for the denitration transitions are 
observed between 375° C and 700° C. This is in rough agreement with the Ca(NC>3)2 
spectra, but is contrary to the findings for Ba(N03)2- It appears as though Ca(N03)2 and 
Ba(NC>3)2 denitration is accelerated by each other’s presence or by the copper oxides in the 
films. It is also plausible that Ca(N0 3 ) 2 and Ba(N0 3 ) 2 complex together in aqueous 
solutions. Calculations of the fractional weight percent changes using the individual 
denitration data predicts a total percentage weight change of 57.56%, in rough agreement 
with the 54.35% weight change observed in the mixed films. It seems reasonable to assume 
that in water-only solvents most denitration should be complete above temperatures of 700° 
C. 

b. The Effect of a Water-20% vol Ethanol Solvent. 

Figure 5 shows DTGA spectra from the 20%-vol ethanol solution. Spectral 
signatures in the 0th and 1st derivative spectra are identical, indicating Cu2(0H)3'N03 
denitration again proceeds unaffected by the presence of either Ca(N03)2 or Ca(N03)2 in 
ethanol solvents. An analytical comparison was not possible. 

Figure 6 shows the Differential Scanning Calorimetry (DSC) measuiemtnts in 
oxygen on the dried residue of the precursors mixed in water and 20% vol ethanol solution. 
Strong endothermic peaks are noticeable at temperatures just below outgassing 
temperatures identified by DTGA. (See Figure 5). 1 his suggests that the individual nitrates 
are breaking apart seemingly unaffected by each other s presence. The lack of any 
exothermic peaks in the spectra shows that no latent heats of crystal fusion are being 
released. Although the use of ethanol may help dissolve Ba(N03>2 and prevent its 
precipitation from solution, it does not promote the formation of fused crystal 
microstructures. 



Furthermore, the addition of ethanol to the solution reduces the liquid surface 
tension to such a degree that boiling during the drying process caused liquid drops carrying 
raw materials powders to spew from the substrate. Ethanol experimentation was 
discontinued at this point because of the adverse effects splattering would have on 
stoichiometric control. 

c. The Effect of a Water-20% Glycerol Solvent. 

Experimentation with glycerol solutions was conducted to try to improve film 
uniformity. Since glycerol is a chain molecule with many hydrogen bonding sites the 
individual nitrate species have a better chance of being physically complexed closer to each 
other if they bind at different sites along the chain. It is intended that the close physical 
proximity of the individual nitrate species provided by the glycerol chain would improve 
uniformity, inhibit single species clustering, and enhance the formation of BaCaCuO 
crystallites. Stoichiometric formation of BaCaCuO crystallites is preferred to CuO 
crystallites segregated from CaO/BaO polymorphous phases. This will help ensure chemical 
integrity during subsequent processing and single phase formation during thallium diffusion. 

Figure 7 shows that the qualitative character of glycerol-complexed BaCaCu- 
denitration in oxygen changes substantially. DTGA spectra below 182 C, (the boiling point 
of glycerol) is due primarily to dehydration. Spectra between 182° C and 250° C is attributed 
to glycerol outgassing. The most notable feature of this spectra is that Cu 2 (OH) 3 'N 03 , if it 
is at all present in the glycerol residue, no longer denitrates by itself. Furthermore, the 
largest relative change in weight occurs below 500° C. The weight change between 250° and 
500° C accounts for roughly 67% of all outgassing above 250° C. This indicates that 
CafN03)2 and Ba(NC>3)2 denitration participates in the copper oxide conversions and that 
BaCaCu complexes are being formed in solution. Roughly one third of the relative weight 
change occurs above 500° which can be solely attributed to Ca(N03)2 and Ba(N03)2 
(complexed) denitration. This data also shows that the microcrystallites formed with a 
glycerol binder are relatively stable in oxygen at temperatures up to —900° C. Changes in 
the first derivative spectra above 900° C suggest that the compound could break down at 
elevated temperature. 

Figure 8 shows DSC spectra under oxygen atmospheres for the mixed precursors 
formed in water/glycerol solvents. Three exothermic peaks are evident at 204° C, —300° C, 
and 341° C. This data, and the marked change in the DTGA spectra for glycerol solutions is 
direct evidence that the glycerol molecule chemically binds the individual metallic nitrate 
precursors, promotes reactivity between the precursors and the formation of microstructural 
crystallites in the sprayed films after thermal processing. An unusual feature of this spectra 
is that structural relaxation is not observed at temperatures slightly below the principal 
outgassing band between 250° and 500° C. 

3. Effect of Substrate Temperature on Film Deposition 

In chemical spray pyrolysis substrate temperature can affect both surface 
morphology and film stoichiometry. Thermal mass transfer, convection currents, vapor 
phase dynamics or evaporative effects and secondary chemical reactions can adversely 
affect chemical incorporation rates into the films. The thermodynamics of this deposition 
technique has not yet been fully characterized so process control has to be obtained through 
empirical study uniquely. To identify suitable substrate temperatures chemical 
incorporation was studied by spraying 223 BaCaCu-niti ate ratios dissolved in all t tree 
solvents at substrate temperatures between 180° and 350° C. 

Energy Dispersive X-ray Spectroscopy was performed on the films to evaluate 
atomic ratios. Figure 9 shows the atomic Ca/Cu and Ba/Cu ratios observed in the films 
sprayed using an ethanol 20%-vol solvent versus substrate temperature. The Ca/Cu and 
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Ba/Cu ratios used in solution are also shown. By simple inspection it is apparent that during 
the spray process both Ca and Cu are not tully incorporated since the Ba/Cu film ratio is 
always greater than the solution ratio and the Ca/Cu film ratio is always less than the 
solution ratio. 

Figure 10 shows copper and calcium incorporation rates against substrate 
temperature between 180° and 350° C. Incorporation rates were calculated assuming 
Ba(NCh)? is incorporated at 100%. Under this assumption a system of three linearly 
coupled equations can be constructed for which there are two unknowns. Percentages of Ca 
and Cu lost can be calculated against Ba incorporation using the ratio °f ^ . 

solution. The accuracy of the calculation can be tested by comparing the ratio of Ca/Cu in 
the films to value obtained through calculation. Figure 1 1 shows the percentage differences 
between the calculated and observed values. 

A very sharp increase in the copper incorporation and an enhancement in the 
calcium incorporation was observed at 260° C. This matches the peak denitration 
temperature for Cu NC>3)2. This would appear to suggest that film incorporation rates for 
the metallic cation > can be quite sensitive to oxidizing conditions and mass transfer 
phenomena at the s ibstrate surface. 

Since subst ate temperature has such a strong influence on the lesultant 
stoichiometry in th. hydroxynitrate films a hot plate was constructed to control substrate 
temperature to wi hin 1° C. This construction employs an Omega(TM) Series 21HH) 
Programmable Temperature Controller capable of ramp and soak temperature cycling. 

Films prepared from glycerol 20%-vol solvents were tar less susceptible to surface 
reactivity. Scanning Electron Microscopy (SEM) revealed that a 260° C substrate 
temperature provided the best surface morphologies. Plates 1-3 show surface 
characteristics of films prepared using water-only, ethanol (20% vol), and glycerol (20 7o vol) 
at that substrate temperature. Both the ethanol and water-only films show signs of jingle 
species agglomeration as indicated by the D1 GA spectra. Although the films prepared with 
glycerol had the best submicron microstructures, (See Plate 4), they remained quite porous. 
(Plate 5), with a high density of weak link interconnections between grain granules. 

Continued r search using glycerol (20% vol)/water solvents sprayed onto substiates 
heated to 260° C h is shown that solution stoichiometries can be reproduced in the films. 
Table 1 shows atomic percents measured in a series of films using semiquantitative Energy 
Dispersive X-Ray Spectroscopic analysis, by Photometries, Inc., Woburn, MA.. These 
results indicate that a strong correlation between solution and film stoichiometries can be 
achieved using these process parameters. Plate 6 is a backscattered scanning electron 
micrograph of sample 223-02. Good chemical uniformity and small microstructures can 
formed using these irocess conditions. 


B, Effect of Substrates on the quality of annealed films. 

Precursor fi ms were prepared on < 100> oriented epitaxially polished ytterium- 
stabilized zirconia (YSZ) and magnesium-oxide substrates. After spraying at a substrate 
temperature of 260° C in air the films were immediately denitrated by therma. bake out 
using a temperature cycle determined from thermal analysis measurements. 

In both instances microcracking was observed in the oxide precursors. 1 his is a 
shortcoming to any solution process where the films essentially have to be dried from a 
slurry and large volumes of reactant by-product is driven from their interior. This results in 
a generally porous texture. Films sprayed onto YSZ substrate had very poor adhesuin to the 
substrate after the denitration process, and were not suitable tor Tl-ditlusion. (See Plate S). 



Most of the YSZ-substrate films were damaged during the diffusion step either by film 
evaporation or flaking. 

Films deposited onto MgO substrates had sharply reduced microcrack densities and 
far superior adhesion. The MgO substrate films did indeed become superconducting after 
Tl-diffusion when suitable processing parameters were used. No increased microcracking 
was observed even after electrical characterization in liquid nitrogen. These results suggest 
very strong suitability of MgO substrates with Tl-Ba-Ca-Cu-O superconducting films 
prepared using solution processes over a range of temperature from 77 K to 1225 K. 

C. Thallium Diffusion Processing. 

The superconductive Tl-Ba-Ca-Cu-O phase is formed from the Ba-Ca-Cu-O 
precursor by diffusing thallium into the film. This can be done with a pellet of thallium 
oxide in a small porcelain crucible. The pellet is of sufficient size so that it does not totally 
evaporate during the short annealing step. The substrate containing the precursor films is 
suspended over the crucible and the assembly is placed in a quartz tube furnace at 850° C 
for 15 minutes. During the process, argon, air, oxygen, or a mixture of gasses are flowed 
over the sample. Only oxygen atmospheres were successful in preventing the sprayed 
precursor films from evaporating off of the substrate. The gas flow rate is adjusted 
empirically to optimize the oxygen content of the resulting film and to maintain sufficient 
oxygen overpressures above the films. 

D. Electrical Properties of tiie Superconducting Films. 

Figure 12 shows the R-T curve for the first superconducting sample obtained. The 
transition is sharp and the zero-resistance temperature is 112 K.. While this is not the 
highest temperature reached by the thallium containing superconductors, it exceeds the 
highest reported value for sprayed films by 15 K. [3] In all five superconducting films were 
produced with zero resistance temperatures ranging between 97 K and 1 12 K. 

Preliminary measurements of critical currents at 77 K in these spray pyrolyzed films 
range between 100 and 250 A/cm 2 . These values are low for high-field magnet applications 
and result in major part from the porosity of the films, excessive weak link interconnections, 
poor grain alignment, and nonuniform crystallographic orientation. Many of these 
deficiencies can b<. reduced by finer tuning the process. 

E. Crystallographic Orientation of the Ceramic. 

An important consideration if the construction of useful magnets is the orientation of 
the material in the films used. To obtain high critical currents needed for high fields, the 
crystal grains in polycrystalline films should be oriented with the c-axis perpendicular to the 
substrate so that the current Hows in the a-b plane. Thus, it was important for us to 
determine whether or not films could be deposited with a preferential orientation. After 
processing, we analyzed our films using XRD and determined that the films grew with the 
c-axis perpendicular to the substrate surface. 

Table II compares the x-ray diffraction patterns of Tl-Ba-Ca-Cu-O bulk powders to 
an epitaxial thin film prepared by rf-sputtering and the T c = 112K film prepared by 
chemical spray pyrolysis and thallium diffusion. The epitaxial film has total c-axis 
orientation in the plane of the substrate, and the powder diffraction has completely random 
orientations. The c-axis oriented patterns are highlighted and the normalized peak 
intensities are given to show preferential c-axis orientation in the plane of the substrate for 
the spray pyrolyzed film. 



F. Silver-doping Experimentation. 

Silver doping into the precursor films can be easily achieved using silver nitrate, 
AgfNQik which dissolves readily into water. Silver incorporation into the Y-Ba-Cu-O has 
been shown to greatly enhance critical current densities, ductility and other mechanical 
properties. For these reasons the feasibility of silver-doping by chemical spray pyrolysis was 
examined in Phase 1. 

Three silver-doped films (10%-, 20%-, and 30%-wt) were prepared and each showed 
negligible microcrack densities after denitration baking, regardless ot the substrate used, i.e., 
YSZ or MgO. Unfortunately in an attempt to phase separate the Ba-Ca-Cu-O precursor 
from the silver, by annealing at 1000 C, the silver-doped films evaporated. An undoped film 
was annealed at the same time and it did not completely evaporate which suggests that silver 
incorporation into the Ba-Ca-Cu-O precursor increases the vapor pressure of the 
compound. 

III. PLANS FOR FUTURE DEVELOPMENT. 

The successful demonstration of a deposition process tor high-T c superconducting 
films which can be applied to the preparation ot large-area thick tilms has now made it 
possible to consider the manufacture of liquid nitrogen-cooled superconducting wire tape. 
A cross-section of a prototype design for this construction is given in Figure I T As 
diagrammed, a stainless steel strength membered copper tape will be used as the base 
structure. Spray pyrolyzed MgO film will be deposited upon it to provide good adhesion 
between the superconductor and the support structure. Researchers at the Naval Research 
Laboratories, Washington, D.C., have successfully demonstrated that spray pyrolysis can be 
used to deposit 100 oriented MgO tilms on sapphire, fused silica, anti silicon. [4] 

Silver-doped TIBaCaCuO films will be sprayed through a mask to form tracks on the 
MgO-coated base ape. The use of silver is being considered to improve the ceramic’s 
ductility and mechanical workability, [5], as well as to enhance the material’s current 
carrying capability. [6] Applying the superconductor in the wire as tracks should also reduce 
A.C. losses and improve thermal and electrical shunting ot the stored energies and stresses 
which develop in the films and can quench a superconducting magnetic coil. A silver layer 
will overlay the superconductor to improve conduction to the copper shunt, and an 
insulating tape will cap the tape to inhibit high-voltage electrical discharging between 
windings when the tape is assembled into an operating coil. 1 his design is being considered 
because it provides sufficient structural support to the mechanical stresses expected from 
5-6 Tesla magnetic fields. 

A major objective in the successful demonstration ot a high-tield magnet will be the 
development of strong flux pinning potentials in the superconducting films. It will be 
necessary to pin magnetic fluxoids at defect sites with sufficient strength to resist the I .orentz 
forces induced by the driving currents. The tree movement ot tluxoids in the 
superconducting wire tape will dissipate energy and pinch-ott critical current densities 
needed to operate in strong magnetic fields. The value ot critical current densities in all 
type-II superconductors is strongly affected by processing parameters and is not necessarily 
an intrinsic property of the superconducting system To achieve this aim significant research 
will be dedicated to optimizing c-axis orientation of the films by adjusting deposition 
parameters, and using other salt reactants or chemical binders in fabricating the precursor 
films. The application of nuclear irradiation [ 7 j and mechanical shock compression jb] tan 
improve flux pinning potentials and will be evaluated in this development as well. 

Finally, once large area films have been deposited onto the tape structure which have 
suitable electrical properties, the spatial characteristics of currents flowing in the>-e tilms 
under applied magnetic fields will be measured. These measuiements will lx trutud to the 
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determination of the optimal topological surface structure for the high-field magnet. This 
work will be done with the MIT Francis Bitter National Magnet Laboratory in Cambridge, 
MA. Once these properties have been assessed, a coil utilizing the wire tape in a pancake, 
Bitter coil, or hybrid pancake-Bitter coil geometry will be designed to generate a >6 Tesla 
field and tested. 
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TABLE I 


DESIRED ATOMIC 
SAMPLE PERCENTAGES 


I.D. 

Ba 

Ca 

Cu 

223-02 

28.57 

28.57 

42.86 

233-02 

25.00 

37.50 

37.50 

243-01 

22.22 

'44.44 

33.33 

253-02 

20.00 

50.00 

30.00 

224-02 

25.00 

25.00 

50.00 

225-01 

22.22 

22.22 

55.56 


MEASURED ATOMIC SURFACE 

PERCENTAGES UNI- 

Ba Ca Cu FORMITY 

30%±2 28%±2 42%±2 GOOD 

28%±2 41%+.2 30%+ 2 POOR 

22%±1 46%+2 32%±2 GOOD 

20%+1 48%+2 32%+2 GOOD 

29% +.2 29% +.2 41% +.2 POOR 

26%±2 25%±.2 48% ±2 POOR 


Table I. Stoichiometric control achieved by spraying nitrate precursors 

dissolved in water and 20%-vol glycerol on substrates heated to 260 C . 


Table II-XRD orientation 


POWDER 

PYROLYSIS 

RF SPUTTERING 

fhkh I/lo 

fhkn I/lo 

fhkll I/lo 

002/5 

002/2 

002/9 


004/2 

004/5 

006/5 

006/4 

006/8 



008/4 


103/5 



105/1 


00JL0/9 

0010/8 

0010/ 9 

107/87 

107/3 



110/2 


200/100 

118/14 



119/5 



205/5 

0012/100 

0012/100 

210/17 

212/10 

214/8 

220/47 

1011/1 

200/2 

- 


1112/1 


300/18 

215/1 

0018/25 

0018/18 

316/14 




X-ray diffraction patterns for T^CaBajCin. The preferential orientations of a spray 
pyrolyzed film prepared during Phase 1 is compared to those ol an epitaxial rl- 
sputtered film, (Yolchikavva, et a!., Appl. Pliys., July 18, 1988), and to the random 
orientations of a powder sample, (Hermann, PHYS. RliV.) LETT., Vol. 60, No. 16, 
1988). 
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Figure 1. DTGA spectra of Cu2(OH)3-NC>3-the residue of aqueous Cu(N03)o dried 
at 170 C. 
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Figure 2. DTGA spectra of aqueous Ca(NC>3)2 dried at 170 C. 
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Figure 5. DTGA spectra of mixed 212-Ba-Ca-Cu hydroxynitrates dissolved in water 
and ethanol (20% vol) and dried at 120 C in air. 
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Figure 6. DSC measurement of 2-1-2-Ba-Ca-Cu hydroxyiiitrate residue mixed in 
water and 20%-vol ethanol solution in an oxygen atmosphere. 
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Figure 8. DSC measurement of 2-1-2-Ba-Ca-Cu nitrate residue mixed in water and 
20%-vol glycerol solution in an oxygen atmosphere. 
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Cu/Cu and Ua/Cu Ratios in Film and Solution versus Substrate Temperature 
1 SO r T — 1 1 1 



Figure 9. Ba/Cu and Ca/Cu atomic ratios in films sprayed from mixed BaCaCu 
nitrates dissolved in water and ethanol (20%vol) versus substrate 
temperature. Ba/Cu and Ca/Cu atomic ratios in solution are marked by the 
solid lines. 
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Cation Incorporation Kates vs. Substrate 
Temperature (Assuming 100% incorporation) 



Figure 10. Percent incorporation rales for Ca and Cu sprayed from mixed BaCaCu 
nitrates dissolved in water and ethanol (20% vol) versus substrate 
temperature. (Calculated assuming 100% Ba incorporation over substrate 
temperatures). 
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Assuming 1 00% Ra Incorporation unci Calculated Cu arid Cu 
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Fiuure 11 Percentage deviation between Ca/Cu ratios measured in films and I Cu/Cu 
8 ratios calculated from Ca and Cu loss est.mates assummg 100 /o Ba l.lm 
incorporation. 
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Figure 12. R-T curve lor TMJa-Ca-Cu-O superconductor prepared by spray pyrolyzing a 
B:i-C 'a-Cu-O precursor film and thallium diffusion to form the superconductive phase. 


Insulating Adhesive Tape 



Figure 13. Cross-section of the prototype design of Tl-Ba-Ca-Cu-O wire tape for use 
in a superconducting magnetic coil. (Not drawn to scale). 





Plate 1 . SEM micrograph of a spray pyrolyzed Ba-Ca-Cu-0 precursor prepared 
using a water-only solution. (Magnification 1250 X). 



Plate 2. SEM micrograph of a spray pyrolyzed Ba-Ca-Cu-0 precursor prepared 
using a 20%-vol ethanol/water solution. (Magnification 1250 X). 
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Plate 3. SEM micrograph of a spray pyrolyzed Ba-Ca-Cu-O precursor prepared 
using a 20%-vol glycerol/water solution. (Magnification 1250 X). 



Plate 4. SEM micrograph of the submicron microstiucture of the grain granules in a 
spray pyrolyzed Ba-Ca-Cu-0 precursor film prepared using a 20%-vol 
glycerol/water solution. Weak links connecting granules limit high critical 
current densities. (Magnification 12. 9K X). 
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Plate 5 SEM micrograph of a crack in a Ba-Ca-Cu-0 precursor film lifting off from a 
vtterium-stabilized zirconia substrate reveals the porosity of the preliminary 
films. The high density of weak links in the films will have to be reduced to 
obtain reasonable critical current densities. (Magnification 3.01 K X). 



Plate 6 Backscattered electron SEM micrograph showing the baseline chemical 
uniformity of the Ba-Ca-Cu-0 precursor film prepared using a 20%-vol 
glycerol/water solution. (Magnification 321 X). 
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Bamam Eaa 


The purpose of this session is to explore the interest of the magnetic bearing community in 
applying superconductors to magnetic suspension systems and to discuss the issues related to these 
applications. As most of you probably know, superconductors have good potential, at least 
theoretically; the current density of superconductors is an order of magnitude or a couple of orders 
of magnitude higher than the current density of copper. Typical values are about 10 Amps per 
square millimeter versus about 100 Amps per square millimeter and flux densities are also fairly 
high. The flux density of copper as used in magnetic bearings is limited by iron to about 2.3 Tesla 
compared to the flux density of superconductors which could be demonstrated up to about 20 Tesla 
for short samples but typically about 5 Tesla normally. So this theoretical potential could be 
translated into some practical applications to improve performance. I think that this is the particular 
motivation for this session. We have a panel that consists of people who 1 suppose most of you 
know. On the right side we have Gerry Brown from NASA Ixwis who was in the 
superconductivity area about 30 years back dealing with superconducting magnets and Dr. 
Abdelsalam from Madison Magnetics who has been working for the past 10 years on MSBS and 
associated systems and on the left side we have Dr. Warren Kelliher from NASA Langley who is 
working in the area of high T c superconductor development and applications, and on the left most 
side we have David Eisenhaure, President of SatCon Corporation who has expertise in Low T c 
superconducting suspension applications. The way in which we thought we could organize this 
panel discussion was to give some time to individual panel members to give us their views about 
the applications and possible potentials of superconductors and then to have audience participation 
for interactive discussion to identify other potential applications. We'll start with Gerry Brown to 
talk about his views. 

Gerald Brown 

As far as applications go, I will note the ones that we are interested in at NASA Lewis 
where cryogenic temperatures are inherently available, such as turbo pumps for perhaps the Space 
Shuttle main engines or orbital transfer vehicles. Liquid hydrogen is the most attractive from a 
temperature viewpoint because its way below the critical temperatures of the superconductors we 
are talking about and also its perhaps the safer application. I'm not sure whether people are that 
ready to accept heavy current windings in a lox pump where every metal present can become fuel 
We're also interested in vibration isolation in space at Lewis but because of the associated 
cryogenic support that you would need with high temperature superconductors I don't really see 
them as being applicable in that area except perhaps if you have a very large coarse isolation stage. 
For instance if you were going to isolate the entire experiment payload bay; get a rough isolation at 
that level where you could afford perhaps the cryogenic support system and then provide finer 
isolation at the experiment level. One of the key problems of high T c superconductors (many have 
already been mentioned), current density, looks as if it will come along. I was talking to John 
Steckly and he mentioned that there had been reports of current densities that rival that of NbTi, 
Niobium Titanium. Structural integrity of the ceramics is often focused on but this is not a new 
problem. It has always been assumed that the inter-metallics like Niobium-Tin and that son of 
compound would fail if you ever allowed the superconductor to go into tension so it’s always a 
prestressed compressive situation that you want this type of superconductor to be in. So that is not 
an entirely new problem. The degree of the problem may be different. AC losses seem to me to be 
a major consideration and these have been beaten for the lower temperature superconductors to the 
extent that AC coils can be operated at power line frequencies. There are about three different 
kinds of AC losses and they all respond to making finer filaments and twisting these filaments with 
ever tighter pitches and transposing these filaments by twisted twist and so forth. Whether you can 
actually manufacture the ceramic materials in this form is certainly not clear, we don't even have a 
big bulk superconductor that's usable yet, but if we expect to carry any substantial component of 
AC current in the superconductors these items will have to be faced up to. On the other hand, it's 
been suggested by some of the people at the University of Wisconsin, Roger Boom and company. 



that perhaps you don’t need to carry the highest frequency components in the superconductor 
itself If you do have a normal sheath then that part of the Fourier analyzed current that is of high 
frequency can perhaps be carried in the sheath and what is of low enough frequency to soak in 
through the skin depths of that conductor could then be carried by the superconductor, so perhaps 
you don’t have to get the entire AC signal or AC current to soak clear into the core superconductor. 
Then finally just harking back to the days of the old superconductors, it could take a long time. 
From the time that people started showing really good current densities m short smples ot 
Niobium Titanium and Niobium Tin it was probably 10 years or more before we had decently 
stable magnets that you could depend upon. So it wont happen over night m my opinion, s ar 
as the thing that's going to work first, I was very pleased to see what Dantam had presented a 
couple of talks ago. Everybody's attracted to superconductors because they can produce so much 
field and since force goes as the square of the field, if you can get five times as much field, say 10 
Teslas as versus 2, it’s very attractive, but the structural problem is substantial when you take that 
big of an increase so I think that easiest first application is going to be just the replacement ot 
copper by superconductor as Dantam was talking about. In the first place the conductor is not 
subjected to very much stress, the relative permeability of the iron or iron cobalt core gives some 
indication of relative stress. In fact I guess it may even go as a square of that but the winding in an 
iron core system doesn’t see that much stress unless you try to bend it around a square comer or 
something. Furthermore, the AC loss situation would be easier to handle there. One advantage 
that Dantam didn't go into in his particular study is that you can push the core materials into the 
saturated range if you have plenty of ampere turns to burn and I think that is an attracUve thing that 
could yield another factor of 2 on the load capacity. One further thing - people that come to this 
discussion from superconductivity tend to think of big coils running at 10 Tesla and no core and 
people that come from magnetic bearings tend to think of just replacing the copper windings. One 
thing about using the cores is that you do not have much in the way of fringe field. Hie distance 
that fringe fields go is going to be more in the order of the gap than of the coil diameters. Many 
applications that we’re interested in such as space applications, may be applications where 
substantial fringe fields are unacceptable. Without a core it’s hard to build a compensating set ot 
windings where the field doesn't project a long distance. If you did go with the pure higher 
temperature windings without a core you can get pressures that rival the best that you can get in oil 
films I forget the numbers exactly, more than 2000 psi I'd say would be easy. At least the field 
strength indicates you could get that much. If you are willing to put coils on the rotors you can not 
only get attraction but also repulsion but again the fringe fields are worse. I have already 
mentioned the AC loss but one final thing before you actually go and design a high temperature 
superconductor winding for liquid hydrogen temperature: I think an analysis should be done to 
determine whether high purity aluminum could do better. Somebody this morning or yesterday 
mentioned that you can get resistance ratios between room temperature and hydrogen temperature 
of 400 or better with very high purity aluminum so you really ought to check where you d be bettei 

off. 


Dantam Rao 

Thanks, Gerry, I would now like to hand it over to Moustafa to talk about his experience. 
Moustafa Abdelsalam 

Mine will be easier than Jerry because I’ll talk about the Magnetic Suspension and Balance 
Systems for example. This is where you have to use superconducting magnets for the support 
system. If you try to use copper magnets the size would become too large and some ot the 
magnets would get away from the model so that you will not be able to use reasonable size 
magnets for the wind tunnel. So I think the first application of low temperature superconductors 
would be in large gap magnetic suspension and levitation and there I think that it would be a 
mistake for somebody to wait for high temperature superconductors to use for those applications, 
because as Gerry said, and I think that everybody agrees, that development just to make a 
superconducting solenoid may be 10 years along the road, not something that will happen in the 



next couple of years. The gain that you will gel there from high temperature superconductors 
versus low temperature superconductors is not that great. You still have to use a cryogen to cool 
the magnet and the structural supports are the same, so unless you can make the high temperature 
superconductor material itself cheaper than low temperature superconductors, the gain there is not 
necessarily that much. You gain some on the refrigeration power but that is most of what you get. 
You know the comment about the aluminum for hydrogen temperature - you can use high-grade 
aluminum. You could actually get down to 2000 resistance ratio for that current density and that, 
with liquid Hydrogen, might be in the range of somewhere between 2 kiloamps per centimeter 
squared and 4 kiloamps per centimeter squared depending on the size and where you can get the 
surface cooling. The only problem that you find there is that the high-grade aluminum is very soft 
so something has to be done to structurally support it. Another application where you have to, just 
by nature, use superconductivity is for energy storage where you can just store energy in 
superconducting coils and for that you can't use any other methods. You know you are forced to 
use, just by the nature of the problem itself, superconducting material. 

Dan tarn Rao 

Thanks, Moustafa. Now 1 will turn it over to Warren. 

Warren Kelliher 

Thank you very much. Can I see a showing of hands of the people who attended the 
Magnetic Suspension Workshop two years ago? Do we have any of the people here? Just a few 
of you. One of the things that was outstanding to me is comparing the papers and the presentations 
that we have today versus what we had there two years ago. Because two years ago the 
announcements of these superconductors had just come out, you had the hype saying that they will 
solve the world's problems etcera. Then we went into the dismay of finding that these materials 
were not that easy to work with and everything else. What we have now is sort of a hard-core 
group of people who are working with the engineering properties of these materials and putting 
them into a fabricated shape so they can be put into magnets, both small and large gap type of 
materials and the work is progressing quite well from what I can see because we've learned an 
awful lot in those two years about what causes the high J c s that are necessary for keeping these 
materials in the superconducting state. We've got up to, in some of the single crystal materials, 1()6 
amperes per square centimeter. It shows the possibilities, even in the bulk, of getting into 
thousands of amperes per square centimeter and that's all we really need for most wire type 
applications. One of the things that I disagree with a little bit on the MSBS is that we do have an 
NTT tunnel here and that handles a cryogen and therefore you automatically have the cryogen to 
run an MSBS system. However, going to a Helium-cooled system causes the cost to become very 
prohibitive (a million dollars a day or something like that). So high T c material will provide a very 
nice MSBS system, at least for the NTF tunnel here, that will be able to support the model and get 
rid of the interferences associated with the sting. One of the persons that attended two years ago 
was Dr. Jack Crowe. I was hoping that he would be able to show up this time but he couldn't. 

He is in Florida state right now and is now Head of the Magnetics Institute. I hope the problem is 
solved there between Boston and Florida but this is again its thought of focusing some of the 
activities in the superconductor materials towards the magnetics field and I see much more of that 
going to take place in the future. The high T c materials are going to find much more practical use 
as a magnetic material in the future. 

Dantam Rao 

Thanks Warren. I think we appreciate your comments on real world applications and how 
far we are today towards these real world applications. Now, I would like to turn it over to David. 



David Eisenhaure 

As advertised, I'm Dave Eisenhaure and I am currently employed by SatCon Technology 
Corporation up in Cambridge, right down the street from the magnet lab that's in question. In fact 
we can talk to the people during lunch hour and get their views on the whole problem. SatCon 
specializes in magnetic bearing and suspension systems; we've worked with both conventional 
magnetic bearings and superconducting bearings and suspensions. In fact the linear optical disc 
bearing in the back of the room is one of our products - just to get in a little advertisement while 
I'm up here and got the mike. What I thought I would do, since I didn’t bring any viewgraphs, is 
tell you a little bit about some of our experiences with superconducting magnetic bearings and what 
drove us down that path and what really created the need. I think our experience may be very 
typical of other applications. I think it was four and a half years ago when we received a contract 
to design an electromagnetic actuation system for laser radar mirrors. For people who aren t 
familiar with them, these are extremely agile mirrors; they weigh several hundred pounds; tney^ake 
two to three hundred kilowatts to drive them; and they require control frequencies as high as 4UU 
Hertz. The typical approach to that problem is hydraulic. That s the competition in this 
application. The problem is that if you lose a drop or two of that hydraulic fluid you min the 
optical surface in the mirror so there's a lot of need to go to a magnetic suspension and drive 
system, or there's a lot of belief that going that way will eliminate some of these problems. 1 he 
difficulty is when you begin to look into that problem you find out that the inertias associated with 
conventional soft iron magnetics suspending and driving these systems are so high that you can t 
meet the system performance specifications. You, basically can't drive the actuators themselves 
even without the laser radar mirror attached. We began looking at superconducting magnetic 
suspensions and drives and these superconducting suspensions and drives were basically 
superconducting fields coupled with either normal conducting or hypercon ducting armatures. 

What we found out by going down that path was that with normal conduct ng armatures we could 
basically beat the specific performance requirements of hydraulics by a little bit and with 
hyperconducting armatures we could beat them by a lot. Since then weve had several additional 
hyperconducting laser radar mirror development programs and we've extended that technology to 
the advanced concept CMG that Jim Downer reported on yesterday. I think one thing to keep n 
mind is the comment that Pat Wolke made during his talk describing the different kinds of actuators 
that are available and which can be built for magnetic suspension and torq.nng systems. I think 
you know that when you're developing superconducting systems you have the same options. 

Every one of those actuators has a superconducting analog and perhaps you have even more 
because it is not necessary to use iron in all of these applications because of the very high current 
densities. I have been to a number of these magnetic bearing and suspension conferences and now 
we're kind of talking about a niche within a niche here because we're talking about 
superconducting magnetic bearings within the general realm of magnetic bearings and sinpensioi s 
and if you look at magnetic bearings and suspensions in general it s kind of a niche techi ology. 

It’s emerging and the ilace that it can attract research and development funds are the areas where 
it's either an enabling technology (you can't do it any other way) or just has an overwhelming 
advantage to conventional technologies. It’s got to be a lot better than the way it s currently being 
done or nobody's goi ng to do it and for magnetic bearings you see a few of those places being 
found right now. It's happening in canned pumps, it's happening in certain specialized <inds or 
compressers. Some people believe it’s happening in attitude control wheels and m some parts of 
the world it's happening in machine tool spindles for extremely high speed very precise machining. 
I think by looking at low that's happening with conventional magnetic bearings and suspensions, 
you can look at the p oblem and say: where can we really use superconducting suspensions - 
where they have a un que place in the world and there are some areas. One area is the applicat ion 
that we found for hig 1 performance CMGs, laser radar mirrors, those kinds of very high specific 
performance require! rents. That same kind of requirement is what s driving the people that at e 
trying to build Gigav att size power systems on space or to put multi-thousand horsepower 
electrical drives in submarines. They want to make them small, compact, and light. On the other 
side of things, there are some people that want to build rotating machinery with microwatt kinds o 


power levels that they're going to send to Neptune or something and maybe the zero resistivity 
requirements of superconductivity will enable that. Maybe it’s hard to build that machinery with 
conventional technology. There's a couple of things that haven't been mentioned here (Dave 
Trumper mentioned it a little bit) that is related to very precise positioning systems. One thing you 
don't usually think about in servo systems, because usually you're being buffeted by everything, 
is the internal noise in these systems. If you're building a very quiet, very precise positioning 
system, what you may find is that the Johnson noise in the actuators is the predominant noise 
source in the system. Well, superconductors give you a mechanism to eliminate that. I think, 
along the same lines in very many precise stable systems, like perhaps electron microscope slides, 
where people are talking about positioning accuracies of an angstrom or better, introducing any 
heat into those systems is a big problem. You don’t have the stability in the materials if you're 
introducing milliwatts or higher levels of heat and superconductivity provides an enabling 
mechanism for that kind of a problem. I think one of the things that a group like this can do is to 
provide some guidance to designers like myself as far as what are the applications that are really 
needed and what are the requirements for those applications. What should we be designing and 
why and that's all I have to say. Thank you. 

Dantam Rao 

Thanks, Dave. Now, I have a few words to say about the experience that we have at MTI. 
We have, of course, had a couple of contracts dealing with superconducting applications from SDI 
and NASA and basically the object of the study was to identify potential applications in SDI 
systems where superconductors could form a good marriage with the system requirements. I is a 
rather difficult study because of the fact that, as you know, most of the SDI systems are politically 
oriented - they change from day to day and what you see today in the newspaper may not be found 
the next day. Within that kind of uncertainty we did in fact study the Phase I architecture of SDI 
systems and found that there is a certain scope where the cryogenic fluids are naturally available 
within the system. Basically, in the space based engines, that is the power engines where the 
cryogenic turbo machinery exists and the cryogenic turbo machinery already has liquid hydrogen 
available as a fuel. That is one potential application. The other applications which we feel are 
there, as was pointed out earlier, the space shuttle main engine hydrogen turbo pump and similar 
turbo pumps that are being planned. The liquid hydrogen is already available and there is a good 
marriage between the requirements of the superconductor with the applications needs. The 
application demands that the bearing should be very stiff and high load. The high load and the 
high stiffness requirements are normally met now by rolling element bearings and it is difficult to 
meet the same requirements by any other bearing other than a very high stiffness superconducting 
bearing. In addition to that I may digress a little at this stage to say that when I looked into the 
applications I saw the basic trends in the world and I saw a divergence in the trends between what 
Japan does and what the United States does. The Japanese, most of the time, focus mostly on the 
commercial multi-unit kind of applications where the research dollars are spread over a number of 
units so that the people will benefit. They probably start to identify an application right at a top 
level where people will use it, mass produce it, and then develop the technology from that state 
whereas in the United States I see that most of the applications are targeted at one-of-a-kind 
applications and the research dollars that the government spends are concentrated on that 
application and there is actually a problem there. Let me assure you that is basically a socio- 
economic issue which probably some other audience member could comment on The other 
applications which we would have in mind - 1 think most of you are aware of the magnetically 
levitated trains where superconductors are used by the Japanese to levitate the entire train. You can 
consider that as basically a 5 degree-of-ffeedom controlled magnetic bearing except for the fact that 
the particular device is driven by a linear motor. Instead, the rotary bearings are driven by some 
prime mover, except for that difference the superconductors are being used there and the advantage 
of the superconductors there is that they open the gaps quite substantially because of the higher 
ampere turns that are generated by the superconductors. They use an eddy current levitation mode 
there and that could be a basic technology which the United States could adopt if the Maglev Trains 
come around here. The third application I think, which probably some of the members of the 



audience are aware of, potentially multi-unit applications, are Cryo-coolers where current Ciyo- 
cooler bearings are right now either rolling element bearings or gas bearings. I believe that there is 
a potential for the high T c superconductors playing a role in the sense that they could levitate the 
rotating shaft permanently. Once it is levitated it could run without contact and that would reduce 
the wear and the power requirements. Right now the gas bearings they use in the Cryo- coolers 
are limited by the very short clearances - roughly they are under about a 12 micron clearance^ You 
could see that the wear and tear potential is very high there with the gas bearings and the high 1 c 
passive bearings could be probably be an application there. The other applications which I think 
some of the members of the audience are already dealing with are the MSBS suspension systems 
and micro-g isolators. I think there we probably need some participation from the members where 
we could extend the applications into potential space devices. With this I'll turn over the floor to 
the members of the audience and invite them to participate in the discussion. 


Pat Wolke -Honeywell 

Just another chance to editorialize little bit 1 guess but to reiterate what Dave said these 
applications have to be enabling technologies, or something where you get a significant 
performance increase and I would be wary about overselling the capabilities of some of these 
devices. Maybe in the realm of spacecraft applications we talked about ball bearings having a 
limited life and wear out mechanisms. One of the key wearout mechanisms is eventual loss ot 
lubrication and I see in superconductivity that you also have effectively a consumable > n tne ^ 
crvogen and if we talk about applications where they've got a lot of liquid hydrogen onboard 
anyway, all for this other stuff, we've got it available for us. Well typically, m control applications 
where you might be i sing magnetic bearings, those are operational for the entire hteot the ^ 
spacecraft and the other consumables are sitting ready for a particular application. They don t want 
to pop the canister for those. Once you pop the can you can t reseal it, it leaks. So that 
consumable that you have is part of your control system and must be included in calculations of the 
lifetime of your device, so we have to be careful about overselling these things, getting people too 
excited and then being basically disappointed at the end. We have to be realistic up front and 
include all of the things that are necessary in the system. 


Dan tarn Rao 

Thanks, Dave, it is your turn to answer the question. 


Dave Eisenhaure 

I get to answer that? Yes, 1 agree and I think we see two kinds of things; one thing is_ 
people are going to send up satellites that they really need and they are going to get data back from 
those; the other is the kind of thing we see at SDI and that is where some grandiose Battle Station 
is going to be built in the sky and that may never happen. 1 guess there are two levels. My 
thought was that superconductivity is a niche within a niche. 1 personally think that tor a lot ot 
spacecraft applications, conventional room temperature bearings without superconductors arc just 
fine. It is not at all that clear you need superconductors to do those things. If you go down the 
path a little bit and you say, "we have some applications where we really need superconductors, 
we have to build a torquer that weighs 5 pounds and produces 4-5000 foot/pounds of torque in a 
direct drive application," maybe what we need is a space-rated refrigerator to really use these 
things. That’s what the people that are trying to put superconducting electric drives into ships 
believe. They believe that one of the critical pieces of equipment they need in order to use these 
guvs is a refrigerator and whether that refrigerator makes 70 degree Kelvin or whether it makes 4, 
ft's a critical part and we don't have it. You can buy it for your lab, but you can t buy one to put in 
a ship and you certainly don’t have one that can fly. Now, maybe magnetic bearings, either 
superconducting or maybe not superconducting are a part of that system, and you know, Ima 
little bit prejudiced and I think that perhaps they are, but you know that s one of the things that has 



to be thought about. If this is going to happen, that's one of the pieces that has to be done and 
someone has to do it and someone has to pay for it. 

Dan tarn Rao 

Thanks, Dave. I think that I’m going to offer a small comment on that. I would tend to 
agree with Pat’s comment that there might be a likelihood of overselling, because of the fact that 
for over the past 50 years of activity in the superconductivity area we have a very few mature 
applications where the dynamic environment is involved. One particular application where the 
experience was gathered over 15-20 years that comes to my mind was the superconducting 
alternator program and the government spent a lot of money on that program and the only 
comments I heard from the funding agencies were that the multimegawatt alternators haven't 
produced a single watt. 

Gerald Brown 

Just one quickie, in some missions the cryogen may not be a consumable. If it's supposed 
to last long enough so that they've found a refrigerator is better than just a certain amount of 
storage, then you're just using power instead of the fluid. 

John Steklv - Intermagnetics General Corporation 

I'd just like to make a comment. When you mention the superconducting alternator, I think 
that this is certainly one area that has been explored. However, magnetic resonance imaging is an 
area where superconductors are being applied on a commercial basis. It's the largest single 
commercial application of superconductivity. For those of you who are older you probably know 
what it's all about, for those younger fellows that don't need to have MR1 scans this is a Computer 
Tomography (CT) type device and you use 1 meter bore superconducting magnets. Now these 
didn’t show up all of a sudden fully developed. You needed to have extremely high uniformity, 
measured in parts per million. You had to have magnets that didn't decay to better than one part in 
ten to the seventh per hour and you also needed to have a refrigeration system that was compatible 
with use in a hospital. Now, I think that this goes back to the system comment, you need all of 
this in order to make it work, plus you have to have it economical. You can't just have somebody 
that wants to use superconducting wire for MRI because it just won't work. You have to use the 
wire that’s available with the cryogenics system that exists now with a cryostat and with the costs 
that it takes to assemble all of this. I think again that this is the reason magnetic resonance imaging 
is a viable application today. 1 think that any other application needs all of these elements. The 
other area that's very successful, again it's not quite as difficult, is high energy physics. They 
make large accelerators; there's 1000 superconducting magnets that are operating at Fermilab that 
have been operating for the better part of a decade and there's the supercollider which is just in the 
process of being started in Texas that will have 12000 superconducting magnets. All very' 
successful applications of superconductivity, but again, it involves cryogenics, it involves 
magnetic design that is particular to that application and the correct economics, and I think that 
when you have that the application will work. 

Dantam Rao 

Would anybody like to add to what John said? 

Moustafa Abdelsalam 

I think I agree with John, that if you have the application, you don’t wait and I think that 
this was part of that recommendation that I will offer. For example, for MRI, if we waited until 
we had high temperature superconducting wires to make the magnets, we would still be waiting. 
And at the same time, we would have lost all the experience that we got along the way. Somebody 
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else might have built them. For space applications, let me add that since we have hydrogen on the 
space station, it would make sense to use high temperature superconductor energy storage. Mnce 
you are not going to pay for the refrigerant, you know for hydrogen you can actually store it in the 
system and you can use the magnetic energy storage system instead of batteries in space. 1 here are 
some studies that use toroids, in that case, you don't have any stray fields on the station and you 
can actually store energy there and it’s weight effective compared to other systems. Another area 
for the high temperature superconductor that I think can probably be done soon, is the leads tor 
low temperature superconductors. A major part of the losses for low temperature 
superconductors, especially for magnets, is the leads for those magnets, I m talking about tor 
example MSBS where you have to have the leads connected all the time for control and then it you 
can make those leads high temperature superconductor then you can decrease a lot of the 
refrigeration power that you need for that. 

Dantam Rao 

Thanks, Moustafa. Regarding John's comment, 1 think that I would fully agree with what 
John said, that the current experience of superconducting magnets is quite extensive, that is we 
have very mature and technologically proven applications; as he pointed out, there are quite a lot ot 
MRI systems which are being used all of the time. The superconducting supercollider is one ot the 
systems that is being projected and an MSBS system could probably be one ot the potential 
applications but what I was trying to point out when 1 selected the cryo- alternator or 
superconducting alternator for my comments was that in my mind when we are dealing with 
magnetic suspensions we are dealing with mechanical devices and when the mechanical devices are 
rotating or vibrating or whatever, they are distinctly difterent from steady-state devices like M 
and superconducting supercolliders. Its dynamic machines that we are dealing with and their 
dynamic environment. The distinct and unique opportunity that the superconducting community 
had to demonstrate the utility of the superconductors in a dynamic environi lent was the 
superconducting alternator and that was the only system I thought could be compatible with the 
magnetic suspension requirements which the superconductors may face in the future if we think o 
inserting the superconductors into mechanical systems. 


Robert Humphris - University of Virginia 

1 have a question for the panel concerning superconductivity in magnetic bearings for 
rotating machinery. In most of the magnetic bearings with rotating machinery today they are using 

switching amplifiers as the power amplifier, to save power again, and of course these switching 
amplifiers work at frequencies from 20 kilohertz up to 1 (X) kilohertz, and I can just imagine the AC. 
losses involved at these frequencies. Is this going to be a real problem? 

Dantam Rao 

Let me comment on that. I think that there are 3-4 ways in which we could see the 
superconductors inserted into magnetic bearings One approach was, which Gerry Brown 
probably feels comfortable with, and 1 also feel comfortable, was immediate insertion by replacing 
the copper by superconductors. That is an approach which could be viable. If you do that then 
you are going to face the problem which you mentioned but there are alternative approaches. That 
is, instead of trying to expose the superconductors to high frequencies you can design the flux 
circuit such that the superconductors will act as permanent magnets in persistent cunent mode and 
the flux circuit is away from the alternating fields in a decoupled mode. If you do that then you can 
avoid the problem. 


Bob Hum phris 


In that case, what would be the advantage of superconductivity over a permanent magnet 

bias? 

Dan tarn Rao 

Okay, there I think the answer is that permanent magnets have a maximum flux-density of 
1.3 Tesla whereas the magnetic materials saturate at around about 2.3 Tesla. There's a difference 
of 1 Tesla, so probably the superconductors could play a role by increasing the bias flux-density. 

Un-identified Speaker 

I'm intrigued by the idea of energy storage that was alluded to a little bit ago. Maybe this is 
slightly out of the field which you're discussing here, can somebody give me a thumbnail sketch of 
what you’re alluding to? 

M to ustafa Ab de i salam 

Okay, there was a study at the University of Wisconsin, on superconducting magnetic 
energy storage for space applications, where you can ship toroids. D-shaped toroids, to the space 
station and use them to store energy. What you do there is like a persistent mode electromagnet, 
where you have the current flowing all of the time in the electromagnet. You are storing the energy 
as a magnetic field and you can draw in and take part of that energy, you can charge and discharge 
the electromagnet. It works much like a battery except that you are storing the energy as a 
magnetic field not as a chemical reaction and actually there is a Defense Nui lear Agency (DNA) 
Phase I contract study on building 20 megawatt hour superconducting magnetic energy storage. 
That's for utility applications where you level the utility power. You charge it over night where 
you don’t have much load for the utility and discharge it during the day when you have a peak load 
and this way you can level the load on the utility. The first phase of this study was just submitted, 

I think, in August, beginning of August of this year and DNA is studying two designs and I think 
that sometime next year DNA will decide which design is going to be built for the next four years 
so that's Phase II of the study. 

Kirk Logsdon 

I had a question for the panel. Given the practical limitations of keeping flux contained in 
the circuit and not spewing all over the place in certain applications, does anyone on the panel feel 
that research on better soft magnetic materials should be pursued in light of the superconducting 
coils that we have potentially coming on the scene some day? 

Dantain Rao 

I don't know, maybe I'll try to answer that question. It depends upon the system which 
you have in mind. If you have in mind aerospace systems, at least the systems that I have come 
across - the SDI space systems, most of them are extremely sensitive to flux lines. Let me back 
up. I think that there are two types of systems in SDI. There are Sensor Systems which sense the 
oncoming missile where they use infrared detectors. Those systems are very sensitive to magnetic 
flux lines. There are Weapons Systems which throw the missiles on the oncoming missile and 
those Weapons Systems may not be that sensitive, like lasers for example. But the Sensor 
Systems are going to be quite sensitive to flux lines. One way of containing the flux lines is soft 
magnetic materials as you were suggesting. The other way could be using the superconductor in a 
diamagnetic mode, which could probably be attractive, I don’t know. Probably there is a need to 
see the system tradeoff between the weight of a diamagnetic shield versus the weight of the 



magnetic material shield. I think that is a systems study, every specific system has to be looked at 
to see which approach looks better. The approach which we thought could be better is to channel 
the flux, to solve that particular problem using iron, at least right in the beg nning. 

Don Rote - Argon ne Nationa l Laboratory 

We did a study a couple of years ago on applications of superconductors and high 
temperature superconductors to transportation, and one of the things that we looked at was the 
question of SMES, Superconducting Magnetic Energy Storage, and what we found was that in 
terms of the amount of energy that you could store per unit weight you simply cannot beat a batte . 
for storing energy. There's no way around it, the mechanical stresses are sufficiently great that > 
the time you include the components necessary to retain the magnet form, and it gets even worse in 
the case of a toroid, that you simply cannot beat the battery. But there are certain cases under 
which you can beat a battery and it therefore becomes the same question that was raised before and 
that is what are the special applications where superconductivity does provide some benefit It 
does not provide benefit in terms of energy per unit mass stored but it does provide benefits in 
terms of how fast you can charge the system, how fast you can discharge it, how many cyclesyou 
can go through. Chemical batteries are very bad, notonously bad, when you try to charge them 
very fast, discharge them very fast, or expose them to many deep cycles. If the superconducting 
coil is properly designed you can charge it fast, discharge it fast, and recycle it many, many times 
provided the mechanical stresses can be properly accounted for. So there may be applications in 
the space environment, as there are on the ground, where you can take advantage in that regard, 
but don't be mislead by the notion that you can beat a battery for energy storage per unit mass, 
because to our knowledge, at least at the present time, it can't be done. 

Dan tarn Rao 

Yes that is an interesting observation, and I think that Moustafa would be the right person 
to answer the questions raised. I'm aware that the sponsor to the SMES project had the object of 
demonstrating the trade-off between the energy storage device versus batteries by experiments. 
Probably Moustafa would be able to add something more to that. 

Moustafa Abdelsalam 

I tend to agree with you. I'm talking about power per unit mass. If you try to get a certain 
level of power from batteries, you will not be able to do it. You will have an energy storage 
magnet to do that. Its not energy per unit mass. I should have mentioned that. 

David Eisenhaure 

One of the things that hasn't been mentioned too much today and which we hear about in 
the newspapers is the need to ground-test large space structures and 1 hear stories about NASA and 
other agencies flying larger and larger space structures and the question of how these things can be 
ground^tested and how you simulate a space environment. There hasn t been much discussion of 
that here and it seems like we’re probably at the right NASA Center to discuss that. 1 was 
wondering if any of the NASA people could comment on what the needs are and whether some ot 
these large gap suspensions would be applicable to that task? 

Dantam Rao 


I think the comment is open to NASA personnel here. 



Nelson Groom - NASA Langley Research Center 


Well, in answer to that question, I think that there jas applications for magnetic suspension 
to ground tests. As you pointed out, ground testing is very important. There are problems testing 
anything on the ground, as you know, but with large space structures, it is even more so. The less 
influence that you have on a structure, the better off you are, so the answei is yes. 

John Murphv. Rockwell 

I'd like to comment on current space vehicle design and currently it s influenced very 
heavily by the propulsion system and currently we think that the HO propulsion system 
(Hydrogen-Oxygen) is the best for specific impulse. There are other servicing fluids that we could 
pick, but for other reasons, the current design is for that. That leads you to hydrogen and oxygen 
in a cryogenic form or a slush form with high density fuel cells that we currently know about 
through Air Force technology. It looks like superconducting technology as a complement to that is 
very applicable to future space vehicle design and I have alluded to this with several of the speakers 
but it is certainly a fertile area that we should explore for future space vehicle design, specifically 
for SEI or interplanetary application where maintenance becomes a big problem and in-flight 
maintenance becomes a reality which you must address. 

Dantam Rao 

I tend to agree with what we tentatively conclude in our SDI sponsored study - that is 
basically liquid hydrogen was one of the most abundant fluids in the SDI systems and maybe the 
applications could target liquid hydrogen temperatures. But instead of targeting an application like 
that before you start experimenting with liquid hydrogen fluid straight away, maybe you could start 
gaining experience by conducting demonstrations using liquid nitrogen at a higher temperature, 
gaining experience then go down to liquid hydrogen. The two-tier approach could be a viable 
approach to develop and demonstrate the technologies. 

Colin Britcher - Old Dominion University 

I have a comment. I guess I'm sort of a believer in technology demonstration exercises for 
their own sake. I think that is a throwback to the early operations of NASA where they built X- 
Aircraft because they were kind of nifty things to build. (They learned an awful lot on the way of 
course.) I guess the LGMSS experiment that was described this morning counts as a technology 
demonstration exercise and I got the impression from Pierre DeRochemont that there was some 
effort to build a high temperature superconducting magnet with the National Magnet Lab, which 
may also qualify. Are there other efforts underway or should there be other efforts underway? Is it 
the right time for people to try to build large magnets that can carry high AC currents, that kind of 
thing? 

Dantam Rao 

If I understand your question, I think the question was— are there applications other than the 
MSBS for the large magnets-am I right? 

Colin Britcher 

No. The question is: should people be trying to do technology demonstration exercises; 
build hardware almost for its own sake, just to show the practicality of using superconductors on 
large scales or for unsteady load applications or for large force applications, that kind of thing. 
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Dantam Rao 


rather 

the continued funding from the government. 

David Eisenhaure 
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should happen. 

John Murphy 

David, to comment that the next flight vehicle that is a good 1S 

the* v the NASP experimental vehicle, that is due to come on line in about \9 y 

spare application for tins new complement technology. 

John Steklev 

1 iust wanted to make a comment on the fact that low temperature superconductors have 
heen developed recently with very low filament sizes using matrices that are lesistive rather t a 

iilSISiilK 

Generally the losses are frequency dependent, only so much per cycle, and these tondi 
typicallyoperate reasonably well up to fields of about one Tesla (Ac operation). 

Dexter Joh nson. NASA Lewis 

I’ve only been on board with NASA Lewis for about 4 weeks now and I’m just t getli 
learning about magnetic nuX'suspension of large 

wmmmsm? 



where magnetic bearings or magnetic actuation could play a part. Some of the other things that 
have been used have been air tables and that type of thing, although those have caused problems. 
I'm sure that magnetic suspension could very well be used and it's a big area that magnetic 
suspension could be applied to. 

David Eisenhaure 

Where did you do your thesis? 

Dexter Johnson 

In the Structural Dynamics Branch, but now it’s Spacecraft Dynamics Branch here at 
Langley and I did it with the University of Buffalo. 

Dantam Rao 

Are there questions? I think we have stretched it a little bit beyond our scheduled time, 
because it's tea time. Anyway I hope you are very, very thirsty. In conclusion, I would like to 
thank all of the audience who participated in the discussion and hope we have increased the 
awareness of superconductors in the magnetic bearing community and that the next time that we 
meet, we will have the opportunity of listening to more presentations where superconductors have 
been used in bearings or suspensions or similar devices. Thanks very much. 
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Response of the system can be adapted by using various 
sensing feedback methods in the feedback control loop. 
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Active Vibration Isolation Mounts 

By 

H. Ming Chen 
Richard Dorman 
Donald Wilson 


Abstract 

Several approaches toward reduction of vibration from operating machinery to 
ground through active electromagnet mounts are discussed. The basic approach 
to active mount design is to take advantage of the force attenuation charac- 
teristics of soft mechanical spring mounts and supplement their limitations 
with an active electromagnetic system. Techniques discussed include 
vibration cancellation approaches, as weLl as a method of altering mount 
stiffness and damping properties at the disturbing frequencies. Analytical 
and experimental results are presented encompassing the magnitude of force 
reduction and the stability characteristics of each technique. 

Description of Figures 

Figure 1 

The advantages and capabilities of an active mount are summarized in th i tabu- 
lation o t F i gu re 1 . 

Figure 2 

The objective of the active mount as summarized in Figure 2 is to isolate 
machinery or platforms from ground. This is accomplished by sensing the force 
transmitted to ground to drive an electromagnet. 



Figure 3 


This figure illustrates a typical mount using an elastomer to support the 
weight and an electromagnet to provide dynamic forces. 


Figure 4 

Vibration control techniques to be reviewed include the use of filters to 
control the stiffness and damping properties of the support (fixed or tracking 
filters), or the use of an Inverse transfer function. 


Figure 5 

The force t ransmi ss i hi 1 i t y curve illustrates the reduction in transmitted 
forces that can be achieved by altering the support stiffness and damping 
properties. This figure also illustrates that the primary advantages of force 
attenuation occur .it the lower frequency range. 


Figure 6 

Figure 6 further illustrates the influence of stillness and damping modifica- 
tions to the t ransi issbility of forces through the mount. 


Figure 7 


Controls of the active mount depend upon adjusting the stillness and damping 
of the electromagnet as a function of the transmitted force as detected by a 
force gage or loan cell. 
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Figures 14, 15 


These figures refer to the technique used in developing an inverse transfer 
function. The test rig is excited with a sinusoidal signal to the electromag- 
net. The amplitude and phase are recorded at the force gage relative to the 
excitation. This information comprises the transfer function across the 
mount. The control technique is to develop a similar amplitude and phase 
function through a series of filter networks. A 180° phase inversion of this 
signal will drive in opposition to the force gage signal and cancel the trans- 
mi t ted force. 


Figure 16 

Figure 15 is a block diagram of the proposed filter network to create the 
inverse transfer function. 


Figures 17. 18 

These figures illustrate the constructed circuit to duplicate the character- 
istics of Figures 14 and 15. Tests are presently in progress to check out the 
performance of this approach. 
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Figures 8, 9, 10 


These figures illustrate a test rig used to demonstrate force attenuation 
techniques . In Figure 8, the mount supporting a vibration shaker is shown in 
the foreground. The control function is performed digitally in a desk-top 
[BM-PC. The control signal is amplified with a linear amplifier shown as the 
box alongside the rig. A close-up of the rig in Figure 9 shows the force gage 
used as the control pickup. In Figure 10, the electromagnet is visible with 
:he mount disassembled. 


Figure 1 1 


This figure shows the computer CRT screen and the software control commands. 


Figure 12 

In this figure, the real time vibration wave is shown above an FFT plot of the 
vibration as picked up on the force gage. In this figure, the force gage ir: 
inactive . 


Figure 13 

This is a repeat pLot of Figure 13 wi Lh the active mount activated. The output 
from the force sensor has been reduced from an amplitude of 0.3151 v. to 
0.0252 v., as noted from the value of Y at the bottom ot the pLots. 
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Mechanical Technology Incorporated 

ACTIVE MAGNETIC ISOLATION MOUNTS 

Objective: 

Isolate Platform from Both Platform Induced 
and Ground Induced Vibrations 

r Static and 
f Dynamic Loads 

Method: 

Sensor Driven, Actively 
Controlled Electromagnet 
to Cancel Dynamic Forces 

^ Dynamic Forces 

Potential Benefit 

• 15 to 1 Amplitude Reduction 

• Eliminate Multiple Rafting 

• Compatible with Current Mount Designs 

• New or Retrofit Capability 

• Applicable to a Variety of Platforms 



FIGURE 2 




Mechanical Technology incorporated 

ACTIVE MAGNETIC ISOLATION MOUNTS 


Approach: 

• Electromagnetic Support in Parallel with 
Existing Mount 

• Elastomer/Spring Mount Carries 
Gravity Load 

• Electromagnet Reacts Dynamic Load Only 


t 


Gravity- dynamic 



Result 

Minimize Size, Weight, and 
Complexity of Active Mount 


FIGURE 3 




VIBRATION CONTROL TECHNIQUES 
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Mechanical Technology Incorporated 

ACTIVE MAGNETIC ISOLATION MOUNTS 


Approach: 

Minimize Dynamic Force Transmission 
by Actively Altering Support Stiffness and Damping 



Frequency Ratio — 



FIGURE 5 
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Mechanical Technology Incorporated 

ACTIVE MAGNETIC ISOLATION MOUNTS 


Approach: 

Control Circuit Alters Stiffness and Damping 


Electromagnet 



Adjustable 

Stiffness 


Result 

Control Parameters to Achieve 
Reduced Elastomer Stiffness 


FIGURE 7 
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FIGURE 12 
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FIGURE 15 






ONO SOKKI CF-920 MINI FFT ANALYSIS SYSTEM 
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STRUCTURAL AMPLITUDE RESPONSE 




ONO SOKKI CF-920 



FIGURE 15 - STRUCTURAL PUASE RESPONSE 




FIGURE 16 H- S07HAN 

rioutui i-o 5/24/33 


PROPOSED COMPENSATION CIRCUIT 
CLOSED LOOP/FORCE SENSOR 
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FIGURE 17 - AMPLITUDE RESPONSE OF CONTROL CIRCUIT 
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FIGURE 18 - PHASE RESPONSE OF CONTROL CIRCUIT 




VIBRATION ISOLATION of SCIENCE EXPERIMENTS in SPACE 
DESIGN of a LABORATORY TEST SETUP 
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Department of Mechanical and Aerospace Engineering 
University of Virginia 
Thornton Hall, McCormick Road 
Charlottesville 
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DISTURBANCE LEVELS 


Quasi-Steady or 'DC' Accelerations 

Restive Gravity Freauency_lMi) 

IE -7 0 to IE- 3 

Source 

Aerodynamic Drag 

IE-8 

0 to IE-3 

Light Pressure 

IE- 7 

0 to IE- 3 

Gravity Gradient 


Periodic Accelerations 

Relative Gr< 

witv Frequency (Hz) 

Source 

2E-2 

9 

Thruster Fire 

(orbital) 

2E-3 

5 to 20 

Crew Motion 

2E-4 

17 

Ku Band Antenna 
1 


Non-Periodic Accelerations 

Re|ative Q 

ravitv Frequency (Hi) 

Source 

IE -4 

1 

Thruster Fire 

(Attitudinal) 

IE -4 

1 

Crew Push-Off 
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VIBRATION ISOLATION OBJECTIVES 

1 . The Problem 

- Vibration Causes 

• Thruster fire 

• Machinery operation 

• Crew motion 

- Desired Isolation Range 

• 10E-5 to 10E-6 g 

• 0 to 10 Hz frequency 

- Passive Isolation Capabilities 

• 10E-2 to 10E-3 g 

• Not good at low frequencies 

2. The S olutio n 

- Active Control Isolation 

• Electromagnetic Actuators 

• Low Frequency Accelerometers 

• Digital Control 

• Effects of Umbilicals 
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INITIAL EXPERIMENT : 1~D ISOLATION 

1. Long-stroke Shaker 

2. Umbilicals 

3. Electromagnetic Actuator 

• Lorentz Type 

• Long Action Magnetic Actuator (LAMA) 

4. Sensors 

5. "Experiment” Mass 

6. Magnetic Supports 

7. Base 
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University of Virginia / NASA Lewis 

LONG STROKE SHAKER 

6.25-in, 158-mm p-p stroke 



SPECIFICATIONS Moo* "3 


APPLICATIONS 
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DESCRIPTION 


’"e APS "4 md APS ’24 
DUAL-MODE ' Powe' Ai^cMrers are 
resigned bpec f'CdiJy to provide drive 
cicwe r for etectrodynarruc snakers r 5e 
utoui stage -ses proven high 'enaoiiity 
cower operational amplifiers arranged 
•iy produce a oaianced output 

’ -e irrot^iers H ave Matures which make 
jPfdue'v suited ‘or studying the 
'*nam»c characteristics pt structures 
’■-ey mav oe operated n e-mer a soilage 
r current amonfier mone seiectapie from 
're front pare: 2nis operating mode 
s ector switch vjcMdtes shaner Or ve 
cowe' pte fr Jpt!on r either t :urre«t 
• . ' t,iqe mode ‘or joseoaton ot 
-Si/nanne oecav n mojC - r “s 


— e .cmpieteiv seif-conratreo ..r-ts are 
packaged n rugged aluminum ercicswres 
suitaote for Pence or rac* mountmg 
Forced air coceng and massive "eat 
sinks for the output devees nsure :c r '‘ r 
uous operation w>th a shaner denser r g 
r ated force ntc oioc*ed resistive or 
reactive oads 

A cur'ent mor tor Signal ava ao e 
me 'ear panei permits TorMonrg y re 
nstantaneous Output jren» amplitude 
arc pnase E.ectromc protect or : r :^\r t 
//lit detect an utput shcr-": -ground o r 
:ve r, oad tone hens arc ■e'-uve tm 
J r ve S : grai 

Mode< *Z-i-£P features an e*tendec 
pc we* .ptior wn-c*' prow ces mgr-er pea* 
power with the same average pewe- jf 
‘"p oas.c Mode '24 


SPECIFICATIONS 

Average Output, "to shaker r eactive mo 
~ eak Output nto snaker 'eacuve oad 
'^rrent Output peak i random r o>se 
; jrrent Output, continuous 
r, “Guency Range 
"i dt :> gra Voitage 
•V *jt mpedance 
v 4i: ^ „ -e ter red to '^a* 'utput 
'■ f*»*n Mon tOf output 

I'L" jl Power 

Ae. r Pane 1 'on rectors 
1 wer Output 
''put 2«r'*?rt Mon-tor 
AC Power 
A/e-ght 
^ize HxyviO 


Model 114 

•25 7 a r ~s 

250 /' A • , 

r> j a peak 
a 0 4 'ms 
; 2000 
2 J ceak 
*00 K jrm 
+ - tB 
25 J "V A 

'22 J 50 50 10C ■■A/ 

220 240 V .prnnai 

/VK3-3’ 5 Oamon. 

3NC %pe 2 J a 
Ota J-r-r Peiectau e 
25 D ” J -q. 

5 22 < ’ ’ < 3 25 r'enes 
jj « 422 « 235 -m 
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DASHPOT SERIES 160 


MOUNT NO HOLE 
Jir « STS' (7 98 x 9 525mm) 
RECTANGULAR OR 
J7S' {9 525mm) ROUNO 
PANEL THCKN€SS 
0* 12S* (3 17mm} 

mounting NUT AND 
lOCKWASHER PROVIDED 

recommenced 

TIGHTENING TORQUE 

4 im Las 

3/r *32 THREAD 
(9 525 mm) DlA 


750 (19 05 mm) 

MAX DlA 

900 ' (22 86 mm) 
MAX DlA’ WITH 
OPTIONAL CASE 
960 ' <24 9 mm 
MAX DlA 



Bore: 

.627* (16.0mm) 

Damping Coafflclant . 

Regular damping: 2.5#/ 
in. /sec. 

Super damping: lO#/in./ 
SGC. 

Maximum Pull Force 

4# (1.8kg.) 

Maximum Friction Force : 
Less than igm. 


Operating Tamparatura 
Rang#: 

-75®C to + 150°C 

Approxlmata Platon 
Waight: 

3gm. 

Cylinder Waight: 1st inch 
8gm 

Each additional inch: 
3.6gm. 
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Long Action Magnetic Actuator (LAMA) Diagram 
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QA-700 

SUNDSTRAND DATA CONTROL’S 
Q-FLEX^ SERVO ACCELEROMETER 


Featui . 5 

Cost-Eflective High Accuracy 

Field-Adjustable Voltage 
Sensitivity and Range 

Consistently Repeatable 
Accuracy and Stability 

Self-Contained Sensor and 
Electronics in One Small 
Hermetic Package 

Better than 1 micro g Threshold 
and Resolution 

Dual Built-In Test Capability 

Wide Dynamic Range 

Internal Temperature Sensor 
Thermal Modeling 
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PERFORMANCE 

Output Range _ 

Bias 

Bias Thermal CoefTcent ^ __ 

Current Scale Factor 

Scaie Factor ^e^ai Cce^'oer 

L'neantv Error 

input Ax.s Misalignment 

Resolution- Threshold 

Freouercv Response 
O-iO ^z 
*r ; -2Q0 hz 
200- 30*; Hz 

Natural tr equencv 
Damping R.V10 

ELECTRICAL 

input Vcitage 
Quiescent Cur r ent ""a* 

■ soiation ;ase to ail pins 
Temperature Sensor Output* 


r30g 



’ 3m A g err 
200 com ' C 'O rr 


g ra< 
• V -30 


’ "-'a* 



1 VOC 

20mA per suoc v 
■^egonms at 50 - CC 
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MODEL: lrntzl DATE: SEP-14-90 



Axisvm metric Se ctio n of Lo rentz A ctu ator 
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LORENT Z ACTUATOR : DESIGN EQUATIONS 

1. Assume permanent magnet operating point 
for maximum energy product : (-HI, B 1). 

2. Compute magnet flux, fm = * Am. 

3. Compute circuit flux, fc = H 1 * Lm / R, 
where R is the circuit reluctance. 

4. Compare fm and fc. 

5. Adjust operating point until fm = c = f, 
when actual operating point has been found. 
(In case of saturation, f = saturation flux 

in saturated circuit segment.) 

6. Calculate air gap flux density, Bg - f / Ag. 

7. Compute force capability, F = i * / * Bg, 
where / is the actuator current and / is the 
total length of coil wire in the air gap. 

8. Change actuator geometry or 
circuit / magnet material until desired 
force level is achieved. 
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LQRENTZ ACTUATOR ; INITIAL DESIGN 


Total length 

s 

3.17 

in 

Magnet outer diameter 

r 

3.20 

in 

Magnet inner diameter 

s 

2.10 

in 

Magnet length 

s 

1.00 

in 

Shell outer diameter 

= 

2.80 

in 

Shell base thickness 

s 

0.30 

in 

Pole-piece thickness 

= 

0.40 

in 

Core diameter 

= 

1.10 

in 

Air gap 

= 

0.17 

in 

Shell-to-core gap 

s 

0.50 

in 

Gap ratio 

- 

1 : 2.94 

Coil wire diameter 

= 

26.67 

mils 

Number of layers 

s 

4 


Total number of turns 

s 

450 


Maximum coil current 

s 

2.52 

A 

Air gap flux density 

= 

0.50 

T 

Saturation flux density 

s 

1.20 

T 

Maximum force generated 

s 

2.33 

Ibf 

Actuator weight (excl. coil) 

s 

4.12 

Ibf 
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I ORENTZ ACTUATOR : CQMPACT__DESJGN 


Total length 


3.17 in 


Magnet outer diameter 
Magnet inner diameter 
Magnet length 


1.95 in 
1.40 in 
1.00 in 


Shell outer diameter 
Shell base thickness 
Pole-piece thickness 
Core diameter 


1.95 in 
0.30 in 
0.40 in 
1.00 in 


Air gap 

Shell-to-core gap 
Gap ratio 


0.17 in 
0.20 in 
1 : 1.18 


Coil wire diameter 
Number of layers 
Total number of turns 
Maximum coil current 
Air gap flux density 
Saturation flux density 


26.67 mils 
4 

450 
2.52 A 
0.45 T 
1.20 T 


Maximum force generated = 1.93 Ibf 

Actuator weight (excl. coil) = 2.18 Ibf 
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LORENTZ ACTUATOR: FINITE ELEMENT ANALYSIS 


Motivation 

- Different geometric configurations 

- Leakage across shell-to-core gap 

- Saturation in core 

- Minimum weight 

- Reasonable cost 

- Circuit materials with different 

saturation levels 

- Effect of current-carrying coil 

on air gap flux 

- Fringing 

Other Considerations 

- Mesh effects (coarse/fine) 

- Nonlinear analysis (B-H curves) 

- High gradient regions 
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Case II : Circuit Material -- HYMU 80; 
No Current in Coil 



MODEL: lrntz2 DATE: SEP-17-90 
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Case III : Circ uit Material -- HYMU 80: 

2.52 A in Coil, at 0 degree Phase 
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Case IV : Circuit Material -- HYMU 80j 

2.52 A in Coil, at 180 de grees Phase 
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CONCLUSIONS 


Design complete 

- Experiment 

- Lorentz Actuator 

Construction in progress 

- Concrete base in place 

- Shaker and its amplifier bought 

- Magnetic supports under construction 

- Data acquisition system being 

developed 

Experiment operational by late 1990 

- Background vibration measurements 

- Testing with "disturbances" 

generated by the shaker 
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I. ABSTRACT 

The work presented in this paper covers the recent developments in the area of non- 
contacting displacement sensors for a 500 Watt-hour magnetically suspended flywheel 
energy storage system. "Pancake" [permanent magnet, PM, and electromagnet, EM] 
magnetic bearings are utilized to suspend the flywheel. The work includes a detailed review 
of commercially available non-contacting displacement sensors and their suitability of 
operating with magnetic bearings. In addition, several non-contacting displacement sensors 
were designed and constructed for this magnetic bearing application. The results will show, 
currently available, commercial non-contacting displacement sensors will not function as 
desired and that an inductive sensor was developed to operate within this magnetic bearing. 


II. INTRODUCTION 

A prototype 500 Wh flywheel energy storage system is being built for the NASA 
Goddard Space Flight Center. This energy storage system is targeted for spacecraft 
applications, for it exhibits high specific energy densities and can be used for attitude 
control. A conceptual view of the 500 Wh energy storage system can be seen in Figure 1 . 
The energy storage system incorporates three key technologies, interference assembled 
multi-ring composite flywheel, high efficiency brushless motor and generator, and magnetic 
bearings. This paper focuses on displacement transducers, which are a vital element in the 
suspension of flywheels, via magnetic bearings. These transducers or sensors detect the 
displacement of the flywheel relative to the stator portion of the magnetic bearing. This 
displacement is referred to as the measurand. In practice, measurement systems seldom 
respond directly to the measurand. More often, for ease in measuring, it is desirable to 



convert from one physical quantity to another by means of a transducer, pie conversion in 
this case is from displacement to a voltage level which is proportional to the displacement of 
the flywheel. Next this voltage signal is fedback through the control system, compared to a 
reference voltage, and a control current is applied to the electromagnetic coils, in the 
magnetic bearings, to center the flywheel. 



In selecting a non-contacting displacement transducer for a magnetic bearing 
application, requirements for these sensors must be specified. The requirements include, the 
nominal linear range, the resolution, the sensitivity, the frequency response, the. size of the 
transducer, the wire bend radius of the sensor’s lead wire, and the power consumption. The 
main performance requirements of the displacement transducer for the 500 Wh energy 
storage system are listed in Table 1. For this 500 Wh energy storage system the 
displacement transducers were to be moved from sensing the outside periphery of the 
flywheel to sensing the inside surface of the flywheel. The initial location chosen for the 
sensor was on the inside of the magnetic bearing. 



Linear Range: 

20 mils 

Sensitivity: 

60 V/in. 

Resolution: 

20 finches 

Power Consumption: 

less than 40 mW 

Size: 

Must fit on inside 
of a magnetic bearing 
0.250" dia. 1.000" lenght 

Operating Environment: 

Must operate in strong 
magnetic fields 

Table 1. Sensor Specifications for 500 Wh Energy Storage System 


III. BACKGROUND ON COMMERCIALLY AVAILABLE TRANSDUCERS 

Some of the types of non-contacting displacement transducers presently accessible 
are inductive, capacitive, and optical transducers. Although there are other types of non- 
contacting displacement transducers, such as radiation, ultrasonic, and air gauging, only 
inductive, capacitive, and optical types were investigated. The general advantages and 
disadvantages [for use in the 500 Wh energy storage system] of the inductive, capacitive, and 
optical sensors are presented next. Inductive sensor have adequate frequency response, 
small size, and have a relatively large linear range. The main disadvantage of the inductive 
sensor was uncertainty of operating in a strong magnetic field. Also the inductive sensor 
produced a sensing area that was tne shape of a cone protruding from the sensor’s probe tip. 
Any conducting mate? ial that would intersect this conical sensing field would affect the 
output of the sensor. The capacitive sensor could be custom sized to tailor fit the 500 Wh 
energy storage system, due to several companies that build custom designed capacitive 
sensors. Problems of low frequency response, very high cost, and stray capacitance fields 
basically prohibited the use of capacitive sensors in the 500 Wh energy system. The optical 
sensors were quite promising with very high frequency responses, large linear ranges, very 
low cost, and very small probe size with the use of fiber optics. The only disadvantage of the 
optical sensor was the problem of non-uniform surface reflectivity of the target surface. 
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A survey of the commercial vendors, was conducted to locate a displacement sensor 
that could be utilized with the 500 Wh energy storage system. From the search, several 
pmential displacement transducers were identified and several sensors were obtained for 
testing with the 500 Wh energy storage system. 


IV TRANSDUCER TESTING 

The experimental testing of each displacement transducer was separated into two 
parts The first consisted of voltage output of the transducer versus target displacemen 
& and the second dealt with performing experimental tests usmg tfie d.sptoement 
transducer in an actively controlled magnetically suspended flywheel system. This test 
wsl^ utihSd onf mag y ne.i : bearing, see Figure 2. fn the Wowing paragraphs the results 
of several of the sensors tested are presented. 







LYWHEEL 



POSITION SENSORS' 


Figure 2. Actively Controlled Magnetically Suspended Flywheel System 

Two different model sensors (KD-2400 and the ^ 

Kaman Instrumentation Corporation of Colorado Springs, CO. The KD-2400 sens( > r was an 
inductive type sensor and was the sensor used in the previous magnetic bearing systems built 
at the University of Maryland. The voltage output curve for this sensor is shown in Figure 3 . 
TOfsensirSble toJneet all of the stick system’s sensor requirements except for the size 
requirement The sensor’s probe was too large to fit on the inside of a magnetic bearing. 
This sensor was tested for the reason that it was readily available and worked ^cesrfully 
nrevinus magnetic bearing systems. The other Kaman sensor, KD-2300-lbU, passed an 
requirements for the stacf system and proved to be very versatile. The voltage output versus 
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displacement curve is shown in Figure 4. The sensitivity was easily adjustable to 60 volts per 
inch. The linear range for the KD-2300-1SU was only approximately 40 mils for an 
aluminum target ana 35 mils for a stainless steel target. An advantage of this sensor was that 
it did not detect ferromagnetic materials. 

The last sensor tested was the SPOT optical sensor designed and built at the 
University of Maryland. This sensor’s voltage output versus displacement curve is shown in 
Figure 5. Optical sensors produce voltage output versus displacement curves that have two 
linear portions, these portions are named the front slope and the back slope. Since this 
sensor was built in-house it proved to be quite flexible. The sensitivity and linear range were 
easily adjustable to match the required specifications of the stack system. 
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Voltaic Output (rolti) 



Figure 5. Voltage Versus Displacement for SPOT Sensor 



Once all the displacement transducer’s sensitivities were calibrated to 60 volts per 
inch, further experiments on the sensors were performed to reveal if these sensors could 
function on the inside of a magnetic bearing. The first experiment that was performed on all 
displacement transducers was to substitute one of the displacement transducer being tested 
with an existing Kaman KD-2400 inductive sensor in an existing magnetic bearing. All of the 
displacement transducers passed this test except for the SPOT optical transducer. The 
SPOT displacement transducer failed the test because of the target’s non-uniform surface 
reflectivity. The target in this case was the outside surface of the magnetic bearing system’s 
flywheel which is constructed of aluminum. The problem of non-uniform surface reflectivity 
was documented in Figure 6. This figure displays the voltage output of the SPOT and 
Kaman KD-2400 transducers versus time for a flywheel speed of 100 rpm. The peak to peak 
displacement of the Kaman KD-2400 sensor was approximately 0.05 volts, which translates 
at 60 volts per inch to 0.8 mils. The peak to peak displacement of the SPOT sensor was 
approximately 0.14 volts or 2.3 mils. In an attempt to eliminate the non-uniform surface 
reflectivity problem, different colored target surfaces and polished target surfaces were 
utilized. All these types of surfaces did not correct the problem. For this reason the SPOT 
sensor was dropped from further consideration. 



All the inductive sensors passed the suspension test. To conduct further experiments 
it was necessary to design and fabricate a mechanical fixture to house the displacement 
transducers within the magnetic bearing. Once the sensors were placed in the magnetic 
bearing two experiments were conducted, a voltage output versus displacement experiment 
and a suspension test. The voltage output versus displacement test was conducted using both 
the Scientific Atlanta and Kaman KD-2300-1SU sensors. For this experiment, an aluminum 
ring had to be placed on the inside of the return ring to provide the sensors with an 
aluminum target to detect. Figure 7. shows the curves created by the Scientific Atlanta 
sensor. There are three curves in Figure 7., one curve for no current applied to the EM coils 
and two curves with different currents applied to the EM coils. These three curves show that 
an inductive sensor can still function within a static magnetic field. For both the Kaman and 
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Scientific Atlanta sensors the sensitivity of the sensor decreased °nee mside the magneU _ ^ 
bearing. This change in sensitivity was due to the large magnetic fields fro 

magnetic bearing The sensors were recalibrated to 60 volts per inch for the next test 0 
suspension with the inside sensor. The flywheel was rotated to a low speed of a few hundred 
revdS^er minute and the output sisals from the Scientific Atlanta arid Kaman -KD- 
24W sensors were compared. The Scientific Atlanta sensor was tested on the inside of the 
magnetk ^ bearing the voltage output of this sensor was compared to the voltage output 

of one of the Kaman KD-2400 sensors, which was used to suspend the flywheel. The output 
signal of the Scientific Atlanta sensor did not match the signal produced by the Kaman 
sensor The reason why these two graph do not correspond was because the ’ Sc^tif c . 
Atlanta sensor was affected by the dynamic magnetic fields produced when the ^magnetic 
bearing was actively suspending the flywheel. To confirm this hypothesis a second 
eSemwTrun rhis experiment produced a strong alternating electromagnetic field 
directly in front of the Scientific Atlanta sensor. The sensor was fuced relative to the target 
which was the iron pin at the center of the EM coil. The output of the 5™ mt0 

and found to vary although the sensor and target were fixed. This test concluded that _ 
inductive sensors were affected by varying magnetic fields and could not be th^N^Sada^ 
the magnetic bearing. This test was also conducted on the Kaman and Bendy Nevada 
sensors and similar results were obtained. 



Since inductive transducers rely upon magnetic effects they are particularly prone to 
interference from external magnetic fields generated by the magnetic bearings actuators. 
Subseauent experimental investigation of one commercial transducer has shown the 
probable cause to be due to the transducer sensor casing. This casing is manufactured from 
a stainless steel which is mildly, but sufficiently, ferromagnetic so that the transducer 
calibration is affected by changes in the saturation level of the steel when it is immersed in 

an external magnetic field. 
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V. UOMD, RMIT & FARE TRANSDUCER DEVELOPMENTS 


From the above experimental work, it was concluded that the displacement 
transducers tested could not be adapted to the inside of the magnetic bearing. Commercial 
inductive sensors will not work because of the alternating magnetic fields and the optical 
sensors have problems with non-uniform surface reflectivity. It is for this reason that effort 
has been devoted to developing suitable position transducers for 500 Wh flywheel energy 
storage system, which meet fairly stringent requirements of reliability, simplicity of concept, 
robustness, and ease of application. In this work attention has been largely concentrated on 
inductive transducers. 

Both the inner and outer bores of the flywheel grow significantly as it spins from zero 
to a maximum speed of 80,000 rpm. In the case of the inner bore the radial growth has been 
estimated to be 0.43 mm [17 mils] at a speed of 80,000 rpm. On the other hand the rotor 
translational motion is limited by the touchdown bearings to +_ 0.15 mm [6 mils] about its 
nominal center position, which it is observed is significantly less than the radial growth. 
Consequently any rotor position measuring system which senses the position, of either the 
flywheel outer rim or its inner bore, must be able to differentiate between displacements due 
to these two sources. To overcome this difficulty differential transducers need to be used. In 
the case of inductive elements positioned to sense on diametrically opposite sides of the 
flywheel inner bore. If these inductors are connected in an electrical bridge circuit the 
bridge balance will be unaffected by radial growth, and will only sense changes in inductance 
due to translational motion. 

Experience has shown that sensitivities in the range 2 to 40 V/mm are achievable 
with acceptable out put signal-to-noise ratios, but the sensitivities are usually limited in 
applications to the range 2 to 4 V/mm. Commercially available inductive transducers 
typically operate with carrier frequencies of 1 to 2 Mhz, and sense physical displacement 
signals having bandwidths up to 10 khz. 

Inductive transducers require a metallic target, and if they are to sense displacement 
on the flywheel outer rim then it must have a metallic surface attached by some means. The 
limited strength of metals combined with the very high surface speed of the flywheel makes 
this a difficult task. While the difficulty of attaching a metal surface to the outer rim may 
possibly be overcome, mounting the position transducers so as to sense displacement on the 
outer rim of the flywheel is not recommended for other reasons. Firstly, there are problems 
in maintaining concentricity between the flywheel outer rim and its inner bore, both during 
manufacture and especially when it is running at high speed. The latter situation arises due 
to the effects of large rotor growth with increase in rotor speed, and the in homogeneity of 
the composite structure causing an eccentricity to develop between the flywheel outer rim 
and the inner bore which changes as a function of speed. Secondly, the transducer sensors 
need to be rigidly and accurately fixed to the flywheel outer support structure which in turn 
must be accurately positioned with respect to the bearing stator; a complicated 
manufacturing problem. Thus for both manufacturing and operational reasons it is sensible 
to mount the position sensor internally so as to measure the displacement of the flywheel 
inner bore relative to the magnetic bearing stator. 


VI. RELATIONSHIP TO MAGNETIC BEARINGS 

The remainder of this discussion will center on how inductive sensors may be used in 
electrical bridge networks to solve the collocation problem by ensuring the effective 
displacement sensing planes are along the bearing actuator planes of symmetry. Four 
possible arrangements for the inductive sense coils of the transducers are shown in Figure 8. 
The most ideal arrangement is shown in Figure 8(a) where the coils are mounted between 
the pole faces, as this enables direct displacement measurement in the bearing planes 
XaXa’ and XbX B '. For practical reasons this arrangement is difficult to implement and so 
will not be considered further. 
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The remaining possible arrangements can measure the displacements in the bearing 
planes of symmetry provided it is assumed that the bearing actuators and the 
motor/generator are rigid bodies. For example, let us consider the alternatives shown in 
Figures 8(b) and (c) where the inductive sense coils are shown by the arrow bars in each 
case. From simple geometry the displacements x; and X2 can t> e calculated from 
measurements xi’ and xi and for the case shown in Figure 8(b) are given by 

M 1 
rTTEf 


0) 




A f 
"■ITS? 


+ .To 


and 

't + A i ] 
T~+ ~T£l 


( 2 ) 


h = X\ 


M ' 


+ x 2 ' 


As long as tf/l is small the error due to uncertainties in tf can be neglected and the 
computed values of x\ and xi can be used in place of their exact values in the respective 
bearing controllers. 



If the motor cannot be considered to be rigid then the arrangement shown in Figure 
8(d) needs to be used. Here the sense coils A h A 2 , coils A 3 , A* coils B h B 2 and coils %, B 4 
are separately connected in a series. Since these coils are symmetrically displaced about the 
bearing planes of symmetry the transducer outputs will be *i and * 2 , where 

J + :r 

(3) .7, - ^ 


and 

( 4 ) 



+ % 

2 


In Figure 9. is shown a simplified schematic of the experimental inductive bridge 
transducer, which can be used with any of the mechanical arrangements shown in Figure 8. 
The sensor inductors are connected in a Maxwell impedance bridge whose output is fed to a 
synchronous demodulator. The output of the demodulator passes through a low pass filter 
which filters the residual high frequency modulation products as well as any extraneous noise 
induced into the circuitry. The filter output is an analog signal whose magnitude is 
proportional to the displacement. 



From the remarks given above concerning the use of stainless steel in commercial 
inductive sensors care was taken to only use non-ferrous materials, such as aluminum in the 
construction of the experimental sensors. Aluminum was also used for the transducer target. 
Experiments showed that using this construction made the transducers insensitive to changes 
in external magnetic fields. 

The inductive sensors describe above will be utilized with the magnetic bearings in 
the prototype 500 Wh energy storage system. The specifications of this inductive sensor are 
presented in Table 2. A view of the location of these sensors within the 500 Wh energy 
storage system in shown in Figure 10. 
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Operating Principle: Differential eddy current induction 


Linear Range: 80 mils. Can be increased to 200 mils 

Sensitivity: 60V/in. Maximum exceeds 130 V/in. 

Resolution: 20 /^inches at 60 V /in. sensitivity 

Frequency Response: 5 kHz (3dB) 

Output Range: ± 10 V (offset can be 

adjusted to zero output) 

Mounting: Sensors remote from signal 

conditioning (0-4 ft) 

Target: Metallic - non-ferrous 

Adjacent Metal: Insensitive to external magnetic fields 

Coil Inductance Range: 50 fxH to 130 /iH 

Table 2. Inductive Sensor Specifications 
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Figure 10. Location of Inductive Sensors Within 500 Wh Energy Storage System 


VII. CONCLUSIONS AND RECOMMENDATIONS 

The research conducted, showed that commercially available displacement sensors 
could not work with the 500 Wh flywheel energy storage system. Optical sensors were too 
sensitive to the non-uniform surface reflectivity and the inductive sensors were affected by 
the dynamic magnetic fields produced by the magnetic bearings. These conclusions led to 
the development of an inductive sensor, which was built and tested with the 500 Wh 
magnetic bearings. The placement of this inductive sensor within the 500 Wh energv storage 
system was altered and several sensors were used in a differential arrangement. The sensors 
will soon be incorporated with the final 500 Wh energy storage system. 
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ABSTRACT 


The development of a five channel electro-optical model position 
measurement system is described. The system was developed for the 
13 Inch Magnetic Suspension and Balance System (MSBS) located at 
NASA Langley Research Center. The system consists of five linear 
photodiode arrays which are illuminated by three low power HeNe 
lasers, an assembly of lenses and mirrors for shaping the beams, 
and signal conditioning electronics. The system is mounted to a 
free standing, pneumatic isolation table to eliminate vibration 
susceptibility. Two distinct channels are used for sensing 
vertical displacement and pitch, providing an angle of attack ( AOA) 
measuring range of 40 degrees with a +.02 degree precision. 


INTRODUCTION 


The 13 Inch MSBS was originally constructed by Arnold Engineering 
and Development Center in 1965 and was relocated to NASA LaRC in 
1979 [1], Model position was sensed by a system of X-ray sources 
and detectors which produced an analog signal proportional to their 
exposure to the beams. Four beams were directed diagonally across 
the test section at the model to detect pitch and yaw, while a 
fifth vertical beam detected axial motion. A suspended model would 
partially block the X-ray beams, leaving a portion of the detector 
unexposed, or shadowed. Although the X-ray system was suited for 
detecting model position through the aluminum test section, the 
potential safety hazard prompted development of an optical system. 
The optical system operated on the same principle of blocking a 
portion of a laser beam (light sheet) with the wind tunnel model 
and sensing position using light sensitive detectors. A new 
subsonic wind tunnel with a transparent test section was designed 
to accommodate the optical system. 

The original optical system design was dictated by existing 
hardware used for mounting of the X-ray system. A series of poles 



surrounding the test section allowed fore and aft sliding of the 
detectors for different lengths of models. Many custom adapters 
were required to mount the optical system components to the 
existing structure. The optical components themselves, were 
limited to very short focal lengths in order to fit in confined 
areas . 

Installation of this system was completed in 1986 and it was 
successfully used for several years with minimal modifications. 
The system did suffer some deficiencies, however . Building 
vibrations were parasitic because of the critical optical align- 
ment. The position sensors were mounted to precision translators 
partially constructed of magnetic materials. Magnetic field 
intensity variations would sometimes vibrate the position sensors 
into a resonance. Optical realignment for various models was 
tedious and time consuming, often consuming an entire day. The 
limited AOA range of ±8 degrees did not meet the requirements of 
the 13 Inch MSBS, initiating design improvements. 


SYSTEM DESCRIPTION 


The original system consisted of five identical channels using 
the same 1-inch position sensors. The new system uses longer 
position sensors for two of the channels to improve the AOA range, 
resulting in a combination of two separate subsystems. In 
addition, a new mounting system for the optics and sensors has been 
implemented. The following is a detailed description of the new 
system components. 

Position Sensors 

Linear, self-scanning photodiode arrays (PDA) were chosen to 
detect model position. A pair of 1-inch arrays are used for- 
sensing position in side and yaw while a third senses axial 
position. Each of the 1-inch PDAs contains 1024 individual 
elements on 25 micron centers. Each device is mounted to a 3 inch 
square factory supplied circuit board. The circuit boards contain 
all necessary clock circuits for operation. Modifications allowed 
the circuits to be operated in a master-slave configuration, 
driving the arrays with the same 500 KHz clock. 

A pair of 2.4-inch arrays are used to detect vertical and pitch 
position. Each of these arrays contains 4096 elements on 15 micron 
centers. The commercially available units are presently driven by 
modified factory circuit boards. A master-slave configuration with 
a 1.1 MHz clock is also used for these units. The 13 Inch MSBS 
system computer provides a 256 Hz scan initiate signal, which is 
synchronized with both of the PDA sampling clocks. 


Light Sheet Generating Optics 

Two different light sheet generating methods are used for PDA 


illumination to accommodate different requirements for each 
subsystem. The method used to illuminate the 1-inch PDAs is less 
critical since these are normally operated in saturation. While 
the device is saturated, it is not sensitive to small intensity 
variations along its aperture. Minor scratches on the test section 
windows and dust on the optics generally do not degrade the system 
operation. 

A 5mw HeNe laser beam is fanned by using a 4mm diameter glass rod 
as a high power cylindrical lens. This quickly expands the beam 
into a wide angle, reducing the optical path length and saving 
space. The sheet is then collimated by a 50mm diameter plano-con- 
vex cylindrical lens with a focal length of 150mm. Since this 
combination loses very little power, two PDAs can be illuminated by 
use of a plate beamsplitter. Figure 1 illustrates how the two side 
and yaw detectors are illuminated in this manner. The axial array 
is on its own axis and is, therefore, illuminated by its own 2mw 
laser and associated optics for simplicity. 

The two 2.4-inch arrays used for detecting lift and pitch are not 
operated in saturation and are very sensitive to intensity 
variations along the aperture. Several different lens combinations 
for beam expansion, including commercially available units, were 
evaluated before satisfactory results were obtained. A double 
glass rod expander, which spreads a section of the first light 
sheet a second time to overcome the Gaussian profile was successful 
but very inefficient since the resulting light sheet is many times 
wider than the diameter of the collimating lens. A plano-convex 
cylindrical collimating lens was believed to be necessary to avoid 
increasing the thickness of the light sheet. A lens was custom 
ground to meet size and focal length requirements without realizing 
the effects of distortion on the system accuracy. A laboratory 
calibration revealed that the collimating lens caused large errors 
due to spherical aberration and was replaced by a spherical 
achromatic doublet of longer focal length. Construction of a 
second light sheet by inserting a plate beamsplitter caused 
diffraction patterns, ghost images, and added to the intensity 
variations . 

Figure 2 illustrates the present light sheet generation for the 
two vertical sensors. The expanding optics are basically Keplerian 
design so that a spatial filter may be used. A 5mw multimode laser 
was found to improve the uniformity of the light sheet due to its 
flatter profile. A single plano-convex cylindrical lens with a 6.4 
mm focal length spreads the beam while a plano-convex spherical 
lens with a focal length of 600mm and diameter of 95mm performs the 
collimation. The inherently small spherical aberration of the high 
f-number lens reduces system errors at the expense of a long 
optical path length. Errors may be further reduced with a 
substantial increase in component cost. Unlike a cylindrical lens, 
the spherical collimating lens causes the light sheet to change 
thickness over the path length. If the PDA is located near the 
focus of the lens, the resulting line will be extremely narrow, 
making alignment more critical. The greater the distance between 



the detector and the lens focus, the thicker the light sheet. A 
small amount of spreading is desirable to make the system less 
susceptible to vibration. 

Several advantages are gained by splitting the laser beam prior 
to expansion using a cube beamsplitter. The cube beamsplitter does 
not rely on uniform metallic deposition and the splitting process 
has no direct effect on the quality of the light sheets. Although 
a duplicate set of components is required, the height of the two 
optical channels may be independently adjusted to the required 
optical axis. By rotating the beamsplitter about the laser beam 
axis, the reflected beam angle is varied without affecting the 
original beam. This feature allows a much higher AOA when equal 
positive and negative angles are not required during the same test. 

Signal Conditioning Electronics 

Two signal conditioning schemes are used for each of the PDA 
types to accommodate two different video signal formats. The 

inch arrays are operated in saturation and produce a serial train 
of pulses proportional in amplitude to the exposure of each of the 
1024 photodiodes sampled during a scan. The 2.4— inch arrays 
produce a sample and hold analog signal. The longer array is also 
more likely to encounter two shadow edges, resulting in two level 
transitions on the video signal. Detecting two shadow edges and 
their location requires more sophisticated electronics. The 
following will describe the two signal processing schemes. 

The pulse video signal is carried by coaxial cable to the signal 
conditioner where a voltage comparator with an adjustable threshold 
level shapes the signal into TTL compatible rectangular pulses. 
The pulses are fed to binary ripple counters which count the total 
number of pixels exceeding the threshold level. The total number 
of illuminated pixels is then retained by latches for the duration 
of the following scan. The actual location of the shadow edge is 
easily determined by the system computer since one end of the PDA 
must always remain illuminated. 

As with the pulse video processing, a voltage comparator is used 
at the input cf the sample and hold video signal conditioner. The 
comparator tr. nsforms the video signal into a clean TTL level pulse 
by triggering on the voltage level changes occurring at the shadow 
edges. The comparator's hysteresis band along with various timing 
delays of other components help reject the diffraction effects on 
the video signal level transitions. Polarity sensitive monostable 
multivibrators are used to detect light— dark and dark-light 
transitions . 

The location of a shadow edge is determined by counting PDA 
sampling clock pulses during each scan. The count is continuously 
fed to two separate registers during the scan. At the occurrence 
of a video level change, the appropriate register is enabled, 
temporarily storing the data until the end of the scan. The data 
is then shifted to the output registers, where it remains for the 


duration of the following scan. An output enable pulse from the 
system computer controls which tri-state register \ ill be accessed, 
requiring only one set of lines for both outputs. A balanced, 
optically isolated, line driver/receiver system is used for signal 
transfer between the computer and signal conditioner. The 256 Hz 
scan initiate pulse from the system computer is used for various 
resetting and clock synchronization. Shown in Figure 3 is a block 
diagram for one of the two channels. 

Mechanical Assembly 

An aluminum framework was designed for mounting the position 
sensors and a portion of the optical components. The mounting 
system basically consists of two arches joined by two sets of poles 
allowing fore and aft translation of the thn e 1-inch PDAs. 
Precision linear translators allow fine lateral positioning of the 
1-inch arrays. The arches were designed to fit between the test 
section and the electro-magnets. Optical rails mount to either 
side of the assembly, supporting the entire optical system for the 
vertical sensing channels. The 2.4-inch PDAs are mounted to custom 
built aluminum translators for vertical adjustment.. The framework 
assembly is attached to a commercially available 36" X 24" optical 
breadboard on which the remaining optical components are mounted. 
The aluminum breadboard provides flexibility of optical layout as 
well as accepting a broad range of optical holders. Ambient 
building vibration and magnetic field induced vibrations are 
minimized by floating the entire system on a pneumatic isolation 
table. The table, which suppresses vibration above 2 Hz at a rate 
of 12 db/octave, physically isolates the position sensors from the 
electro-magnet structure. Figures 4 and 5 are photographs of the 
optical system before and after installation. 


CONCLUSIONS 


A new model position detection system has been installed and 
demonstrated at the 13-Inch MSBS. System AOA range and general 
reliability is superior to that of the original system while 
meeting cost and development time requirements. Although the 
system has not been calibrated in the test environment, separate 
laboratory calibrations of the 2.4-inch and 1-inch systems were 
performed [2], [3]. Overall precisions of +.002 and +.0005 were 
displayed in linear measurements while angular precisions of ±.02 
degrees and +.015 degrees were exhibited by the 2.4-inch and 1-inch 
systems respectively. Mechanical difficulties remain to be a 
primary shortcoming. The physical size of the factory circuit 
boards limit the vertical adjustment range of the 2.4-inch PDAs, 
directly affecting the AOA range. Smaller clock-g ;nerat ing circuit 
boards are currently under development to alleviate this problem. 
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Summary 


A phologrammetric optical position measurement system is currently being built as part of the NASA Langley 
Research Center Large-Gap Magnetic Suspension System (LGMSS). The LGMSS is a five degrec-of-freedom, 
large-gap magnetic suspension system to be built in the Advanced Controls Test Facility (ACTF). The 
LGMSS consists of a planar array of electromagnets which levitate and position a model containing a permanent 
magnet core. Information on model position and attitude is required to control the position of the model and 
stabilize levitation. The optical position measurement system determines the position and attitude of a levitated 
model in six degrees of freedom and provides this information to the system controller. Eight optical sensing 
units positioned above the levitated model detect light emitted by small infrared Light Emitting Diodes (LEDs) 
embedded in the surface of the model. Each LED target is imaged by a cylindrical lens on a linear Charge 
Coupled Device (CCD) sensor. The position and orientation of the model are determined from the positions of 
the projected target images. A description of the position measurment system, tracking algorithm, and 
calibration techniques, as well as simulation and preliminary test results of the position measurement system 
will be presented. 


Introduction 

A Large-Gap Magnetic Suspension System (LGMSS) is currently being built at NASA Langley Research 
Center to test control laws for magnetic levitation for vibration isolation and pointing. A photogrammetric 
optical position measurement system has been designed and is currently being fabricated at NASA Langley as a 
part of the LGMSS. The optical position measurement system will measure the position and attitude of a 
levitated body in six degrees of freedom and supply this information to the LGMSS control system. This paper 
will describe the requirements for the position measurement system, and present the design of the system as well 
as some preliminary test results on a prototype sensor. A high level description of the system will first be 
presented. This will be followed by a discussion of the tracking and calibration techniques. Lastly, details of 
the system design will be presented and some conclusions about the performance of the optical sensing system. 
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Presentation Outline 


Position Measurement System Requirements 
Overview of the Sensing System 
Tracking and Calibration Techniques 
Closeup of the Sensing System 
Experimental Results 



Position Sensing Requirements 


The requirements for the position measurement system arc derived from the requirements for the LGMSS. There 
arc three fundamental requirements for the position measurement system. These requirements are: 1) that the 
sensing system be able to track the levitated body to an accuracy consistent with the accuracy requirements for 
the LGMSS; 2) that the position information supplied by the sensing system be supplied at a high enough rate 
that the LGMSS can stabilize levitation ol the model; and 3) that the measurement made by the sensing system 
not interfere with the function of the LGMSS. 

To meet the three requirements specified above, a photogram metric optical sensing approach was selected. 

Small point targets embedded in the surlacc of the levitated body arc delected by multiple cameras positioned 
about the model. The locations of the projected target images arc determined and triangulalion techniques are 
used to determine the position and altitude of the model from this information. The accuracy requirement for the 
optical sensing system is set by the accuracy requirement for the LGMSS. The LGMSS is required to position 
a levitated model to an accuracy of +/- 0.01 inches in x, y, and z, and +/- 0.02 degrees in yaw, pitch, and roll 
(defined by the Euler angles psi, theta, and phi). The optical sensing system is allowed thiry percent of the error 
budget and thus the required accuracy for the sensing system is +/- 0.003 inches in x, y, and z, and +/- 0.006 
degrees in psi, theta, and phi. The frequency response for the optical sensing system is 20 samples/second. 

This is not a fast update rate, but is fast enough that the LGMSS can control the levitation of the model to the 
accuracy specified for the system. 



Position Sensing Requirement 


Position Measurement Accuracy 

+/- 0.003 inches xcm, ycm, zcm 
+/- 0.006 degrees in psi, theta, phi 

Frequency Response 
20 samples/second 


Noninvasive Measurement 


O verview of the Optical Sensing System 


In figure J is a block diagram of the optical sensing sytem. Eight infrared light emitting diode (LED) targets are 
embedded in the surface of the levitated model. The targets are multiplexed in lime for target identification. As 
each target is flashed on, it is detected by sixteen linear charge coupled device (CCD) array sensors. A 
cylindrical lens, positioned in front of the detector focuses the light emitted by a target as a line of light on the 
detector array. Each array has 2048 photosensitive elements. A voltage signal is output from the array which is 
proportional to the light falling on each element or pixel, (figure 2) The location of each target image along the 
CCD array is determined by calculating the location of the centroid of the light distribution falling on the 
detector. The analogue video signal output of each CCD sensor is digitized and stored in a random access 
memory (RAM) which can be accessed by a high speed digital signal processor (DSP). There are a total of 
sixteen analogue to digital converters (one for each sensor), and sixteen DSPs (Texas Instruments 
TMS320C30s). The location of the centroid of light is calculated by the DSP and stored in a central memory. 
The locations of the target images in the sensors are used to determine the position and attitude of the levitated 
model. Data flow from the DSPs to the central memory is controlled by a 68000 based microcontroller. The 
calculation of model position and attitude is performed in an I860 based array processor. The array processor 
resides on a VME bus and is controlled by a SUN computer. 
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Optical Sensing System 



Figure 1 . Block diagram of the optical sensing system showing signal flow. 
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LIGHT DISTRIBUTION 



Figure 2. Drawing showing the light distribution incident on the Charge Coupled Device (CCD) 
detector and the CCD dector output. 



Sensing Unit 


Eight sensing units arc located symmetrically about, and approximately five feet above the levitated model. 
Each sensing unit consists of two linear CCD sensors oriented orthogonally with respect to one another.(figure 
5) The output of each sensing unit is thus an x and y camera location for each projected target image. Each 
sensing unit is mounted to the sensing system support structure. The sensing unit mount allows the 
orientation of the sensing unit to be adjusted in two angular directions as well as the vcrticle direction.(figure 4) 
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Figure 3, Drawing of a sensing unit. Each sensing unit consists of two CCD cameras oriented orthogonally. 
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DIAGRAM OF SENSING UNIT MOUNT 



F.i JO-CV; ICi page blank not filmed 
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Figure 4. Drawing of sensing unit mount. 



Coordinate System Definition 


^ste^/^nr?s? S -n, m f0f LGMSS a ^P° sition measurement system is a right hand cartesian coordinate 
systcm.(figurc 5) The z axis points up. The position of the levitated model is defined by the x y and z 

T r" °I k b ‘? y c fi ’| Cd ref< ; rcnce framc in laboratory (or fixed) reference frame! The orientation 
of Uic model is defined by the Euler angles y, 0, <p (yaw, pitch, roll). The rotation sequence is yaw, pitch and 
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Transformation Equations 


The position and orientation of the model are determinal by solving a set of nonlinear equations which relate the 
measured x and y camera coordinates of each target to the six parameters, x C m, y c m • Zcm> an< ^ 9- Two 

sets of transformation equations arc lislal below. The first set of transformation equations describe the 
transformation from body fixed reference frame coordinates (xfo, yb, y.t>) to laboratory or fixal reference frame 
coordinates (xf, yf, zf) in terms of the coordinates of the origin of the body frame in the laboratory frame (x cm , 
y cm , and /cm) and the orientation of the body frame with respat to the fixed frame, given by the sines and 
cosines of the Euler angles psi (y), theta (0), and phi (<p). The second set of transformation equations describe 
the transformation from laboratory reference frame coordinates for the jlh target (xfj, yfj, zfj) to the ith camera 
(xjj, yjj) in terms of the x, y, z position of the perspative center of each camera in the laboratory reference 
frame, denoted by Xj c , Y, c , and Zj c , and the oricnuition of each camera with respat to the laboratory reference 
frame. In the second set of equations, the mij arc elements ol the rotation matrix and arc functions of the camera 
pointing angles. The terms xpj and ypj arc the principal points lor the ith camera, and l x ,yi the focal length 
of camera i. The positions and orientations of each camera are determined by an independent survey and through 
calibration. 

Substituting the first set of transformation equations below into the second produces a set of nonlinear equations 
which daribes the dependence of the xij and yij positions of the projected images of target diode j in sensing 
unit i in terms of the position and orientation of the model. The position and orientation of the model are given 
by the six parameters x cm , y cm , z cm , y, 0, <p. With eight targets and eight sensing units, 128 equations are 
thus generated. These equations arc linearized and solved using an iterative tahnique (Newton's Method). 

xfj = x cm + xbjfcosCOjcosOF) + ybj(-cos(<p)sin(y) +sin(y)sin( 0 )cos(y))+ 
zbj(sin(tp)sin(y) + cos(cp)sin( 0 )cos(y)) 

yfj = ycm + xbj(cos( 0 )sin(y) + ybj(cos(<p)cos(y) + sin(<p)sin( 0 )sin(y) 
zbj(-sin(<p)cos(y) + cos(<p)sin( 0 )sin(y)) 

zfj = Zcm + xbj(-sin( 0 )) + ybj(sin(<p)cos( 0 » + zbj(cos(y)cos( 0 » 


xij = x p i + fi x mn(xfj - X x jC) + m i2(yfj - Y xi c ) + mi3(zfj - z xi c ) 
m x3l( x fj - X x j c ) + m x 32(yfj - Y xi c ) + m x33( z fj - z x 

yij = ypi + fiyf ~™2l( x fi - X v jC) + m22(yfj - Y yi c ) + ™23( y -fi - z yi c ) 


jin y 3 1 (x [j - X yi c) + m y 32(yij - Y y jC) + m y 33 (z fj - Z y jC)j 
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Targets 


As stated previously, eight infrared light emitting diode targets are embedded in the top surface of the model. 
Each target consists of an LED chip soldered into a small parabolic reflector.(figurc 6) The wavelength of the 
light is 850 nm. The diameter of the reflector is 0.047 inches. Figures 6-11 show a diagram of the cross 
section of a target, a picture of a target compared to a dime, a graph of the spatial distribution of the light output 
from an LED target, a picture of a cylindrical model with the targets embedded in the surface, a picture of the 
circuit board for the target driver electronics, and a diagram showing the timing pulses. The larger diameter 
s own In diagram 6 is a heat sink. Special bonding techniques had to be developed which would provide 
adequate heat conduction when the diodes are driven at high current levels. The materials which are used to bond 
the LED and leads to the reflector are noted in diagram 6. 

The driver electronics and power supply for the targets are all located on the model. The size of the driver 
electronics package is about 1 inch in diameter and 0.25 inches in thickness. Two nickel cadmium 
(rechargeable) batteries are provided. The dimensions of the batteries are 1 inch in diamter and about 1 inch in 
length. These batteries can provide sufficient power to operate for about 2 hours. 

In addition to the target diodes, two additional strobe diodes, located on the bottom surface of the model arc used 
to time the sensor acquisition of each target and detect the flashing of the first of the eight targets ( or the 
beginning of frame). Diagram 11 shows the timing pulses for the prototype model. The timing diodes are 
separated from the target diodes in wavelength. The photodetectors used to detect these diodes, which are located 
beneath the levitated model, are filtered.. 

During tracking of the model, each of the eight target diodes is turned on in sequence for a period of time which 
may range from approximately 1.25 ms to 3.75 ms. (The length of time is determined by resistor values which 
arc set in driver circuitry located in the model.) During this time, the CCD sensors integrate the light falling on 
Uie array. Following integration, the signal is read out and the centroid location is determined. After all eight 
diodes have been illuminated they are turned off during the time the sensing system is calculating the position 
and attitude of the model. The length of time required to calculate the position and attitude of the model is 1 5 

ms. Following the calculation, the computed position and attitude information is sent to the LG MSS controller 
and a new frame is begun. 
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LIGHT EMITTING DIODE TARGET 
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Figure 6. Drawing of a Light Emitting Diode (LED) target cross section. 
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Figure 7. Picture of a target LED next to a dime. 
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Figure 9. Picture of cylindrical model with targets embedded in the surface. 



Figure 10. Picture of LED target driver electronics circuit board. 
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Ises for beginning of frame and acquisition of a single target, a) Pulse 
by strobe diode for timing, c) Pulses generated by strobe 
of eight diodes. 





Sensor Electronics 


Figures 12 - 14 show pictures of the various sensor componenls. Figure 12 shows the front and back view of 
the camera board. Also shown is the thermoelectric cooler circuitry. Some of Ihc components on the camera 
board are being operated at high clock speeds. As a result, without active cooling, these parts will overheat and 
be destroyed. A thermoelectric cooler is provided to actively cool these parts. A thermoelectric cooler was used 
to avoid inducing vibrations which could disturb the optical measurement. Figures 13 and 14 show pictures of 
the analogue to digital converter and digital signal processor boards, respectively. 
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Figure 12. Pictures of camera board (front and back view) and thermoelectric cooler board. 
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Figure 13. Shown in this picture is the analogue-to-digital (A/D) converter board. The A/D is a 10 bit flash 
capable of a maximum data conversion rate of 20 megasamples/second. 





Tcxas^nstrumenK^ a picture of the digital signal processor (DSP) board. The DSP is 

lexas Instruments TMS320C30. The clock speed is 32 mHz which yields an instruction cycle time of 

raUo^atouM 6 ^" 8 ' ^ ^ UmC l ° pcrform ** calculation of centroid, background, and signal-to 


a 

noise 
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Graphs of Output of Prototype Sensor 

figure 1 5 is a graph showing the digital output of a prototype CCD sensor when imaging the cylindrical model 
ol ligure 9. The integration time of the camera was set to integrate over the entire frame lime of all eight 
targets Figure 16 is a graph of a portion of the digital output of a camera when imaging a single LED target 
The light distribution spans approximately 1 1 pixels. ° 
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Figure 16. Graph of digital output of a camera for a single LED target. 
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Conclusions 


Early simulation results were used to determine the level of signal-to-noise which would be needed to achieve 
the accuracy specified for the position measurement system. No model dynamics were taken into account in the 
simulation. Rather, the model was assumed to be stationary over the total 10 ms time in which the target 
locations were being sensed. The signal-to-noisc ratio required to achieve the specified accuracy for the sensing 
system has been demonstrated. Also, it appears that a sample rate in excess of 20 samples/second can be 
achieved. However, further tests are required before it is known how much faster the system can track. 
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Conclusions 


Simulation and test results indicate that the 
accuracy requirements for the sensing system 
will be met 

Software is being streamlined to improve the 
frequency response, should be possible to 
sample > 20 samples/sec 
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APPENDIX 


LIST OF ATTENDEES 
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